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Zur Theorie der Verbrennungsgeschwindigkeit, 
mit einer Anwendung auf die Ozonverbrennung 


Von Karl Bechert i 


Inhaltsübersicht N 


Es wird eine einfache Theorie der Verbrennungsgeschwindigkeit entwickelt. Die Ak- 
tivierungsenergie E, die im Ansatzfür die Zahl Z der zur Reaktion führenden Molekül- 
zusammenstöße vorkommt, wird als verfügbarer Parameter behandelt, der aus einer 
gemessenen Verbrennungsgeschwindigkeit zu bestimmen ist. Es gelingt dadurch, die 
Diffusion summarisch zu berücksichtigen, ohne daß sie explizit in den Rechnungen er- 
scheint. Die Reaktionen werden durch eine Gesamtreaktion (Bruttoreaktion) ersetzt, 
welche die Ausgangsstoffe mit den Endprodukten verknüpft die Zwischenprodukte (die 
„aktiven‘‘ Teilchen) erscheinen nicht explizit in der Rechnung, ihre Wirkung steckt in 
der geeignet gewählten Aktivierungsenergie E. — Aus Impuls- und Massenerhaltung läßt 
sich eine Beziehung gewinnen für die Strömungsgeschwindigkeit u des Gases an beliebiger 
Stelle, in Abhängigkeit von der Verbrennungsgeschwindigkeit u,, der Schallgeschwindig- 
keit c, im Unverbrannten (Index /) und den Verhältnissen 7'/7',, M/M, der Temperaturen 
und Molekulargewichte an der betrachteten Stelle und im Unverbrannten. Die Formel 
ist verschieden, je nachdem, ob es sich um Verbrennung oder Detonation handelt (Gl. (8) 
und (13), $2). Die Verbrennungsgeschwindigkeit u, gegenüber dem Unverbrannten hat 
bei Verbrennungen eine obere Grenze (Gl. (11), § 2), bei Detonationen eine untere Grenze 
(Gl. (15), § 2). Das Verhältnis dieser zwei Grenzgeschwindigkeiten. hängt nur vom Ver- 
hältnis der Temperaturen und Molekulargewichte im Verbrannten (Index r) und Unver- 
brannten ab (Gl. (18b), (20), (21), $2). Aus dem Energiesatz und der Gleichung für die 
Abnahme der Zahl der Brennstoffmoleküle infolge der zur Reaktion führenden Stöße 
gewinnt man eine Gleichung ((38), $ 5), welche 7',, 7,, M,, M, mit der Reaktionswärme 
pro Mol Brennstoff und der Gemischzusammensetzung im Unverbrannten verknüpft. — 
Von jetzt an wird auf Reaktion zwischen zwei gleichen Molekülen spezialisiert. Die Kom- 
bination der 3 Erhaltungssätze und der Gleichung für die Abnahme der Brennstoffmole- 
külzahl gibt die grundlegende Differentialgleichung (47), $5. Sie ist von erster Ordnung, 
vermöge unseres vereinfachten Ansatzes für die Diffusion, und nichtlinear. — Im Unver- 
brannten herrscht nicht thermisches Gleichgewicht, es finden auch da Reaktionen statt, 
aber in so geringer Zahl, daß dort alle physikalischen Größen praktisch konstant sind ($ 6). — 
Die Differentialgleichung (47) wird umgeformt, als zweckmäßige Variable erweisen 
sich =T, genauer: tT —T, = y, (7 = ir) und 
der Brennstoffmoleküle im cm? des Gases bedeutet, und M das Molekulargewicht des 
Gases. Aus (47) wird so die Differentialgleichung (80), $ 7, in der die Verbrennungsge- 
schwindigkeit u, als willkürlicher, physikalisch noch unbestimmter Parameter steht. 
(80) verknüpft die Brennstoffmolekülzahl mit der Temperatur. Der Verlauf der Integral- 
kurven wird in § 8 besprochen. Gesucht ist die Integralkurve, die den Bedingungen im 
Verbrannten und im Unverbrannten genügt. Die gesuchte Kurve (C) muß also durch 
zwei vorgegebene Punkte gehen, da sie Lösung einer Differentialgleichung erster Ord- 
nung ist, hat sie nur eine Integrationskonstante als Parameter, kann also nur dann die 
Bedingung erfüllen, wenn nicht nur die Integrationskonstante geeignet gewählt wird, 
sondern auch die noch verfügbare Verbrennungsgesehwindigkeit u, (§ 8). Es läßt sich 
sogar eine analytische Näherungsformel für die gesuchte Integralkurve rv, (y) angeben 
(Gl. (95), $ 9), und eine sehr einfache Näherungsformel für die Verbrennungsgeschwin- 
digkeit %,, in der nur elementare Rechenschritte vorkommen (Gl. (102), $9). Das hier 
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angewendete Verfahren zur näherungsweisen Bestimmung des Parameters läßt sich grund- 
sätzlich auch bei anderen Differentialgleichungen verwenden. Das Problem sei so be- 
schaffen, daß es nur bei bestimmter Wahl des Parameters eine Integralkurve gibt, welche 
den gestellten Randbedingungen genügt. Man setzt als Näherungslösung eine Funktion 
an, welche die gestellten Bedingungen erfüllt und möglichst weitgehend die Eigenschaften 
der gesuchten Integralkurve hat. Einsetzen dieser Näherungslösung in die Differential- 
gleichung gibt die (näherungsweise) Bestimmung des Parameters. Je nach der Stelle, an 
der die Näherungslösung die Differentialgleichung erfüllen soll, ergeben sich verschiedene 
Werte des Parameters, die aber, wenn die Näherung gut ist, praktisch übereinstimmen 
müssen. Diese Probe auf die Güte der Näherung fällt bei unserem Ansatz sehr günstig 
aus, drei verschiedene Anpassungsforderungen liefern dieselbe Formel für die Verbrennungs- 
geschwindigkeit ((102), $9). Die Verbrennungsgeschwindigkeit hängt ab von der Akti- 
vierungsenergie, der Reaktionswärme, Wärmeleitfähigkeit, der Gemischzusammensetzung 
im Unverbrannten, vom Radius und der Masse der reagierenden Moleküle, und auf dem 
Umweg über die Temperatur im Verbrannten von den spezifischen Wärmen der Gase, 
und natürlich auch von der Temperatur im Unverbrannten. — Durch Dimensionsbetrach- 
tungen läßt sich zeigen ($ 9), daß die dimensionelle Abhängigkeit der Verbrennungsge- 
schwindigkeit von den physikalischen Größen notwendig aus den Erhaltungssätzen für 
Masse, Impuls, Energie und der Gleichung für die Abnahme der Zahl der Brennstoff- 
teilchen folgt [(103d), oder in umgeschriebener Form (103n), § 9], in uj kommt so ent- 
scheidend die Reaktionszahl Z vor. Der Ansatz (28) für Reaktion zwischen gleichen Teil- 
chen, der in § 3 begründet ist, führt zwangsläufig von (103d) auf die dimensionelle Form 
der Gl. (102). Zur Festlegung der bei der Dimensionsbetrachtung noch unbestimmt 
bleibenden Größen (ob r, oder r, zu wählen ist, bleibt fraglich, und ähnliches) ist aber 
Rechnung nötig, die eben in dieser Arbeit ausgeführt ist und (102) liefert. Man kann sogar 
zeigen, daß auch die strenge Berücksichtigung der Diffusion an der dimensionellen Zu- 
sammensetzung der Verbrennungsgeschwindigkeit nichts ändert ($ 9). Damit ist (103d) 
als allgemeingültig erwiesen für stationäre Verbrennung, wenn die Reaktion nur durch 
gleiche Teilchen zustande kommt. — (103n) wird zunächst mit der Formel von Nusselt 
verglichen, die mit stark vereinfachenden Annahmen und dem physikalisch nicht ein- 
wandfreien Begriff der Zündtemperatur abgeleitet worden ist (§ 10). Die Nusseltsche 
Formel hat trotz ihrer anderen Herleitung den Bau der Gl. (103n), die aus ihr gezogenen 
Schlüsse über die Abhängigkeit der Verbrennungsgeschwindigkeit von der Wärmekapazi- 
tät, von der Wärmeleitfähigkeit 2 und der Reaktionszahl Z sind bindend, weil sie gar 
nicht charakteristisch sind für Nusselts Formel, sondern notwendig aus den Erhaltungs- 
sätzen folgen. Der Einwand, der gegen die Abhängigkeit der Verbrennungsgeschwindig- 
keit von A gemacht worden ist, wird zu entkräften versucht durch eine Betrachtung über 
die Wärmeleitfähigkeit von Gemischen (§ 10). Unsere Formel (103n) ist der Nusseltschen 
darin wesentlich überlegen, daß der unzulässige Begriff der Zündtemperatur nicht vor- 
kommt; er erweist sich als unnötig, sobald die Reaktionszahl Z nicht konstant gesetzt 
wird, wie es bei der. Ableitung der Nusseltschen Formel geschieht. — Bei Reaktion 
zwischen gleichen Teilchen steigt die Verbrennungsgeschwindigkeit monoton mit der 
Brennstoffkonzentration (§ 10). Kurzer Hinweis darauf, daß bei Reaktion zwischen un- 
gleichen Teilchen die Verbrennungsgeschwindigkeit mit steigender Brennstoffkonzentra- 
tion über den stöchiometrischen Punkt hinaus wächst, und dann wieder fällt (Zünd- 
grenzen), wie in der nächsten Arbeit gezeigt werden soll. — Von der Temperatur im Un- 
verbrannten ist uw, gemäß (102) verhältnismäßig wenig abhängig, im Einklang mit der 
Erfahrung, und im Gegensatz zur Nusseltschen Formel ($ 10). — Vergleich von (102) 
mit den Messungen an der Ozonverbrennung zeigt gute quantitative Übereinstimmung 
(Tabelle 1, $10); für Z/k und den Ozonmoleküldurchmesser ergeben sich vernünftige 
Werte ($ 10). u, hängt nicht vom Druck ab, wenn man den Ansatz (28) für die Reaktions- 
zahl macht, die Messungen scheinen dies zu bestätigen (§ 10). — Gl. (80) kann man auch 
graphisch integrieren ($ 10). — Der Temperaturverlauf bestimmt sich aus der Differential- 
gleichung (106), die entweder graphisch oder näherungsweise analytisch integriert werden 
kann (GI. (112), $11). Eine Abschätzung der Flammendicke geben die Formeln (115), 
(115a) $11, die Dicke liegt in der Größenordnung von 5 10-3 bis 10-2 cm und hängt 
von 7',, T,, Reaktionswärme, Aktivierungsenergie, Gemischzusammensetzung und Druck 
im Unverbrannten, und von der Verbrennungsgeschwindigkeit ab. Zum Druck im Un- 
verbrannten ist die Flammendicke umgekehrt proportional, wenn wir den. Ansatz (28) 
benutzen; bei anderem Ansatz für die Reaktionszahl Z ergibt sich andere Druckabhängig- 
keit der Flammendicke. Wie für die Verbrennungsgeschwindigkeit läßt sich zeigen, daß 
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ind- die dimensionelle Abhängigkeit der Flammendicke von den physikalischen Größen sich 
be- zwangsläufig aus den Grundgleichungen ergibt; Gl. (115a), $ 11 ist in der dimensionellen 
Iche Zusammensetzung unabhängig von speziellen Annahmen über Reaktionszahl und von 
ree der Schematisierung der Diffusion. 
Hauptergebnisse der Arbeit 
ngs- 1. Stationäre Wanderung einer Störung ist in einem Gas oder Gasgemisch 
oe nicht mit beliebiger Geschwindigkeit möglich. Massen- und Impulserhaltung ver- 
rn. langen vielmehr, daß die Wandergeschwindigkeit einer solchen Störung in einem 
tase, von zwei „erlaubten‘‘ Gebieten liegt, nämlich entweder kleiner ist als die durch 
ach- Gl. (11) $ 2 gegebene Geschwindigkeit, oder größer als die durch Gleichung (15) $ 2 
"io gegebene Geschwindigkeit. Diese Aussage gilt für Störungen irgendwelcher Art, 
toff. bei denen die physikalischen Variablen zu beiden Seiten der Störung verschiedene 
ent- Werte haben, also für Vorgänge mit oder ohne chemische Reaktion. 
wi 2. Die dimensionelle Abhängigkeit der Verbrennungsgeschwindigkeit «, von 
mmt den physikalischen: Daten läßt sich aus Dimensionsbetrachtungen gewinnen 
aber ($9, Gl. 103d und folgende). Der Grundgedanke ist: Man schreibe alle Gleichun- 
g gen hin, aus denen eine genaue Rechnung die gesuchte Größe liefern würde. Dann 
03d) kann man durch Divisionen und Multiplikationen die Dimensionsabhängigkeit 
urch der gesuchten Größe aus den Gleichungen ohne weitere Rechnung finden. Als 
selt dimensionsloser Faktor kann eine Funktion hinzukommen, die nur von den 
= dimensionslosen Größen abhängt, die man durch Division aus den einzelnen 
Bi Gliedern der Gleichungen und aus den Randwerten bilden kann. Das Verfahren 
pazi- ist allgemein. Ich bin erst nach Beendigung der Arbeit darauf aufmerksam ge- 
gar worden, daß derselbe Grundgedanke, aber eine andere Art, ihn auszuwerten, 
ii. schon bei E. Fues, Uber die Grenzleistung von Dimensionsbetrachtungen'), 
as angegeben ist. In der Literatur wird dieses Verfahren, von den Gleichungen 
chen selbst auszugehen, nach meiner Kenntnis nie angewendet, obwohl es héchst ein- 
vor- fach ist; auch Fues erläutert es nur an 2 gelösten Problemen der Hydrodynamik, 
in wendet es aber nicht auf neue an. 
der 3. Die Berücksichtigung der Diffusion ändert am Dimensionszusammenhang 
== nichts (§.9, Gl. 103n), liefert vielmehr nur die Beziehung zwischen Wärmeleit- 
ünd- fähigkeit, Diffusionskonstante, Dichte und spezifischer Wärme (§ 9, Gl. 103m), 
Un- die aus der kinetischen Gastheorie bekannt ist. Die letztere Beziehung ist eine 
on Folge der Erhaltungssätze und von kinetischen Vorstellungen unabhängig schon 
zu ein Ergebnis der Kontinuumsphysik. 
ftige 4. In der Arbeit wird ein einfaches Schema angegeben, nach dem Verbrennungs- 
Ze geschwindigkeiten angenähert analytisch berechnet werden können. In dieser 
tial- Arbeit wird das Schema für die Berechnung der einzigen genau gemessenen Ver- 
rden brennung verwendet, der Ozonverbrennung. Das Schema ist den bisherigen For- 
115), meln für die Berechnung der Verbrennungsgeschwindigkeit überlegen, es ver- 
=. meidet den von Nusselt verwendeten physikalisch unzulässigen Begriff der 
Un- Ziindtemperatur (§ 10). 
(28) — 
EN 1) E. Fues, Z. Physik 101, 665 (1937). Vgl. auch K. Bechert, Z. Naturforsch., 
Sommerfeldheft 1948. 
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5. Für die Näherungsformel der Verbrennungsgeschwindigkeit wird ein Rechen- 
verfahren angewandt, das auch bei anderen Differentialgleichungen brauchbar 


Bi ist, wenn ein Parameter in der Differentialgleichung so bestimmt werden soll, 
Br daß die Randbedingungen erfüllt werden. Die Güte der Näherung kann in ein- 
facher Weise beurteilt werden (§ 9). 
Es soll der Versuch einer einfachen Theorie der Flammenwanderung gemacht 

Bi; werden). An Theorien der Verbrennungsgeschwindigkeit fehlt es nicht, eher schon 
we an genauen Messungen. Die beste MeBreihe scheint die von Lewis und v. Elbe 
Be zu sein®). Eine Kritik der bisherigen Theorien findet man in dem Buch von Jost ®). 
Per Physikalisch am wenigsten angreifbar ist die von Lewis und v. Elbe gegebene 

(ings Rechnung, in der die Annahme gemacht wird, daß die Summe von chemischer und 


a Wärmeenergie über die Brennzone konstant ist >). Ohne vereinfachende Annahmen 
kommt man nicht aus bei der Vielfalt der chemischen und physikalischen Vor- 
gänge, die in der Flamme eine Rolle spielen. Es ist die Frage, ob die Vereinfachun- 
gen, die man der Rechnung zugrunde legt, physikalisch einleuchtend sind, und 
ob sie sich an den Beobachtungen bewähren. Unsere grundlegenden Annahmen 
stehen in $ 2 und 3; wesentlich für unsere Rechnung ist außer der üblichen Annahme 


7% stationärer Verbrennung, daß wir nur mit einer einzigen Reaktion, der Brutto- 
‘a : reaktion, zu rechnen brauchen, und daß wir die Diffusion dadurch berücksichtigen, 
Br; daß wir die Aktivierungsenergie als Parameter betrachten. Es ist der einzige der 


Erfahrung anzupassende Parameter, der in unserer Theorie vorkommt. Daß die 
Übereinstimmung mit der Erfahrung gut ist, zeigen die Ergebnisse in $ 10. 

Die Betrachtungen bis Gl. (21) einschließlich gelten für Verbrennung und 

j Bi: Detonation; von da an beschränken wir uns auf Verbrennungen. Bis zum Ende 

1 von § 3 rechnen wir fiir Verbrennungen mit beliebigen Reaktionsgleichungen, doch 

sind auch noch die späteren Gln. (32) bis (38) unter den dort angegebenen Voraus- 


By setzungen für solche allgemeinen Verbrennungen gültig. Im übrigen beschränken 
ee? wir uns von § 4 an auf Reaktionen zwischen zwei Molekülarten, und von $5 an 
BR. auf Reaktionen zwischen zwei gleichen Molekülarten. Unter den im Text ange- 
oo gebenen, für solche Verbrennungen wohl fast immer zulässigen Vernachlässigungen 
ng gilt alles Folgende allgemein, nur für die numerischen Rechnungen in § 10 und 11 
ae zur Verbrennungsgeschwindigkeit und Flammendicke bei der Ozonverbrennung 
a f sind die Zahlenwerte der Ozonverbrennung natiirlich wesentlich. 
bs oF § 2. Folgerungen aus der Erhaltung von Masse und Impuls = 
ss Es handle sich um eine Flamme, die bei der Verbrennung eines Gasgemisches 
entstanden ist, und die von rechts nach links fortschreitet, so daß links (/) das un- 
a verbrannte Gemisch ist, rechts (r) das verbrannte. Den Vorgang schematisieren 
te wir als eindimensional (kartesische Koordinate x) und stationär, so daß alle physi- 
- 
> 2) Die Arbeit ist in allen wesentlichen Zügen fertig seit August 1944, ein kurzer Bericht 
> wurde am 29. 4. 1948 bei der’Griindungsversammlung der Physikal. Gesellschaft Rhein- 
an land-Pfalz in Mainz vorgetragen. 
Ps 3) B. Lewis u. G. v. Elbe, J. chem. Physics 2, 283 (1934). 
” 4) W. Jost, Explosions- und Verbrennungsvorgänge in Gasen, Springer, Berlin, 
1989, S. 104—119. 


- ee B. Lewis u. G. v. Elbe, J. chem. Physics 2, 537 (1934). Er 
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kalischen Größen nur von x abhängig sind. Reibung vernachlässigen wir, ebenso, 
wenigstens explizit, die Diffusion, der wir durch geeignete Wahl der Aktivierungs- 
energie gerecht werden wollen. Das Verbrannte ist dadurch gekennzeichnet, daß 
in ihm die Gesamtreaktion vollständig abgelaufen ist, welche die Verbrennung ver- 
ursacht. In ihm herrscht daher konstante Temperatur 7',, die Flamme ist da zu 
Ende, wo 7’ = T, geworden ist. Einen Temperaturwert, zunächst an beliebiger 
Stelle, im Unverbrannten nennen wir 7',. Links herrscht nämlich nicht genau 
konstante Temperatur, weil auch bei niedriger Temperatur theoretisch die betrach- 
tete Reaktion noch möglich ist, wenn sie auch mit sehr kleiner Geschwindigkeit 
abläuft. Praktisch ist allerdings auch links die Temperatur in der Regel räumlich 
konstant. Die numerischen Rechnungen in $9, 10,11 beziehen sich auf die Ozon- 
verbrennung, und für diese wird in $6 gezeigt, daß 7’, wirklich räumlich sehr 
gut konstant ist. Aber die räumliche Konstanz von 7’, int für die Überlegungen 
im jetzigen Paragraphen nicht Voraussetzung. 

Wegen der Annahme der Stationarität ist die Verbrennungsgeschwindigkei t 
zeitlich unveränderlich. Für die Rechnung betrachten wir die Flamme als ruhend ; 
dann strömt das Unverbrannte mit der (konstanten) Geschwindigkeit u, in Rich- 
tung positiver x auf die Flamme zu, das Verbrannte strömt mit der Geschwindig- 
keit u,. Die Verbrennungsgeschwindigkeit, bezogen auf ruhendes unverbranntes 
Gemisch, ist dann — u,. 

Die Massenerhaltung gibt, wenn wir die Strömungsgeschwindigkeit mit w, 


die Dichte mit o bezeichen: 
ou = const = 0,4). (1) 


Den Druck nennen wir p, dann lautet die Impulserhaltung: / er, 
(2) 
Fug * 4 é 
Das Molekulargewicht sei M, die Gasgleichung wird: 


RT 
M ist wie die übrigen Größen p, o, T, u von x Fa Aus den 3 Gleichungen 
kann man eine Beziehung gewinnen zwischen «, u, 7, T b st M,. Einsetzen von 


(3) in (2) gibt nämlich: 


RT 
o( M wt) = ‘+ ui) 
Division durch ou; gibt mit den Abkürzungen: 
RT, T M 
wegen (1): 


T' u 
Y wt(=) u, T 1); 
daher zunächst 


X ı 
ziehung, die man für numerische Diskussion und graphische Darstellung von (5) 


und den folgenden Gleichungen und Ungleichungen gut brauchen kann. 


Ersetzt man in (5) y formal durch 3.3 , so geht ni in 2. z a; eine Be- 
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Welches Vorzeichen zu nehmen ist, zeigt die folgende Überlegung: im Grenzfall 
T—T, muß gelten: M— M,„u— u,; da aber in diesem Grenzfall aus der vorigen 


Formel wird: 
4 Uy 1-4 | 1 | 


so gilt fiir y > 1 das untere Vorzeichen in (5), fiir y < 1 das obere. 


An (5) sieht man, daß eine stationäre wandernde Flamme nur dann 
möglich ist, wenn die Wurzel reell ist, d.h. wenn: 


) 
> by 
j 
| M’ (F- 


Man hat jetzt zwei Fälle zu unterscheiden: 


1. 
27" 
> +2) ap Gr 
= 


gr 
M’ 
d.h.u,< ‚wenn wir die Schallgeschwindigkeit c einführen durch c? = 


Es ist stets *) 7’ > 1, M’ < 1, daher > 1, für den Fall 1. wird so y > 1, 


xRT 
M 
Dies ist der Fall gewöhnlicher Verbrennung (kleine Verbrennungsgeschwin- 
digkeit); für ihn gilt also: 

Ur 


(1 (I u< ;y= 


%ı 


(8) 


Gl. (8) gibt die Strömungsgeschwindigkeit in Abhängigkeit von 7M und von den 


c 
Größen u,, —— , 7,/M,. 
x 


/ 


1 
Die Ungleichung (7) muß für jedes 7’, M’ der Verbrennung richtig sein, also 
auch für die Temperatur 7’, des Verbrannten und das zugehörige Molekularge- 


wicht M,. Wir setzen: 
9 


T’’;M’” ist der größte Wert, den 7’/M’ haben kann. Da die rechte Seite von 
(7) mit wachsendem 7’/M’ wächst, co ist die wirkliche untere Grenze für y 
gegeben durch die Ungleichung: 
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"u 6) 7, M sollen sich natürlich auf einen Punkt rechts vom Punkt 7, M, beziehen. 


en 


geschwindigkeit gegenüber dem Unverbrannten nennen wir Up, die größte, bei 
Verbrennungen stationär mögliche Verbrennungsgeschwindigkeit gegenüber dem u 

Unverbrannten Up. Mit wachsendem W wird U ; kleiner, U, größer. Das ,,ver- ‘3 
4 

botene‘‘ Gebiet von Geschwindigkeiten «,, für die keine stationäre Flammen- 3 

wanderung möglich ist, wird also mit wachsendem jy immer breiter. ¥ 
Die ,,Schwadengeschwindigkeit“, bezogen auf ruhendes unverbranntes Gemisch, ist 4 

aa 

Ug = Uy — Uy. (16) 

Im Falle der Verbrennung ist . a 

(1 + y) 47 >0 ) - (17a) 

7) Allgemein ist u ju, firy21. 


Wir erhalten aus (10): 


(FF -1+3 y ar (7 1)) 

Die Verbrennungsgeschwindigkeit gegenüber dem Unverbrannten 

kann bei stationären Verbrennungen nicht höher sein als die aus (11) 

folgende obere Grenze für u,. 
Der zweite gemäß (6) mögliche Fall ist: 


” 2 ” 
2. 
das gibt wegen (6) 


Die rechte Seite ist wegen 7’/M’ > 1 positiv und kleiner als 1, so daß auch 
y < 1 gilt. Dies ist der Fall der Detonation (große Verbrennungsgeschwindig- 


keit), es.ist u, > VE . In (5) gilt hier das obere Vorzeichen gu. 


Die rechte Seite von (12) nimmt ab mit zunehmendem 7’/M’, die wirkliche Grenze 
füry, diesmals die obere, ist also wieder durch 7',, M „gegeben. Es muß demnach 
sein: 


rT 
das bedeutet: 
u, a — - » (15) 
27” / Tr’ 47” 
‚1-2 ar 1)) 


Die Verbrennungsgeschwindigkeit gegeniiber dem Unverbrannten 
ist beistationären Detonationen größer als die aus (15) folgende untere 
Grenze fiir u,. Die kleinste, bei Detonationen stationär mögliche Verbrennungs- 
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6) 
) 
\- 
af 
L 
) 
) 
l 
‘ 


der oberen Grenze U, von u, entspricht der Wert 4 


> 


i, =U, 1-9 41-9); Uy = 

(17b) 
Im Fall der Detonation ist uy < 0, und RAR: 

u 
us -%(1->-] (170) 


der unteren Grenze U, von «, entspricht 


a 
M’” fi 


Tr’ pr 


€ 
_ 
41 ] M” —2y) mM” (Z ~1) 


_ Welchen Wert u, bei gegebenen physikalischen Bedingungen wirklich hat, 
darüber sagen die bisher gegebenen Gleichungen nichts aus, sie grenzen nur den 
Bereich der stationären Flammenwanderung ein, also den Wertebereich der u,. 

[73 
Wig fest gegeben. Zu diesem YW 
schwindigkeiten U,, U, stationärer Flammenwanderung bestimmen, welche als Höchst- 
wert bei Verbrennung, als Mindestwert bei Detonation mit Massen- und Impulserhaltung 
verträglich sind. 
Wir führen die Abkürzung ein: 


Wir denken uns 


lassen sich aus (17b, d) die Ge- 


und bekommen aus (17b, d) 
_1+8 
Ur, 45 


Wegen & > 0 wächst die Funktion (18b) monoton mit wachsendem £. Das Verhältnis der 


zwei Grenzgeschwindigkeiten U, und U, hängt nach (18a, b) nur vom Verhältnis ge: z 
ı 
der Temperaturen und Molekulargewichte im Verbrannten und Unverbrannten ab. 
N 
In der Regel ist 7”’/M’’ deutlich größer als 1, mindestens etwa gleich 3, so daß 


man die Formel (11) vereinfachen kann zu: 


hy U v3 U, | + ; (fiir Verbrennungen). (19a) 
5: 


Zu gehört die Schwadengeschwindigkeit 
T M = 
[77 3 [73 ce ] [73 
9 ~ r ‘ 
Us U V 4 | Xp (1 2 >) (1 Ib) 


Gl. (15) vereinfacht sich dann zu: 


„2a M”\_ 2¢, M"\ (für Deto- 
4 Un® (1- m) (19) 


V x, V* nationen). 
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Zu Up gehört Be 
1 M" 


J Cr — 


fist in der Regel wesentlich größer als Eins, mindestens etwa gleich 6, so daß (18b) 
umgeschrieben werden kann in: 


(7) 16 1) (20) 


der kleinste stationär mögliche Wert Up der Detonationsgeschwindig- 
ww yi- 77 größeralsdiegrößtestationär 
mögliche Verbrennungsgeschwindigkeit U,; 7” ist das Verhältnis 7',/7, 
der Temperaturen im Verbrannten und im Unverbrannten, es kann (bis auf einige% 
genau) aus dem Energiesatz berechnet werden, s. $5 und 7; M’ = M,/M, ist das 
Verhältnis der Molekulargewichte im Verbrannten und im Unverbrannten, es 
liegt meist nahe an Eins. Hinreichend genau ist für 7’/M” 3: 


keit ist alsoumden Faktor 4 


U T, M, 


Im Folgenden soll nur von Verbrennungen die Rede sein, nicht von 
Detonationen. Wir behandeln die Formel (8) weiter. Da y gewöhnlich sehr viel 
größer als 1 ist (für die Ozonverbrennung von der Ordnung > 10%), kann man die 
Wurzel in (8) entwickeln und erhält: 


T’ An 
T’ 


‘ 


Meist, auch für die Ozonverbrennung, kann man sich mit dem ersten Glied der 
Entwicklung begnügen (bei der Ozonverbrennung wächst 7’/M’ jedenfalls nicht 
über 20), dann gilt: 


uwru,-T/M. 22a) 


Die spezifische Wärme bei konstantem Volumen nennen wir c,, die bei kon- 
stantem Druck c,, die Wärmeleitfähigkeit A, die bei der einzelnen Brutto-Molekül- 
reaktion pro zerfallendes Molekül freiwerdende Reaktionswärme heiße g’, die Zahl 
der pro cm® und sec stattfindenden Bruttoreaktionen sei Z. Mit Bruttoreaktion 
meinen wir die Reaktion ohne Berücksichtigung der Zwischenstufen, also den reak- 
tionsmäßigen Zusammenhang zwischen den anfangs vorhandenen Molekülen 
und den Endprodukten. Bei der Ozonverbrennung z. B. schematisieren wir die 
Reaktion als: 

20,—30, (23) 


ohne Berücksichtigung der dazwischen gebildeten aktiven O-Teilchen. q’ ist hier 
die pro O,-Molekül freiwerdende Reaktionswärme der Reaktion (23). 


In Wirklichkeit wird bei geringer Temperatur der Zerfall des O, in O, und O über- 
wiegen, der proportional der O,-Zahl vor sich geht. Bei höherer Temperatur, im wesent- 
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lichen Teil des Temperaturanstiegs, tritt die Reaktion zwischen O, und O hervor. Dieser 
Zerfall des O, ist proportional der Zahl der O,- und der Zahl der O-Teilchen. Unserer Rech- 
nung legen wir die Gesamtreaktion (23) zugrunde und bekommen. für die Zahl der Reak- 
tionen Proportionalität mit dem Quadrat der O,-Zahl (s. $4) und mit e”?/*?, Der Ex- 
ponent liegt für 7 ~ 7, in der Größenordnung von —30, für 7’ ~ T, in der Größenord- 
nung von —5. Die Zahl der Reaktionen steigt zunächst mit zunehmendem T steil an 
und fällt dann wieder steil ab, weil wegen der stark zunehmenden Reaktionswahrschein- 
lichkeit die Zahl der O,- Moleküle stark abnimmt. In diesem Gebiet ist die Zahl der 05- 
Teilchen größenordnungsmäßig gleich der Zahl der O-Teilchen, so daß unser vereinfachtes 
Schema hier dasselbe leistet wie eine Rechnung mit dem komplizierteren Reaktions- 
schema 

0,20 + O; O + O, > 2 0,. (23a) 
Andererseits bestimmt der Faktor e~*/*? das Verhalten der Reaktionszahl bei geringer 
Temperatur so sehr — er verursacht den steilen Anstieg von sehr kleinen Werten an mit 
zunehmendem 7 —, daß es in diesem Gebiet nicht so sehr darauf ankommt, ob vor der 
Exponentialfunktion das Quadrat der O,-Zahl oder eine kompliziertere Funktion steht. 


Eine weitere wesentliche Vereinfachung ist, daß wir die Diffusion nicht explizit 
in Rechnung setzen, ihr vielmehr durch geeignete Wahl der Aktivierungsenergie 
genügen wollen. 


Die Diffusion hat mehrere Wirkungen. Erstens: Aus der heißen Zone wandern 
Teilchen, vor allem auch aktive Teilchen — bei der Ozonverbrennung O-Teilchen, bei der 
Kohlenwasserstoffverbrennung Radikale und freie Atome — in das unverbrannte Gemisch, 
und verursachen dort weitere Reaktionen. Dies läuft formal darauf hinaus, daß die 
durchschnittliche Aktivierungsenergie pro Molekül erniedrigt erscheint gegenüber einer 
Rechnung, in der die Diffusion explizit berücksichtigt wird. 

Zweitens: Aus dem Unverbrannten diffundieren Teilchen in die Brennzone und er- 
höhen die Zahl der Reaktionen durch Vermehrung der Zahl der reagierenden Teilchen. 
Dies vergrößert den Faktor, der in der Reaktionszahl bei der Exponentialfunktion e”*!*? 
steht, und wirkt also im Sinn unseres Ansatzes, bei dem wir Proportionalität mit dem 
Quadrat der Brennstoffmolekülzahl annehmen. 

Drittens: Durch den ersten Vorgang wird die Temperatur im Unverbrannten er- 
höht, durch den zweiten Vorgang die Temperatur in der Flamme wohl nicht wesentlich 
geändert, denn die Teilchen aus dem Unverbrannten wirken meist „abkühlend“, solange 
sie nicht reagiert haben. Die Temperaturen 7',, 7, werden davon nicht berührt, "denn T, 
ist die Frischgastemperatur, und 7, ist durch den Energiesatz bestimmt, unabhängig 
von E ($5), im Energiesatz kommt es auf die Reaktionswärme an, nicht auf die Akti- 
vierungswärme. Der Temperaturanstieg von 7, bis 7; wird also durch die eben beschrie- 
benen Vorgänge auf ein räumlich breiteres Gebiet verteilt, als man ohne Berücksichtigung 
der Diffusion annehmen würde. Da der Temperaturanstieg wesentlich durch e~7/*? be- 
stimmt wird ($ 11), so bedeutet die Erniedrigung von E zugleich eine Verbreiterung des 
Gebietes ansteigender Temperatur. 


Es ist also zu erwarten, daß unser Rechenschema trotz der starken Verein- 
fachung die wesentlichen Eigenschaften der Verbrennung qualitativ und einiger- 
maßen auch quantitativ richtig wiedergibt. Die Ergebnisse der Rechnung bestäti- 
gen diese Erwartung, wie der Vergleich mit den Meßwerten der Verbrennungs- 
geschwindigkeit der Ozonverbrennung in § 10 zeigt. Bei diesem Vergleich ent- 
nehmen wir Z aus einem Meßwert. Bei komplizierterem Reaktionsschema hätte 
man mehrere Aktivierungsenergien zu unterscheiden, wollte man genau rechnen. 
Unser vereinfachtes Rechenschema, das alle Reaktionsschritte in einem zusammen- 
faßt, kennt nur einen Wert von E. Gemessen ist in den meisten Fällen auch nur 
ein Wert, eben die Aktivierungsenergie, welche die Reaktionsgeschwindigkeit der 
Gesamtverbrennung bestimmt, ganz im Sinne unseres Schemas. Dieser gemeinsame 
Wert ist übrigens in der Regel so ungenau bekannt, daß die dem ersten Vorgang 
entsprechende formale Erniedrigung von E nicht sinnvoll möglich ist. Dann: ist 
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man also darauf angewiesen, E aus einem Wert der Verbrennungsgeschwindig- 
keit zu entnehmen — wenn diese gemessen ist. Ist sie nicht gemessen, so wird man 
den gemessenen Wert von E unverändert übernehmen, und damit die Verbrennungs- 
geschwindigkeit berechnen. Ist unsere Theorie vernünftig, so muß dann die be- 
rechnete Verbrennungsgeschwindigkeit in der Größenordnung der zu erwartenden, 
aus anderen Daten geschätzten liegen. In einer folgenden Arbeit wird das hier 
gegebene Schema (mit einigen durch das andere Problem gebotenen Änderungen) 
auf die Kohlenwasserstoffverbrennung angewendet werden. 


In Differentialform sagt der Energiesatz, wenn wir die Reibung vernach- 
lassigen : 
d d /,aT 
= 2/9) 4. 2 
(ox c,dT +u 2) 4 pu) = 22) +rq’:Z. (24) 


r ist die Zahl der Brennstoffmolekiile — bei der Ozonverbrennung sind es die O,- 
Moleküle —, die pro Reaktion verschwinden. Durch Integration von (24) erhält 
man mit Benützung von (1), (3): 


T 
ut—uj RT RT\ ‚ar /,aT 


Tı 


-¢, hängt hier außer direkt von 7’, noch auf dem Umweg über die Gemischzusammen- 


setzung von 7’ ab. 


$ 4. Die Zahl Z der Reaktionen 


Die nächste Aufgabe ist, Z zu bestimmen. Dazu beschränken wir uns auf den 
Fall einer Reaktion zwischen zwei Molekülarten « und ß, die Molekülzahlen im cm? 
seien ”,, %g, die Molekülmassen m,, mg. Die resultierende Masse nennen wir Mzeı: 


‘ aie 
Mrei = —, 4 
rel = m. + m, (26) 


Die Molekülradien seien a,, ag. Für die Reaktion kommt es auf die Übertragung 
des Teiles ¢,, der relativen kinetischen Energie an, der zur Bewegung in der Ver- 
bindungslinie der zwei Molekülmittelpunkte gehört. Wir nehmen an, daß nur die 
Stöße zur Reaktion führen, bei denen ¢,, über der Aktivierungsenergie E liegt. 
Eine ee Rechnung gibt fiir ungleiche Molekiilarten a, p: pas 


Z = 2 (ax + ap)’ MB | (27) 


Für die Ozonverbrennung hat man nach dem vereinfachten Reaktionsschema (23) 
Reaktion zwischen gleichen Molekülarten. Die Teilchenzahl im cm# bezeichnen 
wir mit n, die Molekülmasse mit m, für gleiche Moleküle gilt dann statt (27): 
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§ 5. Reaktionen zwischen zwei gleichen Teilchenarten. Ableitung der grund- 
legenden Differentialgleichung (47) 


Es soll sich jetzt nur mehr um Reaktionen zwischen zwei gleichen Molekül- 
arten handeln. Die Reaktion verlaufe nach dem Schema 


| 


ove 


dabei bezeichnet A die Atomart, A,, Ay in der üblichen Weise das Molekül,! ı und 
es muß gelten: 


Br (30) 


Die Zahl der Brennstoffmoleküle A, sei n, im cm?, die Zahl der Moleküle A, sein, 
im cm}, 

Wir betrachten eine Schicht der Dicke dx in der Flamme. Die pro sec und cm? 
links einstrémenden A,--Teilchen kommen nicht alle rechts wieder aus der Schicht 
heraus, ein Teil verschwindet in ihr durch die Reaktion. Da bei jeder Reaktion 
2 A,--Teilchen verschwinden, so gilt für ein Gebiet von der Dicke dx und von 1 cm? 
Fläche, die senkrecht zur x-Achse steht: 


(ny = + 2Z de]; (31) 


dabei ist Z die Zahl der Reaktioneır (29) pro cm? und sec, sie ist durch (28) gegeben, 
wenn dort n durch n,, m durch m, ersetzt wird. Aus (31) folgt: an at) ‘ 


d (n, u) 
de 


Das Integral rechts in (25) läßt sich also ausführen, so daß der Energiesatz wird 
(es ist 7 = 2): 


— q’ (n U— Ny, 
(32) 


Dabei ist n,, der Wert von n, fir 7= 7, 


Wendet man die Gleichung (32) auf 7’ = 7’, an, so läßt sich aus ihr eine Glei- 
chung gewinnen, aus der formal wu; in Abhängigkeit von M,, M,, T,, 7’, bestimmt 


werden kann. Für 7’ = 7’, ist nämlich 4 — 0, weil im Verbrannten nach Voraus- 


setzung konstante Temperatur herrscht, außerdem ist ,, = 0, weil nach Voraus- 
setzung die Reaktion vollständig in der Flamme abläuft. Schließlich ist auch das 


Glied mit z fiir 7’= T, zu vernachlässigen, weil der Temperaturanstieg gegen 


die Flamme zu im Unverbrannten sehr langsam ist; man vergleiche dazu die Be- 


für die Ozonverbrennung in § 6. So erhält man zunächst: 


dT 
rechnung von 
dx/T=T; 
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Dieselbe Beziehung gilt allgemein für beliebige Verbrennungen), 
wenn auf der rechten Seite eingesetzt wird: Verbrennungswärme pro Brennstoff- 
molekül, multipliziert mit der Zahl n,, der pro cm? im Unverbrannten vorhandenen 
Brennstoffmoleküle. g’n,,/o, läßt sich umformen; es ist 


wobei p,, den Partialdruck des Brennstoffs im Unverbrannten bedeutet. Das gibt: 


q = pM, ’ 


L ist die Loschmidtsche Zahl. Wir fiihren die Bezeichnugg ein: % ~7 


Lg =@; (85) 
Q’ ist die pro Mol Brennstoff bei der Bruttoreaktion freiwerdende Reaktions- 
wärme. 2,,/?, ist gleich dem Verhältnis der Brennstoffmolzahl zur gesamten Mol- 


zahl im Unverbrannten, genauer: an der Stelle 7’ = 7',; u, kann mit Hilfe von (8) 
durch u,, 7,/M,, T,/M, ausgedrückt werden, das gibt: 


+ 4 [2 (1 + y)? M” Mm’ 
R 1)= Q Pu, 
£: M, M” M, Pi 4 


wobei : 


Dies ist die angekündigte Gleichung zur Berechnung von u,. in A 


Wir beschränken uns auf den häufigsten Fall, auf y > 1, und behalten dement- 
sprechend in (22) nur das erste Glied bei, so kommt: 


Tr 


/ (a) —1)+ M, 1)= M, en 


Danach wäre u, in einfacher Weise aus 7',, 7, M,, M,, Q’, Pı,/P, zu berechnen. 
Es ist aber zu bedenken, daß c, für jede Temperatur zwischen 7’, und 7’, bekannt 
sein müßte, und dazu müßte man die Zusammensetzung des brennenden Gemisches 
für alle diese Temperaturen kennen. Das sieht nicht nach ernsthafter Schwierig- 
keit aus, weil bei den technisch wichtigen Verbrennungen von Luftgemischen der 
Luftanteil bei weitem überwiegt, so daß c, praktisch den Verlauf der spezifischen 
Wärme der Luft hat — bei der Ozonverbrennung gilt dies natürlich nicht. Die 
wesentliche Schwierigkeit für die Anwendung von (37) zur Berechnung der Ver- 
brennungsgeschwindigkeit besteht darin, daß das Glied mit w; sehr viel kleiner ist 
als die übrigen links. Bei der Ozonverbrennung hat man in CGS-Einheiten: 


1 8) Bei Luftüberschuß. 
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Mittelwert ©, von c, ~ 5,10%; 7, ~ 1500° abs; 7’, = 300° abs; M, ~ 40; 
M, = 32, u, 150. Das gibt: 


Tr 
ir) 1) + (FF 1) 5-10. | 
"Zur Berechnung von u, aus (36) oder (37) müßten also die Größen 7',, T,,... sehr 
genau bekannt sein, viel genauer, als sie empirisch bekannt sein können. Deshalb 
haben wir die Beziehung (36) oder (37) formal genannt. 
Die Gleichungen (36), (37) vereinfachen sich nach der eben enone Ab- 
 schätzung zu: 


Tr 173 
RT,/T Q'p 
c,dT 1)= 
/ 


£ Auch (38) gilt, ebenso wie (33), für beliebige Verbrennungen °), dann 
x An Q’ die Reaktionswärme pro Mol Brennstoff, p,, den Partialdruck des 
ete _ Brennstoffs im Unverbrannten, M, und p, Molekulargewicht und Gesamtdruck im 
Unverbrannten. 

Aus Gleichung (38) kann man die Reaktionswärme Q’ pro Mol Brennstoff be- 
rechnen, oder die ursprüngliche Gemischzusammensetzung ,,/p,, oder 7’, oder 
€ „, wenn die anderen Größen in der Gleichung bekannt sind. Die Verbrennungs- 
geschwindigkeit kommt in (38) nicht mehr vor, wir werden sie aus anderen Be- 
ziehungen gewinnen. 


Zur näherungsweisen Berechnung von 7’, setzen wir: 
= ha (39) 


= f die Z w t. Wir ki 
fi (-f) R 


Tr Ty 
Ty Ti 
mit 


Af, ist ein Mittelwert der Zunahme an F RE pro Molekül infolge der Tempezatu- 
erhöhung von 7’, auf 7’,. 
Die linke Seite von (38) wird so: 


1 
Tr RT,ıT RT, fit2+4f 
R 4 
wegen (fi = cy, wird daraus weiter: 
Mm 2 
T 7 
+ 1) —1)[14 2 Gr ; (43) 
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der letzte Ausdruck in der eckigen Klammer ist meist kleiner als 0,05 und kann dann 


weggelassen werden, dagegen liegt i ”, in der Größenordnung von 0,3. 
ı 
Setzt man für einen Mittelwert von c,: 


A T, Pi Cor Mı 


die gesuchte Gleichung für 7',. (45) benutzen wir später zur Vereinfachung. 

Aus (31a) und (32) bekommt man eine Differentialgleichung für n,u in Ab- 
hängigkeit von 7’. Aus der Abschätzung vor GI. (38) folgt, daß für Verbrennungs- 
geschwindigkeiten wu, der Größenordnung weniger Meter in der Sekunde das Glied 
mit u’ — u; in (32) für alle interessierenden 7’ gegenüber den anderen zwei Gliedern 
in der Klammer gestrichen werden kann; erst für Temperaturen T’x10-T, 
würde es die Größenordnung der anderen zwei Glieder erreichen. So kommt aus 
(32): 


dx = - (46) 
RT RTı , dT rer 
0, Uy, (J c, aT + ) (nmu—nıu) + (7 
und aus (31a) und (28) we 


Tı 


Diese Differentialgleichung werden wir noch umformen, so daB sie eine Be- 
ziehung zwischen n, und 7 allein wird, in der «, als noch unbekannter Parameter 
vorkommt. Durch die Randbedingungen: n, = n,, im Unverbrannten, n, = n,, 
im Verbrannten, ist «, festgelegt. Die Aufgabe ist, «, in Abhängigkeit von den 
anderen physikalischen Größen zu bestimmen. 


$ 6. Berechnung von für 


Die Gleichung (31a) zeigt, daß in der Gegend von 7 = 7, kein räumlich unveränder- 
licher Zustand herrschen kann; es kann nicht on und = gleichzeitig Null sein, weil Z 
nach (28) für 7’ = T, nicht Null ist. Aus (24) folgt mit Hilfe der Gasgleichung: 


dT du R/..aT „dam, 


das zeigt weiter, daß dann auch nicht dT/dx, d?T/dx?, du/dx, dM/dx gleichzeitig Nuil 
sein können. Für 7’ = 7’, herrscht eben kein thermisches Gleichgewicht, es finden viel- 
mehr auch bei dieser niedrigen Temperatur Reaktionen O, — O, statt, wenn auch im Ver- 
gleich zur Zahl der Zusammenstöße verschwindend wenige. Die umgekehrte Reaktion, 
bei der aus 30,-Molekülen 2 O,-Moleküle werden, können wir vernachlässigen, und haben 
es in unseren Gleichungen bereits getan. Sie würde die Zahl der O,-Moleküle geben, welche 
bei der Temperatur 7, mit der vorgegebenen O,-Menge im thermischen Gleichgewicht 
sind. Diese Zahl ist aber verschwindend klein gegenüber der O,-Menge, die im unver- 
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brannten Gemisch als vorgegeben vorhanden ist. Die des unver- 
brannten Gemisches ist eben keineswegs die Zusammensetzung, welche dem Massen- 
wirkungsgesetz entsprechen würde. Für 7 = 7, können wir ebenfalls .die Reaktion 
O, — O, vernachlässigen. Die Zahl der O, -Moleküle, welche bei der hohen Temperatur 7, 
mit O, „Mole külen im thermischen G leichgewicht sind, ist zwar wesentlich größer als die 
entspreche :nde Gleichgewichtsanzahl bei der Temperatur T,, aber immer noch verschwin- 
dend klein gegenüber den im unverbrannten Gemisch vorhandenen O,-Mengen. 
Demnach haben wir für 7’ = T,, die Bedingung n, (7,) = n,, = 0 zu erfüllen, sie be- 


deutet Z = 0, (7) = 0, und wegen (31a) weiter (=) = 0. Nach (48) ist dort auch die 


& 
Bedingung räumlich konstanter Temperatur 7 = 7’, erfüllbar. Im Gebiet rechts von der 
Flamme, im Verbrannten, herrscht also der gleiche physikalische Zustand. In der Flamme 
ändern sich die physikalischen Größen räumlich natürlich sehr schnell, links von der 
Flamme, im Unverbrannten, ändern sie sich auch, wie wir eben gesehen haben. Selbst- 
verständlich sind aber bei niedrigem 7, die physikalischen Größen im Unverbrannten 
räumlich nur sehr langsam veränderlich, was wir jetzt auch direkt nachrechnen wollen. 

Einen deutlichen Hinweis darauf gibt schon die Form von Z. Die Exponentialfunk- 
tion e #/*7) ändert sich gewaltig zwischen 7’, und 7',; E/k hat die Größenordnung 10% ° abs, 
so daß sich der Exponent zwischen 7, ~ 300° abs und 7, ~ 2000° abs von —30 bis —5 
ändert. Die Zahl der Reaktionen ist also für 7, ungeheuer viel kleiner als für die Tem- 
peraturen in der Flamme. Daraus folgt, daß auch die anderen physikalischen Größen sehr 
viel langsamer im Unverbrannten veränderlich sein müssen als in der Flamme. 

Zur Auswertung von: (48) brauchen wir dM/dx und du/dx. Es ist 


n,M, +n, M, (M, — M,) | 


Ny + Ng p 


M= 
und wegen (2), (1) 


außerdem nach (22a): 2 
u _TM:, 


<7, mM’ 


u/u, steigt also nicht über die PEN von 10, andererseits ist w,/c, für die Ozon- 
verbrennung von der Größenordnung von 102, so daß wir schreiben können: 


k T (M, —M,) 


(51) 


Mx M, 


“ Diese Formel werden wir öfters brauchen. 
Bei der Differentiation von M nach x gehen wir vom strengen Ausdruck für M aus. 
Wir führen die Abkürzung ein: 


; 


‘ 


es ist x< 1. Aus (49), (50) bekommen wir: 


dM _ mk T (My, —M,)fdn, 1 1dT «a du 
dn, du 
Zwischen den Differentialquotienten besteht nach (31a) ein Zusammenhang. 
Es ist 
an, „ du 
dz" "dz = 
— drücken wir durch — aus und erhalten für die eckige Klammer in (53): 


dx dx 


weg 
folg 
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außerdem ist 
du u,M, dT 1 T dM ee 


die nach (22) noch zu beachtenden Glieder höherer Ordnung können vernachläßigt werden, 
wie sich nach Ausführung der Differentiation zeigt. 

Aus (53) wird wegen (22a): 
dM n,k(M,— M,) T,M 1dT uM, ı =) 


(54) 


dx = Pit M, 


Ta" "m \dzM Mir 


die zwei Glieder mit d7T/dx heben sich wegen (22a) auf, es bleibt (vgl. (51)): 


dx p, M dx M pıu, M, 
wegen M, > M, ist dM/dx immer <0. Jetzt können wir auch du/dx angeben; aus (54) 
folgt: 
du ‚kTT,(M—M), 22) (56) 


a nimmt also nach rechts zu und ist daher immer > 0. Aus (31a) folgt dann, daß dn,/dx 
immer < 0 ist, was ja zu erwarten war. 


Die Wärmeleitfähigkeit 2 ist für Gase proportional ) T; wir schreiben 
A=A)T. (57) 
dM/dx, du/dx setzen wir in (48) ein und erhalten: 


(« dT T Mj ‚ETT(M—-M)M,, z) RTRTM—M),, z) 


j 


dx M? \dz | M, p, u, pu, M,M, 
(58) 


Die von «+ du/dx, also von der Änderung der kinetischen Energie herrührenden Glieder, 
sind klein gegen die anderen in der eckigen Klammer, es ist nämlich 


TiM*c, TiM-8R T, 8M 
und: 


Pas 
uj T Mj kT T,(M,— M,) Pp, u, M, M, „MT 
T; M M, p, u, “RT KET,(M,—M,) * MT, 
Es bleibt, wenn wir noch (51) beachten: Pret 
dT 2Z RT,/M ) A, aT? 
) — + —» — — = — — ri — +2 . 59 


Wir bilden diese Gleichung jetzt für 7’ = 7', und kennzeichnen dies durch Anhängen 
des Index J an die veränderlichen Größen; die Differentialquotienten dT/dx, d’T/dx? 
schreiben wir abgekürzt 7',, T,,. So wird aus (59): 


n,,M, \M, 
RT, /M M 
=k 7, (5-1); 
n,,M, \M, M, 
das ist von der Größenordnung k T,.für die Ozonverbrennung gleich $ k 7, und deshalb 


klainar s . . . . 1: . . 
viel kle iner als q’, weil die mittlere kinetische Energie pro Molekül natürlich viel kleiner 
ist als die Reaktionswärme pro Molekül. Man kann also (60) vereinfachen zu: 


es ist wegen (51): 


01%, ¢,, : 4, + (61) 
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Eine physikalisch vernünftige Schematisierung der Verhältnisse im RE 


ist die Annahme!P): 


Peet = 0, (62) 


sie bedeutet linearen Temperaturanstieg bei 7 = T,. Wir werden zeigen, daß 7',, sehr 
klein ausfällt und daß 7,,, T,im Gebiet 7 < 7, klein bleiben. Wir schätzen die ( irößen- 
_ ordnungen, es ist in CGS-Einheiten: 
4-8,3- 107 
0,10; 0,18 — 1,5-10?; 300; 
Pr 
(m* + ; 
er 


einführen; es ist das Mischungsverhältnis im Unverbrannten, bei der Ozonverbrennung 
_ das Verhältnis der Molzahlen 0,/O,, also: 
106 
"US 3.7,4- 101.3: 108 
Wir setzen E/k x 1,1- 10%, vgl. $9. Das gibt: 


4 


3,1+-1,4- 10716 -3- 10? 


rs 1,2 10"; m, = 48 1,6 10724; 2a, 3: 10-8. 


Daraus ergibt sich: 


2Z,q' 10-16 - 1,22 - 1038 |/ — + & 4,3; 
1q 07 1 1 48° 1,6- 10% 35-10 4,3; 
Da für Z, > 0 auch 7,,,—0 gehen muß, — wenn keine Reaktion stattfindet, ist die 
_ Temperatur räumlich konstant —, so folgt aus den Größenordnungen der vorigen Zahlen: 
— 64 
al 0, Uy Cy ( ) 


Das Glied mit 7, in (61) ist etwa das 6- 10-tHfoche des Gliedes mit 7’, ,, also ganz zu 
vernachläßigen, T,ı hat die Größenordnung 2,5 10%. 

Jetzt kann man auch die anderen Differentialquotienten berchnen und bekommt aus 
(55): M,, —8: 10°, aus (56): u,, 1,7 10-8, aus (31a):  —1,1+ 10". Wir bilden 
die logarithmischen Differentialquotienten: : 


Mei _9. 49-10, _9.10%; 8-10; 


und aus (28): 2,8 10-7. Alle in den Gleichungen vorkommenden Größen ändern 


Me (a sich also sehr langsam bei T = T, und, wie man aus (59) findet, auch im Gebiet 7’ < T.. 


Wir haben praktisch räumlich unveränderlichen physikalischen Zustand fir T<T.. 


ee § 7. Umformung der Differentialgleichung (47) ER 2 
Wir betrachten zunächst den Ausdruck 


10) Gleichbedeutend ist damit die Annahme, daß man das Glied mit T,,,, gegenüber 


dem Ausdruck links in (61) streichen kann, was durch die untenstehenden Gehlenwerte 
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der im Nenner der rechten Seite von (47) vorkommt. Mit den Abkürzungen (vgl. (39) bis 
(44)): 


(66) 


= 


Af — 3 
Af (T — J (f—f)adT 
wird aus (65) en 
f + - m)* (7 —T;) ty, (T—T). 


Als neue Veränderliche führen wir in (47) ein: Far. i. 


und schreiben (47) in eine Differentialgleichung für N, (t) um. So wird wegen (22a): 


N, 
der Nenner auf der rechten Seite von (47) wird so: #152 2 e 
> wie 
ah 4 ge 
das gibt mit den Abkürzungen: Au 
(A T 1 pr E 4 
1 (2:25, +q ) 2 (69) 


für den Nenner der rechten Seite von (47): 
D, + N. 
q (= D,) . 


Für den Zähler braucht man n,, also wegen (68) M/M, in Abhängigkeit von N, und rt. 
Es ist nach (49): 


1325-1), 


M, 1, 
\M, A! 
für p können wir nach (50) p, setzen, für M,/M, führen wir die Abkürzung ein: 
M, 
Ma’ 
für die Ozonverbrennung ist dies gleich ®/,. Es ist: 
1+m*’ 


wegen (63). Außerdem ist: 


Pr mkT NME 


a-it” 
1+ 
daher: 
M, (m*+u)(y_ Ns 1) (m* + Dy’ 
m* +1 rm m* + u 
die linke Seite von (47) wird: 
_ ur d(N,/r) kr? kr 


(62) 
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wir beachten noch (57), setzen zur Abkürzung i 3 
und erhalten aus (47): x 
dN, N, A(m®*4+1? © 
(74) 


(m* + u)? u; Ni E(u— 1) (m* + 1)? 
— Dy t + Dy) 


Aus dem Energiesatz (32) kommt für 7’ = 7',, wenn wir das Glied mit der kinetischen 
Energie weglassen, was wir vor (38) begründeten (vgl. auch (67)): ” 
(T, — T,) —(A T +q' N; (75) 
das bedeutet wegen (69): es 


D,2,;-= D,. (76) 
Wir setzen daher: 
N, = =% + y. (77) 
Andererseits ist nach der Abschätzung im vorigen Paragraphen: _ 
(T,-T)=gq M,c,,(T,—T)) ~ m* +1 (77a) 
daraus folgt: 
auBerdem ist: 
N * 
P= + me) = =», (78a) 
ı 
D, E (m* + 1) 1 >. as 
So’kommt an Stelle von (74): 
dy A (5 + eu 
dy Tr +Yy u \m* + u ) (t, + y) (v—y) 
(tr + (m* + u) (1 =) (80) 


Dies ist die endgültige Form unserer Differentialgleichung (47), aus der wir u, 
mit Hilfe der Randbedingungen zu bestimmen haben, die für » bei 7’, und 7’, gelten. 

Zunächst suchen wir einen Überblick über die Integralkurven von (80) zu ge- 
winnen. 


$ 8. Die Integralkurven von (80) 

Wegen der physikalischen Bedeutung von » und t kommen nur positive Werte 
dieser Veränderlichen ‘in Betracht. Wir stellen das Steigungsmaß der Integral- 
_kurven fest. 

Es ist dv/dy = oo für t = 0, für die Gerade „A“: vy = y und für die Gerade „D“: 
_t 

u—1l 


, sie entspricht = — oo, vgl. (72). 
M, 


(m* + u) (1—% 
Für uns ist immer << 


u—1l 
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Es ist dv/dy = 0 für v = 0 und für die Kurve „B“: - 


.2 * 2 
v (u 1) (v = y):+ A .) ey 0. (81) 
(t, + y)(m* + (1—Z)) 


Die Kurve B geht durch den Punkt y = 0, v = 0 und liegt im ganzen Bereich y > 0° 
y > Ounterhalb der Geraden A. Es ist nämlich der v-Wert v (B) auf der Kurve B nach (81): 


=. 
(B) = ; (81a) 


1+ 


u (m* u) 1— = 
(tr +y) (m*+y) (:-=) 


der Nenner ist > 1, daher ist v(B)< y. Im Punkt y = 0, v = 0 ist die Neigung der 
Kurve B gegen die t-Achse kleiner als 1. Für hinreichend großes y nähert sich die Kurve B 
asymptotisch der Geraden A. 

Fir vy > 0, t > 0 ist dv/dy > 0, wenn v > y, also links von und oberhalb der Geraden 
A,s. Abb. 1. Im Gebietv > 0,7 > 0, y > 0, y > v > v (B), also zwischen der Geraden A 
und der Kurve B oberhalb der Abszissenachse, ist dv/dy< 0; im Gebiet v> 0, t > 0, 
y > 0, v< also unterhalb der Kurve B, aber oberhalb der Abszissenachse, ist 
dv/dy > 0. 


Dies gibt ein qualitatives Bild der Integralkurven, wie es in Abb. 1 dargestellt 
ist, das wir aber in Einzelheiten noch veri. sern müssen. 

Man sieht, daß der Punkt P:»=0,r=r,d.h.y= 0, ein singulärer Punkt 
ist. Gerade durch P muß die Integralkurve gehen, welche für die Ver- 
brennung maßgebend ist, denn für 
T = T, muß n, = 0 sein, was in unseren v 7) 

Veränderlichen bedeutet, [daß » = 0 sein v 
soll für also für y=0O. Wir 
nennen diese Integralkurve die Kurve ,,C“‘, 
sie muß außerdem durch einen vor- ‚v0 A Pe. 
gegebenen Punkt T=T, n=”, 
gehen, also durch einen Punkt: v=»,, / “7 
tT=T, d. h. y= y, den wir den Punkt Q J an 

nennen wollen. Die Integralkurven von 
(80) enthalten eine willkürliche Integra- “WE 
tionskonstante, sie ist festgelegt, wenn 9 2 

die Integralkurve durch einen vorge- (7T=0) Plt=tr,v=0) ry 
gebenen Punkt gehen soll. In C ist also spp. 1. Erster Überblick über den 
keine Konstante mehr frei, wenn € durch Verlauf der Integralkurven von (80). 
P gehen soll. Also kann im allgemeinen C Hilfskurven sind gestrichelt 
nicht auch noch durch Q gehen, weil Q 

(t= Tt, y= u,) beliebig vorgegeben werden kann. Die in (80) vorkommenden 
Größen: r,, 7, A, m*, u, sind als bekannt anzusehen, «, suchen wir. Da u, ex- 
plizit in (80) vorkommt, hängt der Verlauf der Integralkurven von u, ab. Ändern 
wir uw, in (80), so verändern sich zwar die Geraden A und D nicht. Auch die 
Neigung oo auf der Ordinatenachse und auf A bleibt erhalten und die Neigung 
Null auf der Abszissenachse. Es ändert sich aber die Kurve B und die Neigung 
dy/dy der Integralkurven. Durch.Veränderung von u, kann man erreichen, 
daß die KurveC durch den Punkt @ geht. Dies ist nur für einen Wert 
von u, erreichbar und bedeutet die Bestimmung der Verbrennungs- 
geschwindigkeit. 
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Wir beweisen jetzt, deß Q immer über der Geraden A liegt; das ist gleichbedeutend mit 
dem Beweis dafür, daß v, > y, ist. Es gilt nach (77) », = N,,/D,, und nach (76), (69) und 


(75): 
D 0,€ E r _OTzn Txt 
R (N: ug’ )E ı Nıug ~ Niwa’ 


der Nenner ist kleiner als 1, daher ist v, > y,. Der Nenner liegt aber sehr nahe an 1, d. h. 
Q liegt praktisch auf der Geraden A, die Zahlenwerte von § 7 geben fiir den Nenner in (82) 


ganz rechts: 1 —1-10-%. Wir kénnen also setzen: me 


v, = 7 (82a) 
wenn es nicht gerade auf den Unterschied zwischen v; und y,; ankommt. 

Aus der bekannten Größenordnung von Aktivierungsenergien folgt, daß tr, für 7, x 
300 eine große Zahl ist, etwa von der Größenordnung 30, und daß r, in der Größen- 
ordnung von 5 liegt. Es kann übrigens schon aus dem physikalischen Verhalten von Un- 
verbranntem und Flammenzone der Schluß gezogen werden: r, muß so groß sein, daß die 
Zahl wirksamer Stöße Z für 7, verschwindend klein wird, 7, muß dagegen so klein sein, 
daß Z neben den gaskinetischen Stößen in Betracht kommt. -Mathematisch bedeutet 
dies, daß in der Differentialgleichung (80) das zweite Glied rechts erst für r-Werte in der 
Nähe von r, merklich wird, daß dieses Glied aber für größere t vernachläßigt werden kann. 


Wir setzen: = 


und bekommen für hinreichend großes r aus Gl. (80) 


(84) 
woraus folgt: 
(7, + y)* const. (85) 
Alle Integralkurven in dem uns interessierenden Gebiet v > 0, t > 0 sind also in hin- 
reichender Entfernung vom Punkt vy = 0, y = 0 gerade Linien, die vom Punkt O (v = 0, 


t = (0) weg gerichtet sind. Diese Aussage gilt nur in nächster Nähe der Geraden A nicht, 
. weil dort v x y ist und das zweite Glied auf der rechten Seite von (80) merkliche Werte 
annimmt, auf der Geraden A selbst wird es sogar beliebig groß gegen das erste Glied. 
Es ist also zu untersuchen, ob in der Nähe des Punktes Q (v,, y,) die Näherung (85) 
noch brauchbar ist. Wir berechnen dv/dy für v = v,, y = yı aus (80); für », — y, muß der 
genauere Ausdruck (82) benützt werden, im üprigen können wir v, = y, setzen. Unter 
Beachtung von (77), (73), (64), (28), (73), (44) folgt schließlich für »;: 


+2 + Af, 


denn das zweite Glied ist höchstens von der Größenordnung von 0,07, in der Regel kleinen. 
Wir werden noch zeigen, daß diese Abweichung von dem aus (84) folgenden Wert: 


(86) 


v (yx yı) 
l r 


erst in allernächster Nähe der Geraden A stattfindet. In dieser Gegend sind die Integral- 
kurven nach oben ge krümmt, wie man aus (80) mit Hilfe von (82) und (86) findet. Die 
positive Krümmung vj’ > 0 folgt auch anschaulich aus »; & 1 und aus den Tatsachen: 
Q liegt oberhalb der Ge raden A, links und oberhalb von A ist »’ > 0, die Integralkurven 
sind für y > 1 Strahlen von vy = 0, r = O aus (außer in nächster Nähe von Q). 

pe hee erkennt man, daß in dem Gebiet hinreichend großer y, in dem die Integral- 
kurven Geraden sind, der Faktor e"” dies bewirkt. Im Gebiet hinreichend großer y stimmt 
die Kurve B praktisch mit der Geraden A überein, wie (81) oder (81a) zeigt. Die Kurve € 
geht also von Q aus nach links von der Geraden A weg in Richtung fallender », ist zunächst 
nach oben gekrümmt, verläuft dann längs eines Strahles von» = 0,7 = 0 aus und miindet 
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Die Steigung der Kurve C in der Nähe des Punktes P gewinnt man aus (80). Es 
wird füry< 1, r<1: 


m* 4- 1\2 A m* + 11? 
m* + m* + (87) 
T, 


Diese Differentialgleichung kann elementar integriert werden, man findet: 


B, ist die Integrationskonstante; sie ist Null, weil y = 0 sein soil für » = 0, und weil 
B>0. Das Integral entwickeln wir wir in eine semikonvergente Reihe und bekommen: 


By = Bv—1! (Bv)? + 2! (By)? —3! (Bo)! + —.... (89) 


v = ( verträgt sich mit y = 0, man sieht auch an der Gleichung (80), daß» = 0 eine strenge 

Lösung ist, die t-Achse ist also Integralkurve. Bricht man die Reihe (89) dort ab, wo die 

Glieder zuzunehmen beginnen, so erhält man ein Polynom, das eine hinreichend genaue 

Darstellung von y = f(v) in der Nähe des Punktes P gibt. Die Umkehrung » = f (y) 


ist eindeutig in der Nähe von P, sie lautet: 
og 


€ P 
Es gibt außer der y-Achse nur eine weitere Integrälkurve durch den singulären Punkt P; 
sie ist offenbar die gesuchte Kurve C. In der Umgebung von P ist sie näherungsweise 
durch (90a) gegeben. Da wir (90a) aus der Näherungsgleichung (87) gewonnen haben, 
müssen wir, um C in der Nähe von P genauer zu bekommen, mit dem Ansatz: 


y=y+t... 


(90b) 
in die Differentialgleichung (80) eingehen. Man findet: 
> 
‚v=y+By+B —1)y°+...; (91) 
(m* + u) (1 _ =) 


Bist in (87) erklärt. (91) ist die gesuchte Reihenentwicklung der Kurve C in der Nähe des 
Punktes P(y=0,»=0). Die Kurve C berührt demnach die Gerade A im Punkte P 
und ist in der Nähe von P nach oben gekrümmt. Mit steigendem y nimmt die Krümmung 
ab, der Faktor von 4? in (91) ist negativ. Nach (73) erhält man mit den in § 6 angegebenen 
Zahlen A ~ 5-107 und B »1,0, wenn r, x 5 gesetzt wird. 


Keine der Integralkurven oberhalb A überschreitet die Gerade A, weil in diesem Ge- 
biet überall »’ > 0 ist und weil »' — co geht auf A. Dies gilt also auch für die Kurve €. 
Die Integralkurven unterhalb A überschreiten alle die Gerade A und schmiegen sich ihr 
bald nahe An. Allgemein muß sich jede Integralkurve oberhalb A, sobald sie der Geraden A 
hinreichend nahe gekommen ist, der Richtung von A um so mehr anschmiegen, je größer y 
ist. Bei den höher gelegenen der Integralkurven, die aus dem Gebiet unterhalb A kommend 
auf A zulaufen, geschieht das Umbiegen bei A, zuerst in die Senkrechte und dann in die 
Richtung von A nach außen, so plötzlich, daß sie nur aus zwei Geradenstücken zu bestehen 
scheinen: der Geraden von außen aufO zu bis zum Schnitt mit A, und der Geraden A vom 
Schnittpunkt an nach außen, die beiden Stücke bilden einen sehr großen, fast 2 betra- 
genden Winkel miteinander. Man bekommt so als Verbesserung zu Abb. 1 die genauere 
Darstellung der Abb. 2. In der Abbildung sind zur Erläuterung einige der Integralkurven 
numeriert, damit der Verlauf deutlich wird. Längs der Geraden A fallen die von unterhalb 
A kommenden Integralkurven für größeres y alle praktisch mit A zusammen. 


Aus den Gleichungen (32), (67), (68), (69), (68a), (76), (77) schließen wir auf: 


ug t Dy v—y, 


4 3 
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das läßt sich wegen (69) und (77a) umformen in: 


m 
AT, = — 
T 


folglich ist 7,» O (108) at. Da T, in der Flamme, also für T Z 1800° abs d.h. 


t S 6, sehr hohe Werte hat von der Größenordnung von 105, so wird v — y.~ O (1) für 
diesen t-Bereich und zwar ist v > y auf der Kurve C, also 7’, > 0, wie es sein muß. Für 
y < 1, aber nicht zu klein, etwa y ~ !/,, überwiegt das zweite Glied auf der rechten Seite 
von (80) das erste stark, bei y ~ !/, etwa um den Faktor 20; » kann man in dieser Gegend 
aus (91) berechnen. Die Kurve C entfernt sich in diesem Bereich von der Geraden A, denn 


(92) 


L 


Abb. 2. Schematische. Darstellung der Integralkurven von Gl. (80): ausgezogene Linien. 
Hilfskurven sind gestrichelt. Die Kurve C ist dick ausgezogen 


Pfr=t, 


ihre Krümmung ist positiv, s. (91). Auch für kleiner werdendes 7’, also größeres x und y, 
behält der Temperaturgradient 7’, zunächst noch hohe Werte; y — y wird daher zunächst 
immer größer. Geht man zu immer größeren y-Werten, so wird das zweite”Glied bald 
kleiner als das erste, weil » — y steigt und vor allem, weil e-v stark abnimmt. So kommt 
man in das Gebiet, wo t und damit y so groß wird, daß das zweite Glied in (80) gegen das 
erste vernachlässigt werden kann; denn vy — y ist dann immer noch von der Ordnung 1, 
während r in die Größenordnung von 20—40 kommt. Dies gibt das gradlinige Stück der 
Kurve C, dessen Richtung die eines Strahles vom Punkt O » = 0, r = 0) aus ist. Die 
Berechnung von »; zeigt, daß in (80) für großes y das zweite Glied erst dann neben dem 
ersten wieder merklich wird, wenn v ~ y geworden ist. Dann erst verläßt die Kurve ( 
die genannte Strahlrichtung, um sich der Geraden A anzuschmiegen. Wenn r schon prak- 
tisch gleich r, ist, und» — y y yın, won < 1, aber 


=en (10°) 


— wegen des sehr kleinen Wertes von & ist das möglich —, so ist die Kurve C immer noch 
geradlinig in der Strahlrichtung von O aus. Erst wenn » — y = y, e geworden ist, wird 


di 
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das zweite Glied in (80) neben dem ersten merklich und die Kurve C biegt um in die Rich- 
tung von A. Zugleich wird 7,» T,,, wie man aus (92) findet; der Punkt Q ist damit 3 
erreicht. 

Wegen ex O(10-°) liegt @ da, wo der gradlinige Teil der Kurve C 2 \ 


die Gerade A trifft. Daraus ergibt sich, daß C im gradlinigen Teil die Neigung 
Die gesuchte Kurve C hat also folgende Eigenschaften: Sie geht von P (vy =0, 
y= 0) aus, und steigt bis zum Punkt Q (v, = y, = 1,—r,). Sie berührt die 
Gerade A (y= y) im Punkt P und hat dort die Gestalt: »= y+ By? +.., wo 
B durch (87) gegeben ist. Oberhalb P ist ihr Anstieg zunächst wesentlich durch 
das zweite Glied von (80) gegeben. Sie schmiegt sich bald dem Strahl an, der von 
0 (v= 0, t= 0) aus nach dem Punkt Q gezogen ist, und fällt bis@ mit diesem Strahl 
zusammen. Da t,, tT, bekannt sind, kennt man die Punkte P,Q. Die Kurve C 
geht nur dann wirklich durch Q, wenn u, geeignet gewählt ist, sonst trifft sie die 
Gerade A in einem von Q verschiedenen Punkt. Wir wenden uns jetzt der Berech- 


nung von u, ZU. 
> § 9. Analytische Berechnung der Verbrennungsgeschwindigkeit $ 
Als Näherungsfunktion für C setzen wir an: 
ve + by) + + by, (94) 


dieser Ansatz erfüllt die Bedingungen: gerade Linie für „> 1, (b., b, sind so zu 
bestimmen, daß die Gerade ein Strahl von O durch Q wird), gradliniger Anstieg 
bei P (die Berührung der Geraden A im Punkte P gibt zwei weitere Bedingungen 
für die Konstanten), gute Darstellung des qualitativen Verlaufes zwischen P und 
dem Strahl von O aus, vor allem in dem wichtigen Gebiet nicht zu großer y. Das 
letztere erkennt man so: in dem eben genannten Stück überwiegt das zweite 
Glied rechts in (80); »—y ist dort von der Größenordnung von y, (t,-+ y)"! 
kann man nach y/t, entwickeln, die beiden Quadrate im Zähler und Nenner (ohne 
1/u;) stammen von (M/.M ,)*, das nie sehr von 1 verschieden ist. Das zweite Glied 
in (80) hat also im betrachteten Bereich praktisch die Form: const - y e~¥. Aus 
den eben genannten Bedingungen findet man: 


Ty ‘Ty 
b=1,( —)= 4, bs, b, = — = —1,),, 
b,—6b,+6,=1, 
dies gibt endgültig als Näherungsformel für die Kurve C: Me 


Vo (1— bs (1 + T,))— 5s T,] e’+ bs (y + T,); b,= 1— 


Es wäre nicht schwer, den Ansatz (94) so zu ändern, daß die im vorigen genannten 
Bedingungen für vg erfüllt bleiben, und daß bei y <1 die genauere Näherung 


(91) gilt. Ein brauchbarer Ansatz dazu ist 
a +ay+a; = - 
vo = ay Ao (95a) 


aya+l 
= 5 


=> 


ing 1, 
k der 
Die 
‘dem, | 
rve 
= 
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Das im folgenden beschriebene Verfahren, wu, mit der EN (94) 
zu bestimmen, läßt sich ebenso auf (95a) oder andere zweckmäßige Ansätze an- 
wenden. Wir gehen darauf nicht ein, sondern bleiben bei dem einfachen Ansatz , 
(94), (95). 


Nach unseren Überlegungen können wir erwarten, daß die Funktion (95) eine 
brauchbare Näherungsdarstellung der Kurve C sein wird, und zwar derjenigen 
Kurve C, die zum richtigen Wert von «, gehört. Denn für beliebiges, aber fest 
gewähltes wu, geht die Kurve C nicht durch den Punkt Q, entspricht also dann nicht 
unserer Funktion v,. Nur für einen ganz bestimmten, eben den richtigen, Wert 
von u, geht C durch Q und wird also nur für diesen Wert von u, durch v, näherungs- 
weise dargestellt. »., hat nach unserem Ansatz das für die „richtige‘* Kurve C 
notwendige Verhalten genau für y—0 und y— y, hat auch im Zwischengebiet 
ungefähr die richtige Form. Hätten wir also die „richtige‘‘ Kurve C, so müßte 
sie bei P undQ mit». übereinstimmen (in der Näherung, in der wir die Bedingungen 
bei P und @ erfüllt haben), aber auch im Zwischengebiet muß sie sich gut an vg 
anpassen. Dies ist eine Forderung an die Differentialgleichung, 4, muß einen solchen 
Wert haben, daß »., der Differentialgleichung möglichst gut genügt. Wir haben 
nämlich in vg weh Parameter mehr frei, daher bedeutet jede über die bei Pund @ 
schon e ‚rfüllten Bedingungen hinausgehende Forderung, v., solle Eigenschaften der 
richtigen Kurve C haben, eine Festlegung des Parameters «,, der in der Differential- 
gleichung noch wählbar ist. 


Dies ist ein Verfahren, das &rundsätzlich auf Parameterbestim- 
mung auch in anderen Differentialgleichungen anwendbar ist: Das 
Problem sei so beschaffen, daß es nur bei bestimmter Wahl des Parameters eine 
Integralkurve gibt, welche den gestellten Bedingungen genügt. Man setzt als 
Näherungslösung eine Funktion an, welche die gestellten Bedingungen erfüllt und 
möglichst weitgehend die Eigenschaften der gesuchten Integralkurve hat. Ein- 
setzen dieser Näherungslösung in die Differentialgleichung gibt die (näherungs- 
weise) Bestimmung des Parameters. 


Bei unserem Problem haben wir 2 Möglichkeiten: 


I. Wir benützen die Entwicklung (91), die ja aus der Differentialgleichung (80) 
gewonnen ist, und verlangen, daß sie auch in höheren Gliedern mit der Entwick- 
lung von v,, übereinstimmt. Damit ist erreicht, daß die Näherungskurve »,, auch 
in der weiteren Umgebung des Punktes P, auf die es ja physikalisch wesentlich 
ankommt, mit einer Kurve C übereinstimmt, die zu dem so berechneten u, gehört 
und daher nicht notwendig durch Q geht. Da aber »., nach unserer Überlegung 
die „‚richtige‘‘ Kurve C gut darstellen wird, so wird das aus der Forderung I be- ° 
rechnete wu, nahe mit dem richtigen «, übereinstimmen. 


II. Wir verlangen, daß die Näherungslösung »,, an einer physikalisch wichtigen 
Stelle die Differentialgleichung (80) streng erfüllt. Das heißt, wir wählen eine 
Temperatur 7', die im Gebiet starker Umsetzung liegt, also näher an 7',alsan 7’, 


= 27.+7 
etwa die zugehörigen r- und y-Werte gilt dann: T x 57 
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der Gegend von 0,6 bis 0,8. Wir setzen, da wir ja % willkürlich wählen können — 
nur zweckmäßig muß die Wahl sein — ' 
An dieser Stelle ¥ soll ».. der Differentialgleichung (80) streng genügen, das soll 
unsere Forderung sein. Zu dem so berechneten uw, gehört dann eine Integralkurve 
durch den Punkt %, v. (y), welche in diesem Punkt die gleiche Neigung hat wie 
vo. Sie wird weder genau durch P noch durch Q gehen, aber annähernd, wenn 
yo eine gute Näherungskurve ist. 
Wir leiten jetzt die Formeln für die Verbrennungsgeschwindigkeit ab, die aus den 


Näherungsannahmen I oder II folgen. 
I. In der Nähe von P gilt nach (95): 


ul, _ 2\ te 
Vergleich mit (91) gibt 2 Bedingungen fiir uj, weil wir 2 Glieder vergleichen. Es ist eine 
Probe auf die Güte unserer Näherung, daß die beiden Bedingungen natürlich an- 
nähernd denselben Wert von u, gebei? müssen. 
Der Faktor von y? gibt wegen (87): 
‘ m* + 1\2 
2A +n) ~ 


u) = 


und der Faktor von y° 


‘ m* + 1\2 
% 34 e 1— 
2 T 


= 
Das Verhältnis der zwei Werte von u? ist: u, ¥ 
2 
(98e) 
e 
3 2 (u—1) 1 


egen t, > 1. Das letzte Glied in der Klammer kann weggelassen werden, es ist < 0,03, 


das zweite Glied ist < — , die Klammer ist also wenig veränderlich (vgl. die Zahlwerte 


x r 
in Tabelle 1!). Zahlenmäßig gilt: ut 
x 1,1 bis 1,2. (98d) 


Ui(1) 
Das Näherungsverfahren hat also diese erste Probe gut bestanden. 
Wir können feststellen: Methode I gibt als Formel für die Verbrennungsge- 
schwindigkeit: 
V24,: 4(2a,)? V A, B* 
m*+u 2 mı 
u = —; 
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Wir werden noch zeigen, daß die Methode II praktisch dieselbe Formel gibt. 
Nach (98e) ist uw, wesentlich durch e”’r/? bestimmt, es ist groß für möglichst 


kleines 7,—= 7; , also für hohe Temeparur im Verbrannten und kleine Akti- 

vierungsenergie, %, wächst monoton mit 7’, und auch mit 7’,, der Temperatur im 


Unverbrannten. Für 7, gilt nach dem Energiesatz (77a): 
E Cy, 01 E fi+2+4Af, m 2 1 
m*+i 
(99) 
Bei festem 7, = cn ist 7, um so kleiner, je größer 4 E ist, je kleiner die Zahl 


der Freiheitsgrade des Moleküls, auch die Zah} der durch Temperatursteigerung 
anregbaren Freiheitsgrade, und je größer m* ist, das Molzahlverhältnis Sauer- 
stoff/Brennstoff. Die Aussage über die Abhängigkeit von m* hängt aber wesent- 
lich an der Voraussetzung einer Reaktion zwischen gleichen Molekülen; für Reak- 
tionen Brennstoff-Sauerstoff, wie sie bei der Kohlenwasserstoffverbrennung eine 
Rolle spielen, gilt sie in dieser Form nicht. g’/E ist übrigens immer > 1, weil die 
Reaktionswärme q’ pro Molekül bei Verbrennungen größer ist als die Aktivierungs- 
energie E pro Molekül, die freiwerdende Energie größer ist als die hineingesteckte, 
sonst entsteht keine Verbrennung. 

II. Die zweite Methode der Berechnung von u, verlangt Erfüllung der Differential- 
gleichung (80) durch vy, für y=% = r,/3. Wir berechnen aus (95): 
ve (t, + y) —% [y®(—1 + + 1)) + yt, (—1 + + 2)) + 
Also wird mit der Bezeichnung 7, = v, (9): 


4 1 1 7 4 T Ie 
Vo (tT + ¥) —% =e 5 (1 + 3 =| 9 =): 


(100a) 


und: 
E 00€ 


. . . . eae . -Tr/3 Tr 
dieser Ausdruck ist übrigens, wie zu verlangen, immer positiv, weil e 3 +1)<1 


fiir alle positiven t,. Wir brauchen noch: 


- [(1- (F +2) 1 =) 1 e 3 t 1 . 


(1004) 
So ergibt sich durch Einsetzen in (80): Par 
2 + : 2 01) 
u ; (1 
= | 3(u— DIR 
-tr/3 (Tr 
Im +n) (1-«"?G +1)) 
das Verhältnis der zwei Werte von u? wird: 
ae = Tr 4 = tr 


m* + u 


alle | 
Eins. 


V 


die | 
bre} 
ansa 
Offer 
in d 
Pun 

V 
ihrer 
gült 
keit: 


Mole] 
stoff) 
2a, 
wärn 
dime 
Einh 

D 
Wir \ 
mäß 
schor 
Ener 


ist n 


= 
0, ¢ 


mens 


ran 
= 
4 
. 
= 
a 
IE 
| 
- — 
1 ) 
3 (1 1 
3 ) 


01a) 


= 


K. Bechert: Zur Theorie der Verbrennungsgeschwindigkeit 219 


alle Klammerausdrücke sind wenig veränderlich und liegen dem Wert nach nahe bei 
Eins. 


Wir können feststellen: 


Ui 


x 1,2 bis 1,4; (101b) 


Ui (1) 


die beiden Methoden geben also praktisch dieselbe Formel fiir die Ver- 
brennungsgeschwindigkeit, was wiederum zeigt, daß unser Näherungs- 
ansatz gut ist. Die Übereinstimmung ist so gut, wie man sie nur erwarten kann. 
Offenbar kann man allgemein die Güte der Näherung für die Parameterbestimmung 
in dieser einfachen Weise beurteilen, daß man den Parameter an verschiedenen 
Punkten bestimmt. 


Wir benutzen im folgenden den für numerische Rechnung einfacheren Wert 


un. Wir fügen den Faktor / 2* 1,4 hinzu, weil die Näherungsformel für «711, 
ihrer Ableitung nach etwas besser sein wird als die für %(r, und erhalten als end- 
gültige Form der Näherungsformel für die Verbrennungsgeschwindig- 


keit: 
‘ a\/, /3 
TI 


Wir stellen die Bezeichnungen nochmals zusammen: t, = E/kT',; t, = E/kT ; 
T,, T, Temperatur im Verbrannten, Unverbrannten; E' Aktivierungsenergie pro 
Molekül Brennstoff; k Boltzmannkonstnte; m* Molzahl O,/Molzahl O, (Brenn- 
stoff), beides im Unverbrannten; u Verhältnis der Molekulargewichte 0,/0;; 
2a, m, Molekiildurchmesser, Molekülmasse des Brennstoffs; g’ Reaktions- 


wärme pro Molekül Brennstoff; 2, Wärmeleitfähigkeit des Gemisches// T. Alle 
dimensionierten Größen sind im gleichen Maßsystem zu messen, z.B. in CGS- 
Einheiten. 

Der erste Faktor in (102) ist 2/ A,, vgl. (98e), er bestimmt die Dimension. 
Wir wollen nun zeigen, daß der Bau der Formel (102) für «;, also die dimensions- 
mäßige e “Abhängigkeit zwischen u, und den physikalischen Bestimmungsstücken, 
schon aus rein qualitativen Überlegungen ohne viel Rechnung folgt. Nach dem 
Energiesatz ist (vgl. (77a)): 


06, = (77a) 
wenn wir ihn für die beiden Grenztemperaturen 7',, 7’, hinschreiben. Andererseits 
ist nach dem Energiesatz dimensionsmäßig (vgl. (25)): 

Lepr T—rq fZdx; (103a) 


r= Zahl der Brennstoffmoleküle, die pro Reaktion verschwinden. Die neben 


0,%,c,, 7 auftretende Größe rq’ f Zdx scheint darauf hinzuweisen, daß die di- 


Zd 
mensionslose Größe 


eine Rolle spielen könnte; da aber (vgl. (31a)) 
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ist, so bekommt man aus dem Verhältnis nur die dimensionslose Größe 
q 


7)” die nach dem Energiesatz gleich 1 1 ist. Aus (31a) folgt 
(7, — 


d (n, u) 
~rZ. (31a) 


Durch Division erhält man dimensionsmäßig aus (103a), (31a): HE 


Jew Sor 1 
rZ 


weil n,  n,, und u wu,. Das gibt wegen (77a): ur 


T) _rZ AM + 2. (1034) 
auf der rechten Seite kann noch ein dimensionsloser Faktor f stehen, der von di- 
mensionslosen Größen abhängen kann, also von 1,, t,,m*, u. Nach der kinetischen 
Gastheorie ist 
T 


(2 a)? m 


(103e) 


bis auf einen Faktor der Größenordnung Eins; 2a, m sind mittlere Werte für 
Moleküldurchmesser und Molekülmasse der Gasmischung. So wird aus (103d), 
wenn wir ”,, = Pı,/kT', beachten: 


tit, (2a) Ymg' pi, tit,V t 


Die Gln. (103d, f) zeigen, daB die Formel fiir u, wesentlich vom Ansatz für die 
Stoßzahl Z abhängt. Für diese Anzahl Z der pro cm? und sec zur Reaktion führen: 
den Stöße konnen wir allgemein ansetzen (vgl. (28)): 


(103g) 


n, ist wie bisher die Zahl der Brennstoffteilchen im cm?; es ist angenommen, daß 
die Reaktion nur durch die Teilchen 1 zustande kommt; der allgemeinere Fall 
wird in der nächsten, Arbeit behandelt; der Faktor 7? gestattet, auch das Mit- 
wirken innerer Freiheitsgrade bei der Reaktion zu berücksichtigen t); für s Schwin- 
gungsfreiheitsgrade, die bei der Reaktion mitspielen, kommt ein Faktor r?°/(2 s)! 
zum Ansatz für Z in (28), wenn es sich um Zweierstöße handelt; g ist ein Faktor 
von der Dimension cm**~* sec=!. Man erhält so aus (103f): 


F ist dimensionslos und hängt von Tp, Ty, m*, ab. sea 


Bei unserem — (28), der ein spezieller Fall von (103g) ist, RE „=. 


B=0; 9 ~ (2a) \; = , der weggelassene Faktor ist von der Größenordnung 1. 


11) Vgl. die Zusammenstellung der Formeln dazu bei H. J. Schumacher, Chemische 
Gasreaktionen, Verlag Steinkopff, Dresden 1938, S. 20. > 
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öße Setzen wir (103g) in rg ein mit dem angegebenen Wert von g und berück- 
sichtigen, daß A = AT , so wird für « = 2, ß= 0: 


| 
la) ; > u (2a) A, E? ( (103i) 


k's Vm q 


was der Form nach mit (102) übereinstimmt. Daß im Exponenten eher r, 
Se) als 7, zu stehen hat, wird man daraus vermuten, daß erst bei hoher Temperatur die 

Zahl der Reaktionen stark ansteigt; aber ohne Rechnung kann man natürlich 
weder diese Vermutung prüfen noch die Funktion f/r wirklich angeben. Aus 


(103i) wird für unseren Ansatz (28) wegen A, u ——— 
3d) (2 a)? Ym 


(1-3) 


. ist von der Größenordnung —- ; — 2 7 ist sehr genau gleich Ey,), - 10!*/.M, also von 


)3 
9 der Größenordnung 10!*/M a wo M das Molekulargewicht der Gasmischung 
bedeutet. 

= (103d)-und alle daraus abgeleiteten Gleichungen bis (103k) folgen also dimen- 


sionsmäßig aus den Grundgleichungen. Dies ist selbstverständlich, denn alle 
einer naturwissenschaftlichen Betrachtung zugrunde liegenden Gleichungen können 
dazu benutzt werden, dimensionsmäßige Zusammenhänge aufzustellen; und weil 
3f) in unserem Problem keine Apparatdimensionen vorkommen, muß es möglich sein, 
die Verbrennungsgeschwindigkeit dimensionsmäßig aus den Grundgleichungen 
allein zu gewinnen. 


die Bis jetzt sind wir von der Annahme nicht abgegangen, daß wir die Diffusion | 
m vernachlässigen oder besser: summarisch berücksichtigen können, ohne die Glei- 

chungen zu benutzen, welche die Diffusion beschreiben. Wir wollen jetzt zeigen, 
3g) daß der Zusammenhang (103d) auch noch gilt, wenn die Diffusion be- 

rücksichtigt wird. Die Gleichungen, welche diesen allgemeinen Fall regieren, 
daß stehen z. B. im Buch von Jost}*) (die Impulserhaltungsgleichung ist dort noch 
Fall hinzuzufiigen). Erhalten bleiben die Beziehungen (77a), (103a, b); an Stelle von 
Mit- ‘ (31a) tritt dimensionsmäßig: 5 
vin- 

' d dn, d(n, u) 

DE) + (1031) 

D ist die Diffusionskonstante. Da also p> ~ n;,U,, so folgt aus (103a, 1): 
| x 

A 

3h) (103m) 

eine Dimensionsbeziehung, die aus der kinetischen Gastheorie bekannt ist. Aus 
= 2; (103a, 1) ergibt sich so wegen (77a): S u 
g 1. _ArZ(T, —-T)) +1)_ (103n) 

1 Cpr Rıı M1 

sche 12) W. Jost, Explosions- u. Verbrennungsvorgänge in Gasen, Springer, Berlin, 1939, 
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das ist wirklich Gl. (103d). Die letzte Form der Gl. (103n) ist besonders durch- 
sichtig. Die Gleichungen (103d, f) gelten also für die dimensionsmiäßige 
Abhängigkeit der Verbrennungsgeschwindigkeit von den physika- 
lischen Größen ganz allgemein, sie folgen für stationäre Verbrennung allein 
aus den Erhaltungssätzen, wenn die Reaktion nur durch die Teilchen 1 zustande 
kommt; dagegen ist der dimensionslose Faktor f in seiner Form nur durch wirk- 
liche Rechnung zu finden. (102) ist ein spezieller Fall von (103d, f) und ist bedingt 
durch unseren Ansatz (28) für die Stoßzahl Z. ed ths 


10. Vergleich mit der Erfahrung 


Zunächst vegleichen wir unsere Formel für die Verbrennungsgeschwindigkeit 
mit der von Nusselt!?) aufgestellten. Die Formel von Nusselt ist aus den An- 
nahmen gewonnen: Es gibt eine ,,Ziindtemperatur“ 7’, oberhalb deren das Gemisch 
sich entzündet, unterhalb deren es aber nur langsam reagiert; zwischen 7’, und T’, 
ist der Temperaturanstieg allein durch Wärmeleitung von der Flamme her be- 
dingt; von 7’, bis 7’, steigt die Temperatur linear; die Reaktionsgeschwindigkeit 
ist temperaturunabhängig im ganzen Bereich von 7’, bis 7',. In unserer Bezeich- 


nung lautet die Formel von Nusselt: 2 
A(T, — _ (104) 


uy = ; 


s ist die Zahl der Moleküle des Endproduktes der Reaktion, die bei der Einzelreak- 
tion entstehen (für die Ozonverbrennung ist das Endprodukt O,, und s= 3, Be- 
zeichnung wie in (29)); 0, € sind Mittelwerte der Dichte und spezifischen Wärme des 
Gemisches; n,, ist die Zahl der Moleküle des Endproduktes der Reaktion pro cm? 
im Verbrannten. Wegen p, ~ p, (vgl. (50)) und wegen (71) folgt: 


T 
Me, = M1 (m* + 1). 
Das gibt der Gleichung von Nusselt die Form 
7 
— Ts %,(m* + 1) (104a) 


Das ist aber die Form (103n), die ja nach unserer Überlegung allgemein für 
alle Gleichungen gilt, welche die Verbrennungsgeschwindigkeit beschreiben. Es 
folgt durch Division der Gln. (104a), (103n): 


es bleiben also Verhältnisse von Temperaturen und ein Molenbruch übrig, eine 
dimensionslose Funktion, die bei unseren Dimensionsbetrachtungen zu den Gln. (103) 
immer hinzukommen kann. 

Die Schlüsse über die Abhängigkeit der Verbrennungsgeschwindigkeit von 
0, €, Z, die aus der Nusseltschen Formel gezogen worden und die an der Erfah- 
rung bestätigt worden sind, gelten folglich allgemein, gelten für jede Theorie der 


13) Wir folgen der Darstellung, die W. Jost in seinem oben zitierten Buch gibt, S. 108. 
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Verbrennungsgeschwindigkeit und sind nicht charakteristisch für die Voraus- 
setzungen, unter denen die Formel von Nusselt gewonnen worden ist. Wir nennen: 
1. Ersetzt man in einem brennbaren Gemisch ein chemisch inertes Gas.durch ein 
anderes von etwa gleichem A, aber verschiedenem 5 €, so muß dieVerbrennungs- 
geschwindigkeit für das Gas mit größerem o € kleiner sein. (Unsere Formeln (103) 
haben statt 5© den Ausdruck 9, ©,,, wofür die Begründung in unseren Uber- 
‚legungen zur Energiegleichung (77a) liegt.) 2. Nach (104) und allgemein nach 
(103n) muß das Verhalten der Verbrennungsgeschwindigkeit allein durch A be- 
stimmt sein, wenn in einem Gemisch ein inertes Gas durch ein anderes von gleichem 
0 ¢ ersetzt wird. 

Dieser Schluß ist nicht durch die speziellen, sehr schematisierenden Annahmen 
der Theorie von Nusselt bedingt, ist vielmehr allgemeingültig und folgt mit 
Notwendigkeit aus den Grundgleichungen. Gegen den Schluß ist eingewendet 
worden 4), daß He eine über 6mal größere Wärmeleitfähigkeit hat als A bei gleicher 
Wärmekapazität, und daß daher die Verbrennungsgeschwindigkeit in Gemischen, 
denen He zugesetzt ist, wesentlich größer sein müßte, als in Gemischen mit A- 
Zusatz, während die Messungen nur eine geringe Erhöhung der Verbrennungsge- 
schwindigkeit zeigen. Beobachtet ist in einer Mischung von 10%, CH, mit 90% 
einer „Atmosphäre“ , die zu 21% aus O, und zu 79%, entweder aus H oder A bestand, 
daß die Verbrennungsgeschwindigkeit in der Mischung mit A 59 cm/sec betrug, 
in der Mischung mit He 82 cm/sec. - 

Mir scheint der Einwand nicht zwingend. Die Wärmeleitfähigkeit A in den 
Formeln (103n) und natürlich auch in (102) meint die Wärmeleitfähigkeit des Ge- 
misches. Bei einer Rechnung mit einer Formel für die Wi ärmeleitfähigkeit von Ge- 
mischen, die mir aus einer einfachen gaskinetischen Überlegung zu folgen scheint, 
bekam ich für das theoretische Verhältnis der Verbrennungsgeschwindigkeiten 
bei den angegebenen Mischungen 2,1 statt des beobachteten 82/59 = 1,4, also 
einen 50%, zu hohen Wert. Die Wärmeleitfähigkeit ist bei diesen Gomme leider 
nicht gemessen worden. 


Eine Verallgemeinerung der üblichen elementaren Ableitung der Wärmeleitfähigkeit 
aus kinetischen Betrachtungen führt auf die Beziehung: 


Ay 


a, + a, Ns (“ + my, 
1+; Er 2a, ) | ms + ny 2a, | m; m 


( ay ms + a, + 
2a, y= 2a, ) | 
Ay 
Ni (* + mk 
= nk | mi 
dabei sind die 4; die Wärmeleitfähigkeiten der Einzelgase, n; die Molekülzahlen im cm’, a; 
die Molekülradien, m, die Molekülmassen. Bei der Berechnung für die ge nannten Methan- 


mischunge n sind die Werte verwendet worden (in A): 2acm, = 3,34, 2ao, = 3,0, 
2ape = 2,0, 2ay = 2,9. Für alle 2, ist gemäß der kinetischen Gastheorie angenommen 


worden, daß 2, = 4, V 7 sei; als 4), wurden verwendet die Werte (in CGS): CH al ». 


0, 10%, He 6,1 + 10%, A 0,92 + 108, 


am Vel. die Darstellung bei Jost, S.110/111. 
Ann. Physik. 6. Folge, Bd. 4 
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Einen Fehler von 50% wird man nicht zu groß finden, wenn man bedenkt, daß die 
Verbrennungsgeschwindigkeitsmessungen in der Regel nicht sehr genau sind, und daß die 
7; nicht genau proportional / T sind. Qualitativ ist aber an der Formel für A leicht zu 
erkennen, "warum die Messungen der Verbrennungsgeschwindigkeit in Mischungen von 
Wasserstoff mit „Atmosphären‘“ aus 21% O, und 79% He oder Ne oder A ,,keinen wesent- 
lichen Unterschied zwischen den Edelgasen“ erkennen lassen; A, ist für H, besonders 
groß und m besonders klein; der Beitrag des Wasserstoffs zu A wird dann so groß, daß sogar 
der Beitrag des He keine wesentliche Änderung von gegenüber den Werten für die anderen 
Edelgase bewirken kann. 

Übrigens zeigt die Formel für 7, daß die Wärmeleitfähigkeit einer Gasmischung sich 
nur dann einigermaßen nach der Mischungsregel aus den Wärmeleitfähigkeiten der Kom- 
ponenten berechnen läßt, wenn die Molekülradien und die Molekülmassen der merklich zu 
/ beitragenden Molekülarten nicht sehr voneinander verschieden sind. Leider fehlt mir 
experimentelles Material, um diese Behauptung zu prüfen. 

Die Formel: von Nusselt benützt den Begriff der Zündtemperatur. Jost 
legt in seinem Buch ausführlich dar, daß dieser Begriff eine Behelfsannahme mit 
wesentlichen Mängeln ist. Immerhin gestattet dieser Begriff (wie bekannt), aus 
(104) auf die Existenz von Zündgrenzen zu schließen. Wenn nämlich das Gemisch 
nicht stöchiometrische Zusammensetzung hat, so liegt 7’, unterhalb des Höchst- 
wertes, den 7’, haben kann, während die Zündtemperatur durch Abgehen von der 
stöchiometrischen Zusammensetzung sicher steigt (das Gemisch zündet schlech- 
ter), die Differenz 7’,— 7’, wird also kleiner und damit wu}. Wenn u, so klein ge- 
worden ist, daß die Verbrennung nicht mehr wandert, ist die Ziindgrenze erreicht. 
Für den Fall der Reaktion zwischen gleichen Teilchen, den wir in dieser Arbeit 
betrachten, und auf den sich die Formel (102) bezieht, gibt es den Begriff stöchio- 
metrischer Zusammensetzung nicht. Nach (77a) steigt 7’, mit n,,, also mit der 
Konzentration des Brennstoffes, dabei nimmt t, ab, und nach (102) wäthst u, 
mit abnehmendem 7,; w, wächst mit wachsendem n,, also immer weiter. Ein Ab- 
nehmen der Verbrennungsgeschwindigkeit zu beiden Seiten einer bestmöglichen 
Zusammensetzung ergibt sich erst, wenn man Reaktionen zwischen ungleichen 
Teilchen betrachtet; in der nächsten Arbeit soll gezeigt werden, daß dann theore- 
tisch bei Luftüberschuß die Verbrennungsgeschwindigkeit mit steigendem Brenn- 
stoffzusatz steigt, daß dieser Anstieg noch ein Stück über den stöchiometrischen 
Punkt hinaus anhält — wie es beobachtet ist — und daß im Gebiet des Brennstoff- 
überschusses die Verbrennungsgeschwindigkeit nach einer Spanne des Ansteigens 
wieder fällt. Die Behelfsannahme einer Zündtemperatur ist dazu nicht nötig, wir 
haben sie ja auch in dieser Arbeit nicht gebraucht. 


Aus (102) folgt, daß uw, nicht sehr mit der Temperatur 7’, veränderlich ist; denn 
nach (99) ist 7, nicht sehr veränderlich bei veränderlichem r7,, weil der zweite Aus- 
druck im Nenner von r, neben 1/r, in der Regel 4- bis 5mal größer ist als 1/r,, vgl. 
die folgende Tabelle 1. Diese verhältnismäßig geringe Veränderlichkeit von u, 
entspricht der Erfahrung (Jost, S. 112), und widerspricht der Formel von Nusselt, 
nach der «, für 7, = 7'z sogar beliebig groß werden müßte. In Wirklichkeit ver- 
hält sich w, völlig regelmäßig, auch wenn man weit über die üblichen Zündtempera- 
turen hinausgeht. «, wird erst dann stark von 7’, abhängig, wenn die Brennstoff- 
konzentration so klein ist (großes m*), daß r, in die Größenordnung von t, kommt 
(vgl. hierzu den Bau der Formel (99) für r,!), in diesem Gebiet nimmt «, mit ab- 
nehmendem 7', (und daher abnehmendem 7',) schnell ab und wird unmeßbar klein 
(Zündgrenze). Mathematisch hört diese Abnahme erst dann auf, wenn in (102) im 
Nenner der Ausdruck / 1— “7 sehr klein wird, also in nächster Nä 
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ist aber zu beachten, daß (102) eine Näherungsformel ist, und daß die ihrer Ab- 
leitung nach zuverlässigere, aber kompliziertere Formel (101) keine solche Null- 
stelle im Nenner auf der rechten Seite hat. Unsere ganze Näherungsrechnung 
machte übrigens immer wieder davon Gebrauch, daß r, wesentlich größer ist als 
T,, Wie es ja für Flammen stets zutrifit. Mit anderen Worten: Der Grenzfallz, —r, 
hat keine physikalische Bedeutung und ist nach der Ableitung der Formeln auch 
mathematisch nicht zulässig, weil sie für diesen Fall nicht zuständig sind. 


Wir wollen nun zeigen, daß (102) die Messungen an der Ozonverbrennung 
recht gut quantitativ darstellt. Aus den in $3 genannten Gründen haben wir E 
als Parameter zu betrachten, der aus einem Beobachtungswert von u, zu bestim- 
men ist, nur die Größenordnung von £ ist von vornherein aus der Atomphysik 
bekannt » 1 eVolt = 1,6 - 10-12 erg, E/k ~ O (10°). A, könnten wir berechnen, 
wenn uns der Molekülradius von O, bekannt wäre, sicher ist er größer als der 
des O,, q’ ist gemessen = 2,35 - 10-12 erg; Ap = 1,5 -10?. Mit dem Molekiilradius des 
0,: 2a ~ 3 - 10-8 und E/k = 1,1 - 10 würde man bekommen: A, ~ 6,4 - 10" (=) 
Angesichts der empirischen Unsicherheit des A,-Wertes wird man besser auch 
A, als Parameter beniitzen. 


Wir bestimmen also E/k und A, aus 2 gemessenen Werten von u,. Ars den ~ 
Messungen von Lewis und v. Elbe sehen wir: ie 


für m* = 3,054; T, = 427° abs; T, = 1343° abs ist u, = 158; ets 
für m* = 1,016; 7, = 468° abs; 7’, = 2044° abs ist u,= 747. 00 


Aus den angegebenen Zahlen finden wir mit Hilfe von (102) E/k = 0,947 - 104° 
abs; A, = 1,52-10" (cm/sec)?. Diese zwei Zahlen entsprechen gut der theore- 
tischen Erwartung für E/k und A, und geben eine erste Bestätigung dafür, daß die 
Formel (102) die Verbrennungsgeschwindigkeit der O,-Verbrennung darstellen 
kann. Weiter müssen aber jetzt alle übrigen gemssenen u,-Werte mit den 2 eben 
bestimmten Parametern aus (102) berechnet werden können. Das Ergebnis zeigt 
die folgende Tabelle 1. Die eingeklammerten Werte sind zur Berechnung von E/k 
und A, verwendet worden, sagen also nichts über den Vergleich zwischen Beob- 
achtung und Rechnung aus. In der letzten Spalte stehen zum Vergleich die 4 von 
Lewis und v. Elbe berechneten Verbrennungsgeschwindigkeiten. 

Der theoretische Wert von A, wäre mit dem soeben gefundenen Wert 9,47 - 10° 
für E/k: A, = 4,9- 10", wenn für 2 a, wieder der sicher zu kletne Sauerstoffmolekül- 

4 A» 

durchmesser eingesetzt wird. Aus dem Verhältnis a = ne — 1,76 würde der 
Ozonmolekiildurchmesser 1,76mal größer herauskommen als der von O,, das ist 
wohl zu groß, aber eine überraschend gute Übereinstimmung in der Größenord- 
nung. Die gemessene Aktivierungsenergie der Reaktion 20,—30, beträgt 
24 + 1Cal pro Mol; unser aus den Messungen mit Hilfe von (102) berechneter 
Wert ist 19,0 Cal, wieder ganz gute Übereinstimmung; es ist zu bedenken, daß 
unser Wert einiges kleiner sein muß als 24 Cal. Doch ist auf Quantitatives kein zu 
großer Wert zu legen bei den nicht geringen Vernachlässigungen unseres Ansatzes 
und unserer Rechnung. Wesentlich ist vielmehr, daß (102) allen qualitativen An- 
forderungen genügt, einigermaßen auch den quantitativen, und daß für E/k und 


4, sich vernünftige Werte ergeben. 
15* 
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Tabelle 1 


| 
a. henh. | %ı ber Fehler = u, von 
me | | | Gi. | | 7, | w((ber.) —| Lewis |Fehler 
8 | ne (102 uy (beob.) | berechnet Lewis 
| | cm in % von : 1m % 
| 
3,054 | 300 | 1239| 550 | 71,7] 31,6 | 7,65 30 | 253 | 360 
327 1259| 74,3 85,3| 29,0 | 7,54 15 
361 | 1289 | 100,0 | 108 | 26,2 | 7,35 > 
401 | 1322 | 134,2 | 136 | 23,8 | 7,16 1 
427 | 1343 | 158,0 |(158) | 22,2 | 7,05 (0) 451 | 185 
451 | 1362 | 178,0 | 191 | 20,1 | 6,95 7 
1,497 | 301 | 1680 | 141 238 | 31,5 | 5,65 69 
¥ 307 | 1685 | 187 249 | 30,9 | 5,63 33 
Lt 333 | 1704 | 224 284 | 28,5 | 5,56 27 
ee 376 | 1738 | 282 361 | 25,2 | 5,45 24 
2 408 | 1763 | 352 405 | 23,2 | 5,38 15 
441 | 1790 | 435 467 | 21,5 | 5,30 7 
477 | 1818 | 563 539 | 19,9 | 5,21 —4 
1,016 302 | 1922 | 160 369 | 31,4 | 4,95 130 333 | 108 
312 | 1928 | 287 394 | 30,4 | 4,91 37 
350 | 1958 | 379 465 | 27,1 | 4,85 23 
376 | 1976 | 495 520 | 25,2 | 4,80 5 
cf 411 | 2003 | 615 609 | 23,0 | 4,73 — | 
453 | 2035 | 733 710 | 20,9 | 466 | —3 | 
468 | 2044 747 (747) | 20,3 | 4,63 | (0) 664 —11,1 
486 | 2063 | 780 806 | 19,5 | 4,59 +3 


Nach (102) hängt die Verbrennungsgeschwindigkeit w, nicht vom Druck ab, 
* 

= a ist nach (49), (70), (71) gleich M,/M,, also von p, nicht abhängig. In 

Tabelle 1 sind die gemessenen Druckwerte nicht angegeben; von Lewis und 

v. Elbe sind die folgenden Drucke gemessen (in mm Hg): 


denn’ 


8,054 427 | 1848/2560 
- 1,016 302 | 1922 | 495 
1,016 468 | 2044 | 3760 


Man sieht, daß p, sehr stark variiert; soweit man die Gl. (102) als durch die 
Messungen gemäß Tabelle 1 bestätigt ansehen kann, wird man auch sagen können, 
daß die Druckunabhängigkeit der Verbrennungsgeschwindigkeit bestätigt ist. 


Aus (77a) folgt wegen (71), (49) 


Cor Pp; +1) (M,m* + (105) 


in unseren Rechnungen haben wir c,, als Konstante behandelt, was es nach den 
Definitionen (66), (67) nicht genau ist. Wir können aber die Güte dieser Näherungs- 
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annahme an der Tabelle 1 prüfen. Nach (105) muß 7',— 7’, bei festem m* konstant 
sein, wenn ¢,, wirklich konstant ist. Dies ist recht gut erfüllt, wie Tabelle 3 zeigt. 

Die geringe Abnahme von 7,— 7, mit zunehmendem 7’, bei festem m* spiegelt 
die Temperaturzunahme von ¢,, wider. Für 32c,,/R erhält man die Werte der 
4. Spalte; für 5 Freiheitsgrade läßt die kinetische Gastheorie den Wert 3,5 erwarten. 
Der höhere Wert » 4 entspricht der An- 


regung weiterer Freiheitsgrade durch die Tabelle 8 
Man kann die Verbrennungsgeschwin- | 7,—T, 
digkeit auch durch graphische Integration 
der Differentialgleichung (80) bestimmen. 3,054 300 | 908 3,98 
Wie in $8 gezeigt wurde, hat man die 401 921 4,05 
| 451 911 4,09 
Integralkurve zu suchen, welche vom 
Punkt Q (v,, ausgeht und in Pw=0, | | 
y= 0) mündet. Je nach dem Wert von 477 1341 493 
u, wird die von Q ausgehende Integral- 1,016 | 302 1620 418 
kurve links oder rechts oberhalb von P 15692 | 
vorbeilaufen, nur für einen ganz be- | 486 yn 


stimmten Wert von %, mündet sie in P. 


Die Rechnungen dieser Arbeit waren urspriinglich mit einer Formel fiir die StoBzahl Z 
gemacht, nach welcher beim StoB die gesamte Relativenergie tibertragen werden kann, 
nicht nur der Teil ¢,, der Bewegung in der Verbindungslinie der zwei Molekülmittelpunkte. 
Die Formeln für diesen Fall unterscheiden sich in «, von denen, die wir hier benutzen, um 


einen zusätzlichen Faktor Yr + 1, und es war bei y = 1 angepaßt worden, statt bei 
y = T,/3, vgl. (96). Für diesen Fall war durch graphische Integration gezeigt worden, daß 
die u,-Werte nach dem analytischen Näherungsverfahren und dem graphischen Verfahren 
sehr gut übereinstintmen; maximale Abweichung war 2%; E/k und A, waren wieder als 
Parameter behandelt. Wir können also mit Sicherheit behaupten, daß die graphische Be- 
stimmung von a, auch bei unserem Ansatz für Z mit großer Genauigkeit dieselben Werte 
geben wird, wie unsere analytischen Näherungsformeln. Damit sind zugleich unsere 
Näherungsformeln erneut als mathematisch gut erwiesen. 


Die rechte Seite ist als Funktion von y, daher auch von 7’ bekannt, wenn » (y) 
durch graphische oder näherungsweise analytische Integration nach den in $9, 
10 beschriebenen Verfahren bestimmt ist. 


Wir rechnen mit der Näherungsformel (95) für v. Beachten wir, daß y, =r,—r, 
ist, so wird: 


und wegen T = und t=1t, + y: 
d 


Ei: + =) Ag 


> 


§ 11. Der Temperaturverlauf in der Flamme a> 
Aus (92) folgt fiir den Temperaturanstieg wegen p,, = naeh u (78a) er 
(m +1) B(1—) Ve 


a x 
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1 
D, 
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| 
t. 
)) 
n 


Nach (108) ist der Anstieg der Temperatur wesentlich bestimmt durch e-v = e #7 - const, 
weil dieser Ausdruck sich in dy/dx besonders stark ändert. Wir setzen x = 0 für 


für 


T, — Ty) 27,7 

TG, rt Tı Tr 

ral 2 


und erhalten als Gleichung zur Berechnung des Temperaturverlaufes in der Flamme: 


Y 
J dy PıT, (110) 
(tr 


(m* + 1) (1—=)% 


Das Integral werten wir näherungsweise aus, indem wir setzen: re oe 


1 
damit ist das Verhalten des Ausdrucks — für y> 1 und für y < 1 richtig erfaßt; im 


ganzen Bereich der y- und y,-Werte geht der Unterschied zwischen der Näherungsformel 


— (110a) und dem genauen Wert von n nicht über 25% des genauen Wertes hinaus. 
J 


Wir setzen: 
, g pt 
eund erhalten 


1 1 Ya 2 
m=) + Ya-Vr yı —2 
E 3% 3 in(! t+ Vt, Va = 
Yyı VT Vt, Yr- Yr 


dabei ist y=r—r,; r=E/(kT). Wir haben damit x in Abhängigkeit von T, 
also den in der Flamme, 


(112) 


Für eine Abschätzung der Dicke der Flamme setzen wir: 


2 T,=T,—T,= T,(1- 7): T,=7,(1+ 7), (113) 


und betrachten den Abstand zwischen den zugehérigen Koordinaten x, und 2, 
als Maß für die Flammendicke. Mit der Abkürzung: 


6= 7 ,/T, = (114a) 
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K. Bechert : Zur Theorie der Verbrennungsgeschwindigkeit 


Nach einer einfachen Rechnung bekommt man: 
73 
1 
1+ 1 1 
| 1-6 


( 10+) Yan] 


n 

Für die Verbrennung mit den Daten: m* = 1,016; 7’, = 302; 7, = 1922; p, = 
495 mm Hg; w, (ber.) = 369, erhalten wir aus (111) und (115): K ~ 9,0 und als 
Flammendicke: 2,— x, & 3,7-10-% em. Zur Feststellung der Größenordnung 
genügt die Formel: 

2 41 4| T, 

Kr ly 6 of T, ( ) 

An Stelle von 3,7 - 10”? cm gibt diese Formel 3,9 - 10”? cm für das soeben gegebene 

Beispiel. 

Bei der Verwendung der Gln. (115) oder (115a) ist zu beachten, daß sich nach 
unserer Definition (113) die Temperaturen 7',, 7, um etwa 20% von 7',, 7’, unter- 
scheiden. Setzen wir 

ö 


10 


an Stelle von 6, so liegt 7’, um etwa 2%, unterhalb 7',, und 7’, etwa 2%, über 7’. 
Die Formeln (115), (115a) gelten unverändert auch für 0’. Mit dieser De Gaidion 
wird die Breite der’Flammenfront für unser Beispiel nach (115a), wenn wir sie 
zum Unterschied jetzt mit x, bezeichnen: x, — x, = 1,2 10? cm. 

Wegen des Faktors 1/K ist die Flammendicke gemäß (111) umgekehrt pro- 
portional zum Druck p, im Unverbrannten. Für diesen Schluß ist die Druckunab- 
hängigkeit der Verbrennungsgeschwindigkeit gemäß (102) entscheidend, daher 
auch der Ansatz Z ~ nj, der aus unserer Annahme der Reaktion zwischen gleichen 
Teilchen folgte. Bei anderer Druckabhängigkeit der Stoßzahl Z wird die Ver- 
brennungsgeschwindigkeit druckabhängig gemäß (103f, h), und die Flammen- 
dicke nach (111), (115a) anders vom Druck abhängig als gemäß 1/p,. 

Wie für die Verbrennungsgeschwindigkeit, läßt sich auch für die Flammen- 
dicke zeigen, daß unsere Gleichungen ihrer dimensionellen Form nach allgemein- 
gültig sind, wenn wir den Ansatz für die Reaktionszahl Z offen lassen. Die summa- 
rische Behandlung der Diffusion, die für die Durchführbarkeit unserer Rechnung 
und die Bestimmung des dimensionslosen Faktors in (115), (115a) wesentlich ist, 
beeinflußt nicht die dimensionelle Form der Gln. (115), (115a), und daher auch 
nicht die qualitative Abhängigkeit von den physikalischen Größen. Nach (1031, m) 
und (77a), (71) ist nämlich “dimensionsmiBig: 
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~ 


Andererseits ist nach (111): \ et 
ı (1—F) m + 1) 
K t/a Py Ty (1165) 


was mit (116a) übereinstimmt bis auf den Umstand, daß in (116a) der durch Di- 
steht, während 


Vz 

(116b) 7 zeigt. Unsere Formel (115a) fiir die Dicke der Flammenfront kann somit 

in ihrer dimensionellen Abhängigkeit von den physikalischen Größen als gesichert 


angesehen werden. Die velle Berücksichtigung der Diffusion würde diese Abhängig- 
keit nicht anders ergeben, als (115a) behauptet; uw, ist allgemein gegeben durch 
(103d,f). Spezielle Annahmen über die Reaktion wirken nur auf die Reaktions- 
zahl Z, die in uw, in der Form Vz vorkommt. Dagegen ist für das quantitative Er- 
gebnis die Form der geschweiften Klammer in (115a) von Bedeutung, die für unser 
Zahlenbeispiel etwa den Wert 10 hat, wenn wir die Breite x,— a; berechnen wollen. 


mensionsbetrachtungen nicht näher zu bestimmende Faktor 


Den Temperaturverlauf kann man auch graphisch aus (106) gewinnen, die 
rechte Seite ist eine bekannte Funktion von 7. Man hat also = { (7) und 


findet zu jedem dT das zugehörige dx = dT'/f (7'), so läßt sich die Kurve 7’ = T(z) 
graphisch bestimmen. Zu jedem Wertesatz 7',, 7',, m*, p, gibt es einen eigenen 
Temperaturverlauf. Qualitativ ist es aber immer derselbe: praktisch unmerkliche 
Steigung im Unverbrannten, steiler Anstieg im Verbrennungsgebiet, Übergang 
zur räumlich konstanten Temperatur 7’, im Verbrannten. Je höher die Verbren- 
nungsgeschwindigkeit und die Temperaturen im Verbrannten und Unverbrannten, 
desto kleiner die Flammendicke. 


Nach dem Verfahren dieser Arbeit lassen sich grundsätzlich alle Verbrennungen 
behandeln, bei denen es wesentlich nur auf den Stoß zwischen zwei Molekülarten 
ankommt. Der größere Teil der Ergebnisse dieser Arbeit setzt gleiche reagierende 
Moleküle voraus. Für ungleiche Reaktionspartner wird das Rechenschema kom- 
plizierter. Dieser Fall wird in einer weiteren Arbeit behandelt und auf die Ver- 
brennung von Kohlenwasserstoffen angewendet werden. Auch für diesen Fall läßt 
sich ein analytische Formel für die Verbrennungsgeschwindigkeit aufstellen. 


ee Institut fiir theoretische Physik der Universitit. ey 


E (Bei der Redaktion eingegangen am 13. September 1948.) Pate 
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Die Fluoreszenzanregung von festem und flüssigem Naphthalin, 
Diphenyl und Phenanthren durch Alphateilchen, 


it schnelle Elektronen und Gammastrahlen u ” 
rt Vom Lieselott Herforth und Hartmut Kallmann ee 
> 
w. In folgenden ei Versuche über die Anregung und die Lichtausbeute von 
n. 


festem, kristallinem Naphthalin, Diphenyl und Phenanthren durch Alphateilchen, j 
ie schnelle Elektronen und Gammastrahlen beschrieben. Ihre physikalische und 
technische Lichtausbeute relativ zu ZnS—Ag werden angegeben. Die Absolut- 
ausbeuten betragen größenordnungsmäßig: für x-Teilchen 0,3%, 1,0% und 0,6% ; 


r) für y-Strahlen 5%, 7,5% und 11%. Die Abhängigkeit der Lichtausbeute von 
n der Temperatur wird bis zum Schmelzpunkt herauf bestimmt. Beim Schmelzen 
le geht die Fluoreszenzausbeute auf einen verschwindend kleinen Bruchteil zuriick; 
ig beim Erstarren der Schmelze geht sie bei allen drei Substanzen wieder in die Höhe, 


L- ‚erreicht aber nur etwa den halben ursprünglichen Wert. Pulverisieren der er- 
Ly starrten Substanz bringt die physikalische Ausbeute wieder auf den alten Wert 
zurück. Die Gesamtausbeute sowie das Verhalten im flüssigen wie im wieder er- 


n starrten Zustand werden theoretisch diskutiert. N 

I- 1. Einführung 


Im Anschluß an die Fluoreszenzmessungen von kristallin-festem Naphthalin 
durch Korpuskularstrahlanregung!) haben wir die Anregung von Diphenyl und 
Phenanthren durch «-Teilchen, schnelle Elektronen und y-Strahlen untersucht 
und die physikalischen und technischen Ausbeuten der drei organischen Stoffe 
zunächst in kristallin-festem Zustand miteinander verglichen. Es wurde mit der 
gleichen Elektronenvervielfacheranordnung gemessen, wie wir sie in allen voran- 
gehenden Arbeiten) verwendet haben und wie sie die Abb. 1 weiter unten wieder- 
gibt. Genau wie bei Naphthalin zeigten auch verschiedene Diphenyl- und Phenan- 
thren-Sorten verschieden starke Lichtausbeuten. Reinigung durch Umkristalli- 
sieren in Methanol verbesserte die Lichtausbeute in jedem Falle. Dabei erwies 
sich das Umkristallisieren in Methanol für die Fluoreszenz günstiger als in Tetra- 
chlorkohlenstoff. Anregungsversuche an Naphthalin, welches in Tetrachlor- 
kohlenstoff umkristallisiert und sublimiert war, zeigte für Leuchtstoffschichten, 
> 200 mg/em?, eine sogar schlechtere technische Ausbeute als nicht umkristalli- 
siertes, während für die physikalische Ausbeute die Art des bei der Herstellung 


1) I. Broser, L. Herforth, H. Kallmann und U. Martius, Z. Naturf. 8a, 6 (1948). 
2) I, Broser u. H. Kallmann, Z. Naturf. 2a, 439 (1947); 2a, 642 (1947); u.l. c.!). 
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verwendeten Lösungsmittels keine Rolle zu spielen scheint. Dies zeigt, daß durch 
Umkristallisieren in Tetrachlorkohlenstoff Spuren von Verunreinigungen ange- 
reichert werden, ‘die eine zusätzliche Absorption oder Streuung des Fluoreszenz- 
lichtes im organischen Leuchtstoff bewirken aber nicht eine Verminderung seiner 
physikalischen Lichtausbeute. Auch alle benutzten kristallin-festen Diphenyl- 
schichten wurden, genau wie früher die Naphthalinschichten, durch Sublimation 
unter Druck von 5mm Hg gewonnen. Diese Art der Herstellung von Schichten 
erwies sich für Phenanthren als nicht zweckmäßig, da die sehr locker aneinander- 
sitzenden Kristalle zu sehr groBer Streuung fiihrten. Es erwies sich bei Phenanthren 
als besser, die in Methanol umkristallisierte Substanz in Schälchen zu füllen und 
durch Zusammendrücken gleichmäßige Schichten herzustellen. Die physikalische 
Ausbeute ist bei Verwendung von nur umkristallisierter Substanz genau so groß 
wie von sublimierten Schichten. Die 


Einstrahlung . technische Ausbeute wird durch die 

Leuehtstoff Glasplatte Art der Herstellung beeinflußt 
Fluoreszenzlich' Multiflex- (Streuung). 

Öalvanometer Wir erläutern noch einmal das 


Meßverfahren: Es wird sowohl mit 
x-Teilchen, die nur sehr wenig in die 
. Substanz eindringen, als auch mit 
Abb. 1. Elektronenvervielfacheranordnung einer harten y-Strahlung, die die ganze 
bei kontinuierlicher Einstrahlung Substanzmenge durchdringt und prak- 
tisch nur sehr wenig geschwächt wird, 
angeregt und die Intensität des Fluoreszenzlichtes als Funktion der Massen- 
belegung gemessen. Abb. 1 zeigt das Schema der Versuchsanordnung. Aus der 
2 Intensitätskurve bei Anregung mit «-Teilchen in Abhängigkeit von der Massen- 
A belegung ergibt sich dann die Streuung des Fluoreszenzlichtes, während aus der 
y-Strahlmessung auf eine etwa vorhandene wahre Absorption zu schließen ist. 
Abweichung der Intensitätskurve von der geraden Linie bei Anregung mit y- 
Strahlen bedeutet wahre Absorption, wobei die wirkliche, die Absorption verur- 
sachende Schichtdicke (wahrer Lichtweg) infolge der Streuung wesentlich größer 
ist als die gemessene Schichtdicke®). 


Ry Ry Rg Ry 


+ 
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2. Die Fluoreszenzanregung durch Alphateilchen 


Unsere früheren Anregungsversuche von Naphthalin®) waren mit einem a- 
und f-Teilchen aussendenden Radiumpräparat ausgeführt worden, und der An- 
; teil der allein durch «-Teilchen erregten Fluoreszenz muBte mit Hilfe von Ab- 
Bes. sorptionsmessungen ermittelt werden. Jetzt wurden alle Messungen mit einem 
reinen Alphastrahler, etwa 1 mC Polonium, durchgeführt. Es wurden die An- 
regungsversuche von sublimiertem Naphthalin wiederholt und die früheren Mes- 
sungen bestätigt. Die Lichtstreuung im Naphthalin war so gering, daß sie noch 
innerhalb der Fehlergrenzen der Meßgenauigkeit lag. In Abb. 2 ist die Intensität 
des Fluoreszenzlichtes bei «-Teilchenanregung in willkürlichen Einheiten für 
_ Naphthalin, Diphenyl und Phenanthren in Abhängigkeit von der Schichtdicke 
in mg/cm? aufgetragen. Die Ausbeute steigt bei Naphthalin bei Schichtdicken, die 
der Reichweite der «-Teilchen entsprechen, auf einen Maximalwert an, der dann 


3) ].c. unter ?). 
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bis zu hohen Schichtdicken annähernd bestehen bleibt. Dieser große Bereich 
konstanter Lichtintensitäten zeigt (da immer von der der Photozelle abgewandten 
Seite bestrahlt wurde, s. Abb. 1), daß eine Lichtstreuung im Naphthalin nicht vor- 
handen ist, während bei den beiden anderen Substanzen die Fluoreszenz-Inten- 
sitäten bei x-Teilchenerregung nach Erreichen eines Maximums infolge der Streu- 
ung des Lichtes in dem mikro- 
kristallinen Material sehr er- 
heblich absinken. Naphthalin 
zeichnet sich vor allen anderen 
Substanzen also dadurch aus, 
daß es fast keine Streuung auf- 
weist, obwohl der Napthalin- 
block an sich trübe aussieht. 
Diese Trübung scheint aber nur 
an der äußersten Oberflächen- 
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\ Sp 
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Intensität [willkürlich Einheiten] 


yer! 


schicht merklich zu sein. Eine 0 700 200 300 400 500 [mg fem?) 
mit der Schichtdicke anwach- ——> Schichtdicke 
sende Absorption ist auch fiir Abb. 2. Fluoreszenzstrahlung bei «-Teilchen- 


Ti - or. anregung von kristallin-festem Naphthalin, Diphe- 
nyl und Phenanthren in Abhängigkeit von der 
anden. 1p eny un enan- Schichtdicke 


thren verhalten sich infolge der 

Streuung in der technischen Ausbeute also wesentlich ungiinstiger als Naphthalin. 
Thre physikalische Ausbeute, die durch den Anstieg der Kurve im Nullpunkt bzw. 
die Höhe des Maximums gegeben ist, ist bei «-Teilchenerregung aber größer als 
fiir Naphthalin. 

Beim Vergleich der Ausbeutekurven ist zu berücksichtigen, daß die verwendete 
Photozelle für die Fluoreszenzstrahlung wegen der verschiedenen Wellenbereiche 
verschieden empfindlich ist. Während Naphthalin und Diphenyl nahezu im gleichen 
Wellenbereich (um 3850 Ä) emittieren, emittiert Phenanthren in einem etwas 
langwelligeren Gebiet (um 4500 A). Da die Photozellen-Empfindlichkeit für dieses 
Gebiet um den Faktor 1,9 geringer war, mußte die gemessene Phenanthren-Kurve 
(in der Abb. 2 punktiert gezeichnet) um diesen Faktor vergrößert werden, um mit 
den anderen Substanzen verglichen werden zu können. 

An Naphthalin und Diphenyl wurden auch einzelne Lichtblitze gemessen®). 
Die Größe der Lichtblitze betrug bei Naphthalin etwa !/,, und bei Diphenyl !/, der 


> Größe der Szintillationen in 
* ZnS—Ag. Im Gegensatz zu 
my! ,| ZmS—Ag ist aber die Abkling- 
X 60 1 zeit außerordentlich kurz und 
Q 4 sicher < sec. 

3. Die Fluoreszenzanregung 
I gemessen. durch schnelle Elektronen 

R Die Anregung durch schnelle 


m 20 30 90 500 Kiektronen erfolgte wie früher 


—— Schichtdicke . . 
Abb. 3. Fluoreszenzstrahlung bei Betateilchen- 
anregung von kristallin-festem Naphthalin, Dipheny] 5) Nach Messungen von I, 
und Phenanthren in Abhängigkeit von der Schicht- Broser u. H. Kallmann, über 
dicke die gesondert berichtet wird. 
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_ In Abb. 4 ist die Fluoreszenzerregung durch y-Strahlen des Diphenyls und 
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schwachen x- und f-Teilchen aussendenden Radiumpräparat. Zur Ausschaltung 
der x-Teilchen wurde zwischen Präparat und Leuchtstoff eine 5,2 mg/cm? 
dicke Papierfolie gelegt. In Abb. 3 ist die Fluoreszenzstrahlung für ver- 
schiedene Schichtdicken für Diphenyl und Phenanthren aufgetragen. Zum 
Vergleich ist wieder die entsprechende Kurve für Naphthalin eingezeichnet, 
Für schnelle Elektronen ist die physikalische Ausbeute von Phenanthren (auf 
gleiche Empfindlichkeit korrigierte Kurve) am besten und die von Diphenyl 
1,5 mal besser als die von Naphthalin. Jedoch schon für Schichten > 300 mg/cm? 
übertrifft die technische Ausbeute ‚von Naphthalin diejenige von Diphenyl. Die 
technische Ausbeute von Phenanthren wird schon für Schichten > 150 mg/cm? 
sehr schlecht. Dieses Verhalten ist wieder durch die in Diphenyl und besonders 
stark in Phenanthren auftretende Lichtstreuung bedingt. 


4. Die Fluoreszenzanregung durch Gammastrahlen 


ra 


Phenanthrens und zum Vergleich die des Naphthalins in Abhängigkeit von der 
Schichtdicke aufgetragen. Während die Naphthalinkurve bis zu Schichtdicken 


ae von 500 mg/cm? und wie aus 
90 anderen Messungen hervorgeht, 
auch noch darüber hinaus ein 
2 60 s lineares Ansteigen aufweist, 
= ols zeigt der Verlauf der Diphenyl- 
NS F .| Kurve, daß Diphenyl schon für 
| Prenantiven Schichten von etwa 300 mg/cm? 
S die eigene Fluoreszenzstrahlung 
& | | abserbiert. Bei Phenanthren 


0 0 200 300 400  500(mgfm* tritt diese Absorption schon 

———> Schichtdicke bei sehr diinnen Schichten ein. 
Abb. 4. Fluoreszenzstrahlung bei Gammastrahlen- Diese Absorption macht sich 
anregung von kristallin-festem Napthalin, bei diesen Substanzen nur des 


Dipl lundP thren in Abhängigkeit d : 3 


bei ihnen der wahre Lichtweg 
durch Streuung auBerordentlich stark vergréBert wird (siehe weiter unten 
Abb. 10). 

Auch bei Diphenyl wurden die durch einzelne schnelle Elektronen ausgelösten 
Lichtblitze gemessen. Es zeigte sich entsprechend der Messung bei kontinuier- 
licher Erregung, daß die Intensität der einzelnen Lichtblitze größer ist als in Naph- 
thalin, solange man nicht mit zu energiereichen Elektronen arbeitet. Für energie- 
reiche Elektronen, z.B. Höhenstrahlen, bedeutet daher die Verwendung von 
Diphenyl auch bei Messung einzelner Lichtblitze keine Verbesserung gegenüber 
Naphthalin. 


5. Messung der spektralen Verteilung des Fluoreszenzlichtes 


Durch Einschalten verschiedener Jenaer Farbfilter in den Strahlengang des 
Fluoreszenzlichtes wurde die ungefähre Lage der Fluoreszenzspektren der Leucht- 
stoffe ausgemessen. Anregung durch «-, ß- und y-Strahlen führte zu gleichen 
Fluoreszenzspektren, vorausgesetzt daß die verwendete Schicht so dünn war, daß 
weder Streuung noch Absorption auftreten konnte. Naphthalin und Diphenyl 
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zeigten fiir diinne Schichten gleiche Spektren. Fiir dickere Schichten wird 
fir Diphenyl der Anteil kurzwelligen Lichtes am Gesamtspektrum immer 
geringer. In Abb. 5 ist fiir 


Naphthalin und Diphenyl 

der Anteil der Fluoreszenz- & 40 

strahlung mit einer Wellen- & Lal 

linge >4150Ä relativ zur ed 

gesamten Fluoreszenzinten- ~2 2 

sitit in Abhängigkeit von 

Für Naphthalin bleibt dieses s! 

Verhältnis konstant, während & 


es für Diphenyl infolge der Schichtdicke 


stärkeren Absorption des Abb. 5. Relative gefilterte Fluoreszenzstrahlung 

kurzwelligen Anteils mit wach- nei Alphateilchenanregung von kristallin-festem 

sender Schichtdicke ansteigt. Naphthalin und Diphenyl in Abhängigkeit von der 
Schichtdicke 


6. Ausbeute-Messungen 
In der unten folgenden Tabelle sind die physikalische und die technische Aus- 
beute für die drei Substanzen bei Erregung mit a-Teilchen, Elektronen und y- 
Strahlen angegeben. In den beiden letzten Spalten sind außerdem noch die Ab- 
solutwerte angeführt, die einer demnächst erscheinenden Arbeit*) entnommen 
sind. Auf die Diskussion dieser Werte kommen wir im letzten Absatz ausführlich 
zu sprechen. 


Tabelle: Physikalische, technische und absolute Ausbeuten 


Emissions- physikalisch | technisch absolut 
| x ß y x B y : y 
ZuS-Ag....| 4500 | 100 100 100 100 | 100 | 100. 0,280! 0,135 
Naphthalin. . . | 3850 1,1} 37 37 1,4 233 1000 0,003) 0,05 
Diphenyl. . . . | 3850 3,6| 55 55 4,2| 233 580 0,010) 0,075 
Phenanthren . . 4500 2,1| 104 81 2,5! 198 241 | 0,006) 0,21 


7. Abhängigkeit der Fluoreszenz von Temperatur und Aggregatzustand 
Um die Abhängigkeit der Lichterregung der untersuchten Substanzen auch 
von Temperatur und Aggregatzustand zu bestimmen, wurde eine Versuchsanord- 
nung, wie sie Abb. 6 schematisch 


3-Emstrahlung zeigt, verwendet. Dicht über dem 


dinner | | | | | | | Elektronenvervielfacher wurde in 
Papdechel Ihermoelernen! einem Glasschälchen die gereinigte 
Substanz in einen elektrischen Ofen 
NZZ Maphthalın in Glasschale eingebracht und von oben mit y- 


Huoreszenz/icht 


Strahlung erregt. Die geringe Be- 


einflussung des Elektronenverviel- 
Elek/ronenvervielfacher fachers durch das Erhitzen wurde be- 

Abb. 6. Versuchsanordnung zur Messung der 6) I. Broser, H. Kallmann u. 
Flucreszenz hei verschiedenen Temperaturen U. Martius, Z. Naturforschg. (z. Z. 


des Leuchtstoffes im Druck). 
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relative Intensitat [in Yo) 


sonders bestimmt und bei den Messungen in Abzug gebracht. Eine MeBkurve mit 
Naphthalin zeigt Abb. 7, in der die relative Intensität bei der Temperatur ¢ zu 
der bei 20° C aufgetragen wurde. Man sieht, daß etwa von 40°C an ein etwas 
stärkeres Absinken der Fluoreszenz auftritt und erst in der unmittelbaren Nähe 
des Schmelzpunktes die Intensität außerordentlich stark (bei Naphthalin prak- 


tisch auf Null) zurückgeht”). Bei 


V9? tieferen Temperaturen steigt die 
700 nn Fluoreszenzausbeute nur noch 

3 wenig an und ist bei etwa—30°C 

So = um etwa 14% größer als bei 

Be R Zimmertemperatur. Gleiche Resul- 
af ® a tate wurden mit Diphenyl und 
2 90 De Phenanthren erzielt. Auch bei 
| diesen Substanzen sank die Fluo- 
Sw x reszenz-Intensität in der unmittel- 
| TN baren Nähe des Schmelzpunktes 

120 10 700] auf sehr kleine Werte herab. Wir 


führten sodann mit der flüssigen 
und der wieder erstarrten Substanz 
Messungen durch, und zwar haben 
wir die Fluoreszenz-Intensität in 
Abhängigkeit von der Schichtdicke 
aufgetragen (für alle drei Substanzen Abb. 8, 9 und 10). Kurve I gibt immer die 
Intensität der Ausgangssubstanz an, die durch Auskristallisieren bzw. Subli- 
mieren gewonnen worden war. Kurve III gibt die Intensität der flüssigen Sub- 
stanz und Kurve II diejenige der aus der Schmelze wieder erstarrten Substanz 
an. Es zeigte sich, daß bei allen drei Stoffen die Intensität der Substanz im 
flüssigen Zustand gegenüber der 


——— /emperaltur 

Abb. 7. Relative Intensität der Fluoreszenz von 
kristallin-festem Naphthalin in Abhängigkeit von 
der Temperatur 


Intensität der festen Substanz T 90 Th 
außerordentlich klein war (bei & 

Naphthalin sogar Null). Hier- „7? 

bei sei erwähnt, daß das geringe RS 60 ee 

Leuchten von fliissigem Di- - 
phenyl und Phenanthren in I 0 
demselben Wellenlängengebiet 

liegt wie das Fluoreszenzlicht 

der entsprechenden festen Sub- X — a a 7 
stanzen. Es handelt sich also ond 


ine Fl —— Schichtdicke 

Abb. 8. Fluoreszenz-Intensität in Abhängigkeit von 
zenzanregung vonnderFiuSSI&- der Schichtdicke für aufsublimiertes (la), zer- 
keit vorhandenen Verunreini- mörsertes und auskristallisiertes (Ib), erstarrtes (II) 
gungen. Die physikalische Aus- und flüssiges (III) Naphthalin bei y-Strahlenanregung 
heute der erstarrten Substanz - 
lag bei allen drei Körpern unterhalb (etwa um den Faktor !/, kleiner) der der 
Ausgangssubstanz. Sublimiertes Naphthalin gab mithin die doppelte Fluoreszenz- 


7) Vergleiche Versuche bei optischer Anregung von festen, flüssigen und erstarrten 
Kohlenwasserstoffen und ihren festen Lösungen mit Naphthacen, die zu ähnlichen Er- 
gebnissen führten, von A. Miller u. C. A. Baumann, Canc. Res. 8, 217 (1943). Siehe 
auch Hausser u. Mitarb., Z. physik. Chem. (B) 29, 417 (1935), 
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Intensität als aus der Schmelze erstarrtes Naphthalin. Bei Phenanthren wies 
die erstarrte Substanz eine wesentlich größere Durchsichtigkeit, d.h. eine sehr 
viel kleinere Streuung auf als die aus der auskristallisierten Substanz herge- 
stellten Schichten. So ist der Umstand zu erklären, daß die Kurve II der Abb. 10 
die Kurve I bei etwa 200 mg/cm? Substanz überschneidet. Kurve I weist die 
größere physikalische Licht- 
ausbeute auf. Aber schon bei 
Schichtdicken von 200 mg/cm? 
ist infolge der sehr starken 


Ss 


S 
N 


Streuung der Schicht eine so 
starke Absorption vorhanden, 


Intensität [milkürhche Einheiten] 


daß nahezu eine Sättigung der 0 7: 
Lichtintensitat eintritt. Phen- 

anthren hat gegenüber Naph- 2 

thalin den Nachteil, daß man tr z Ä 
bei ihm keine durchsichtigen 0 m 200300 


Schichtdicke 

Abb. 9. Fluoreszenz-Intensität in Al-hängigkeit von 
der Schichtdicke für aufsublimiertes (I), erstarrtes (II) 
und flüssiges (III) Diphenyl bei y-Strahlenanregung 


kristallinen Schichten erhält, 
während das bei Naphthalin 
durch Sublimieren der Fall ist. 

Wir haben nun versucht 
festzustellen, worauf diese bei allen drei Substanzen auftretende Verringerung 
der physikalischen Lichtausbeute bei der erstarrten Substanz beruht. Um zu- 
nächst sicher zu sein, daß diese geringere physikalische Lichtausbeute nicht 
durch Inhomogenitäten in den aus der Schmelze erstarrten Schichten hervor- 
gerufen wird, haben wir, obwohl eine solche Annahme schon durch das lineare, 
regelmäßige Ansteigen der Fiuores- 


40 zenz-Intensität der Kurve II in Abb. 8 
NR a I und 10 mit wachsender Massen- 
> y ’ belegung widerlegt wird, Schichten 
von 50, 100 und 150 mg/cm? Massen- 
x Pe belegung aus der erstarrten “Masse 
70 hergestellt und sie mit x-Strahlen an- 
S 4 geregt. Bei x-Strahlenerregung wissen 
R | wir aus den Naphthalin-Messungen 
8 . (siehe die Naphthalin-Kurve der 

: Abb. 2), daB eine Abhangigkeit von der 


700 


Massenbelegung überhaupt nicht vor- 
Schichtdicke 


handen ist. Es zeigte sich nun bei 
diesen Messungen mit x-Teilchen- 
erregung an den erstarrten Schichten, 
daß ebenfalls keine Abhängigkeit von 
der Schichtdicke (Massenbelegung) 
vorhanden war, und daß die Ausbeute nur etwa 50%, der Erregung bei sub- 
limiertem Naphthalin betrug. 

Weiter haben wir noch folgende zwei Versuche durchgeführt: Es wäre möglich, 
daß die Kristalle in den erstarrten Schichten in ganz bestimmter Weise an der Auf- 
lagefläche anwachsen, und daß sie in verschiedenen Richtungen verschieden stark 
emittieren. Wir haben daher an einzelnen erstarrten Blöcken die Lichtintensität 
in verschiedenen Richtungen gemessen, ohne nach irgendeiner Richtung einen 


300 
Imgjem“] 
Abb. 10. Fluoreszenz-Intensität in Abhängig- 
keit von der Schichtdicke für auskristalli- 
siertes (I), erstarrtes (II) und flüssiges (III) 

Phenanthren bei y-Strahlenanregung 
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merklichen Intensitätsunterschied feststellen zu können. Schließlich haben wir 
noch die erstarrten Schichten nach ihrer Messung fein gepulvert und das Kristall- 
pulver zu Schichten zusammengepreßt und dann gemessen. Die physikalische 
Lichtausbeute stieg hierbei auf den alten Wert an. In Abb. 8 sind diese Mes- 
sungen durch die punktierte Kurye Ib wiedergegeben, und zwar ergab sich dann 
genau dieselbe Kurve, gleichgültig ob wir die Messungen mit Naphthalin vor- 
nahmen, das aus dem sublimierten, festen Zustand zu Pulver zerrieben war, oder 
aus Naphthalin, das aus der Schmelze erstarrt war. Zu der gleichen punktierten 
Kurve kommt man auch bei Verwendung von Schichten aus umkristallisierter 
Substanz. Daß die punktierte Kurve im Gegensatz zu der Kurve Ia nicht linear 
mit der Schichtdicke ansteigt, sondern eine Sättigung anzeigt, beruht darauf, daß 
sich die Absorption in der gepulverten sowie auskristallisierten Substanz wegen des 
durch die große Streuung sehr stark vergrößerten wahren Lichtweges schor: be- 
merkbar macht. 


Zusammenfassend kann man also sagen, daß die Herabminderung der physi- 
kalischen Lichtausbeute bei der aus der Schmelze erstarrten Substanz weder auf 
eine Orientierung der Kristalle oder eine bevorzugte Lichtaussendung in einer be- 
stimmten Richtung, noch auf irgendeine Beeinflussung durch Fremdsubstanzen 
zurückzuführen ist. Wir möchten noch besonders darauf hinweisen, daß die Aus- 
beuteverminderung auch nicht durch eine Verschiebung der Frequenz des Fluores- 
zenzlichtes in ein für die Photozelle weniger empfindliches Gebiet vorgetäuscht 
wird. Wir haben uns überzeugt, daß eine Frequenzverschiebung nicht stattfindet. 
Auf die vermutliche Deutung des Effekts kommen wir im nächsten Absatz zurück. 


8. Diskussion der Meßergebnisse 


- E- a) Anregungszustand der emittierten Elektronen 


Die von uns untersuchten Fluoreszenzspektren schließen sich unmittelbar, um 
einige 100 Ä nach Rot verschoben, an das jeweilige Absorptionsspektrum an (nicht 
nur bei den hier erwähnten Stoffen, sondern auch noch bei einer Reihe anderer 
Stoffe®), die konjugierte Doppelbindungen enthalten). Sie ähneln außerordentlich 
den von Haußer®) und Mitarbeitern bei einer großen Anzahl von Kettenmolekülen, 
die konjugierte Doppelbindungen enthalten, beobachteten Fluoreszenzspektren. 
Die dort gefundenen Fluoreszenzspektren sind, da sie im festen Zustand, in der 
Lösung und im Dampf beobachtet wurden, zweifellos den Molekülen der unter- 
suchten Substanz selber und nicht einer Beimengung zuzuschreiben. Wegen der 
außerordentlichen Ähnlichkeit der hier beobachteten Fluoreszenzspektren mit 
denen von HauBer, insbesondere wegen der engen Koppelung zwischen Emission 
und Absorption, besteht kein Zweifel, daß es sich auch bei den hier untersuchten 
Fluoreszenzspektren um die Fluoreszenz der Moleküle des Grundstoffes infolge 
ihrer Doppelbindungselektronen handelt. Alle unsere Messungen sprechen dafür, 
vor allem die Reinigung des Stoffes durch Umkristallisieren. Natürlich würden 
diese Messungen allein nicht genügen, um das Leuchten von Verunreinigungen aus- 
zuschließen, da unsere Reinigungsmethoden nicht die letzten chemischen Mög- 


* lichkeiten ausschöpfen; aber die außerordentliche Übereinstimmung mit den von 


Haußer und anderen beobachteten Spektren und der Umstand, daß solche 


8) Dissertation Lieselott Herforth (1948). 
®) Haußer, R. Kuhn u. E. Kuhn, Z. physik. Chem. (B) 29, 417 (1935). 
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Spektren auf Grund der Theorie iiber die Doppelbindungselektronen zu erwarten 
sind, bestätigen diesen Schluß. Wir betonen dies hier besonders, weil vielfach 
gegenteilige Behauptungen aufgestellt werden. 


Nach der Theorie rührt die Lichtemission von der Anregung der beiden 
lockerst gebundenen Elektronen her, die im Energieband der Doppelbindungs- 
elektronen !%) (1 Elektron pro Kohlenstoffatom) zusammengefaßt sind. Die energe- 
tische Lage dieser Elektronen ist in zahlreichen Untersuc hungen berechnet 
worden. Nur die Hälfte aller Terme dieses Bandes sind von Elektronen besetzt 
(und zwar doppelt besetzt); das Absorptionsspektrum entspricht einem Elektronen- 
übergang von dem obersten besetzten Zustand in den tiefsten unbesetzten Zustand 
dieses Energiebandes, wie es in Abb. 11, dem Termschema (schematisch) 
des untersten Energiebandes 
der Doppelbindungselektro- 
nen des Naphthalins, durch 
Pfeile angedeutet ist. 

Bei diesem Übergang 
werden außerdem noch 
schwingungen angeregt. Die —_Absorptiorsspektrum 
Emission erfolgt von einem gehörender Elektronen- 


übergan 
angeregten Elektronenzu- 
Abb. 11. Termschema des untersten Energichandes 
stand, bei dem die Energie 


2 : der Doppelbindungselektronen fiir Naphthalin 
der Kernschwingungen durch (schematisch) 

Wärmestöße bereits abge- 

führt ist. Innerhalb des Gebietes des Fluoreszenzspektrums ist gar keine oder 
nur eine äußerst schwache Absorption beobachtet worden. Die Rotverschiebung 
des Fluoreszenzspektrums, die bei allen Systemen, die konjugierte Doppel- 
bindungen enthalten — nicht. nur bei den Polyenen, sondern auch bei poly- 
eyelischen Systemen —, auftritt und die vielfach in keiner Beziehung zu den 
angeregten Kernschwingungen der C-Atome steht, wird nun folgendermaßen 
gedeutet: Zu dem oberen Elektronenzustand gehört außer einem Singlet-Zustand 
auch noch ein Triplet-Zustand. Es könnte nun entweder die Emission von dem 
Triplet-Zustand ausgehen, der in der Absorption wegen der geringeren Über- 
gangswahrscheinlichkeit Singlet-Triplet-Zustand praktisch nicht zu erscheinen 
braucht, dann müßte sie wegen der geringen Übergangswahrscheinlichkeit eine 
lange Abklingdauer haben, oder es könnte nach Mulliken diese Rotverschiebung 
des Emissionsspektrums als eine Anregung von Verbiegungsschwingungen des 
Moleküls mit relativ kleinem hy aufgefaßt werden, so daß also die Emission zu 
ihrem wesentlichen Teil in einen Zustand führen würde, der in der unangeregten 
Substanz nicht vorhanden ist. 


unbesetzte 
Zustände 


doppelt 
Zustände 


Innerhalb des Fluoreszenzspektrums ist bei uns nur eine verschwindend kleine 
Absorption (siehe besonders die Kurven Ia und II der Abb. 8, I und II der Abb. 9 
und I der Abb. 10) vorhanden. Es ist also in keinem Teil des Fluoreszenzspektrums 
eine irgendwie erhebliche Absorption festzustellen. Man kann nun in bekannter 
Weise aus dem Absorptionskoeffizienten die Oszillatorenstärke der absorbierenden 
Elektronen berechnen und daraus wieder die demselben Übergang zugeordnete 


1%) Hückel, Grundzüge der Theorie ungesättigter und aromatischer Verbindungen 
(1940), Verlag hemie, G.m.b.H., Berlin. 


Ann. Physik. 6. Folge, Bd.4 


A | | 
wir 4 
stall- a 
sche 4 
Mes- | 
lann 
vor- 
oder i 
Tten 
erter | 7 | 
near | 
daß | 
ı des 
be- 
1ysi- | 
auf | 
be- | 
nzen 
Aus- 
res- 
scht 
det. | 
ück. 4 | 
x d 
| | | 


DE Annalen der Physik. 6.Folge. Band4. 99 | Herfo 


‘ Lebensdauer des angeregten Zustandes. Tut man dies, so erhält man fiir das | alle 
Absorptionsspektrum eine Oszillatorenstärke von etwa ¥/; pro Doppelbindungs- f Mole! 
elektron. Bedenkt man ferner, daß in Wirklichkeit nicht alle n- Doppelbindungs- dere 
elektronen, sonder nur zwei, nämlich die beiden lockerst gebundenen, die optische | absie 


; Absorption verursachen, so würde sich daraus eine Übergangswahrscheinlichkeit | lichk 
Bor bzw. Lebensdauer des angergten Zustandes von etwa 108 sec ergeben. elekt 
Würde nun ein Teil des Emissionsspektrums Übergängen entsprechen, die in e 
den unangeregten, normalerweise vorhandenen Grundzustand führen, so könnten _— 
wir aus der Größe der gemessenen Absorption in diesem Spektralbereich wiederum “a 
auf die Ubergangswahrscheinlichkeit des Emissionsprozesses schließen. Führt B- 
man diese einfache Rechnung für unsere Messungen durch, so ergibt sich unter Dop} 
Zugrundelegung der geringen Absorption in dem Emissionsgebiet eine Lebens- Nun , 
dauer der angeregten Zustände von etwa 10-* sec. Nun kennen wir aus unseren Basti 
Messungen an einzelnen Lichtblitzen in diesen Substanzen eine obere Grenze dieser En 
Lebensdauer. Die Lebensdauer ist sicherlich < 10 sec*). Es ist auch keine Spur rm , 
eines langwelligen Nachleuebtens vorhanden. Es muß also auch der kurzwellige er = 
Teil des Fluoreszenzspektrums offenbar in Zustände führen, die gegen den normaler- ß 
* weise schon vorhandenen Grundzustand angeregt sind. Anders sind die geringe a 
; Absorption und die kleine Lebensdauer im Emissionsgebiet nicht miteinander in ay 
Einklang zu bringen, wenn man an der Annahme, daß der Grundstoff das Licht na : 
aussendet, festhält. 
elekt 
ausb 
b) Anregungsmechanismus der Fluoreszenz Zusti 


\ 
Wie schon oben ausgeführt, beruht die Aussendung des Fluoreszenzspektrums | "8° 
auf einer Anregung der beiden lockerst gebundenen Doppelbindungselektronen haup 


; im Molekül. Nun hat die Messung der absoluten Ausbeute bei y-Strahlanregung Selbs 
a (siehe die Tabelle in Abschnitt 6) ergeben, daß Phenanthren eine Ausbeute von Elek 
‘ über 10%, aufweist, d.h. 10% der gesamten bei der Bestrahlung mit y-Strahlen dazu 


an Phenanthren abgegebenen. Energie wird im Gebiet der Emissionsbande mit | ?ust 
einem Maximum bei 4500 A als Licht ausgestrahlt. Diese Ausbeute ist nur etwa bind 
um den Faktor 3 geringer als die Lichtausbeute des besten ZnS bei Anregung durch | 3eset 
a-Strahlen. Will man die Größe dieser Ausbeute verstehen, so muß man folgendes 5%, 
bedenken: Es handelt sich bei unseren organischen Substanzen um Moleküle, | !n de 
die nur relativ wenig von den Nachbarmolekülen gestört werden. In der eigenen Aush 
festen Substanz, in den Lösungen und im Gaszustand sind Absorptions- und Emis- bind 
sionsspektrum nicht wesentlich voneinander verschieden. Die Wechselwirkung | Zusta 
mit den benachbarten Molekülen verändert das Spektrum nur wenig; es tritt nur der | 
eine geringe Rotverschiebung ein. Bei den bekannten anorganischen Leucht- | '*t- | 
stoffen hingegen handelt es sich um Übergänge, die im wesentlichen durch den | Wen! 
Aufbau des festen Körpers selbst erst hervorgebracht sind. In unserem Falle Erkl 
können wir also hier in gröbster Näherung so rechnen, als ob es sich um die freien noch 
Moleküle handelt. Nun rührt das beobachtete Fluoreszenzspektrum nur von dem auf « 
Übergang der beiden lockerst gebundenen Doppelbindungselektronen her. Bei sche: 
der Anregung durch «-Teilchen und Elektronen werden aber nicht nur wie bei Kop 


der Lichtanregung diese beiden lockerst gebundenen Elektronen angeregt, sondern : N 
alese 


*) Anmerkung bei der Korrektur: Nach neueren Messungen u unserem eer 
von der Größenordnung 10-8 bis 10-7 sec. 
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alle n-Doppelbindungselektronen und außerdem noch die übrigen Elektronen des 
Molekiils. Beim Phenanthren sind das 14 Doppelbindungselektronen und 52 an- 
dere Elektronen, wenn man von den zwei K-Elektronen jedes Kohlenstoffatoms 
absieht, weil diese auch bei Elektronenanregung nur mit geringerer Wahrschein- 
lichkeit angeregt werden. Betrachtet man zunächst nur die Doppelbindungs- 
elektronen und bei diesen auch nur die Anregung, soweit sie in den energetisch 
tiefsten unbesetzten Zustand des Energiebandes der Doppelbindungselektronen 
führt, so ist die Energie, die bei der Anregung durch x Teilchen- oder Elektronen- 
stöße auf diese Elektronen übertragen wird, natürlich erheblich größer als die 
mittlere Energie der optischen Absorptionsbande (einmal, weil auch noch andere 
Doppelbindungselektronen als die beiden lockerst gebundenen angeregt werden 
und sodann weil in demselben Energieband auch noch unbesetzte Zustände vor- 
handen sind, die energetisch merklich höher liegen als die tiefsten unbesetzten 
Zustände, in die die optische Absorption erfolgt). Man kann die Energie, die bei 
dem hier betrachteten Prozeß im Mittel auf die Elektronen übertragen wird, un- 
gefähr aus der Lage der durch die Rechungen von Hückel!!) bekannten Energie- 
niveaus abschätzen und kommt dabei zu einem Wert, der sicherlich doppelt so 
groß ist wie die mittlere Energie der langwelligen Absorptionsbande. Nun ist die 
mittlere Energie der Emissionsbande selber schon wesentlich kleiner als diejenige 
der Absorptionsbande, was einen weiteren Energieverlust ergibt. Wenn also die 
Anregung nur in die Zustände des untersten Energiebandes der Doppelbindungs- 
elektronen stattfinden würde, so könnte man höchstens mit einer Energie- 
ausbeute von 25—30% rechnen. Es werden außerdem noch energetisch höhere 
Zustände angeregt, die nicht dem untersten Energieband dieser Elektronen 
angehören. Es kommt hinzu, daß die Zahl der Doppelbindungselektronen über- 
haupt nur den fünften Teil der Gesamtzahl aller angeregten Elektronen ausmacht. 
Selbst wenn man nun annimmt, daß die Zustände höherer Energie oder auch die 
Elektronen, die von dem Molekül völlig losgelöst werden, ihre Energie, sobald sie 
dazu ausreicht, immer wieder in der Hauptsache zur Anregung tiefer Elektronen- 
zustände verwenden, so kann trotzdem auf das unterste Energieband der Doppel- 
bindungselektronen nur weniger als !/;, der gesamten Energie, die im Körper um- 
gesetzt wird, überführt werden. Man käme also zu einer Ausbeute von höchstens 
5%, wenn man wie oben berücksichtigt, daß nur etwa 25—30% durch Anregung 
in der fraglichen Fluoreszenzbande zur Emission kommt. Die tatsächlich gemessene 
Ausbeute ist aber größer. Offenbar gehen also alle Elektronen des Doppel- 
bindungs-Elektronenbandes, wenn sie angeregt werden, immer so in den Normal- 
zustand zurück, daß als einziger Übergang unter Lichtaussendung praktisch 
der Übergang stattfindet, der unserer schmalen Fluoreszenzbande zuzuordnen 
ist. Übergänge zwischen höher und tiefer gelegenen Zuständen kommen also, 
wenn überhaupt, nur sehr selten vor. Aber diese Annahme allein reicht zur 
Erklärung der Höhe der Energieausbeute nicht aus. Es scheint vielmehr auch 
noch von den anderen Elektronen ein erheblicher Teil der Anregungsenergie 
auf das Energieband der Doppelbindungselektronen übertragen zu werden; zwi- 
schen diesen und den anderen Elektronen im Molekül besteht also eine merkliche 
Koppelung. 

Man kann die Diskussion dieser Ausbeutefragen aber auch so führen, daß man 
diese Ausbeute direkt mit der gemessenen Ausbeute an anderen Leuchtstoffen 


11) Hückel, 1. c.!). 
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rrr Diese beträgt 25—30% der Gesamtenergie bei «-Teilchenerregung im 
ZnS. Hier liegen die Verhältnisse nun ganz anders. Bei diesen Stoffen ist die 
spektrale Lage von Absorptions- und Emissionsbanden im wesentlichen durch die 
Struktur des festen Körpers bedingt. Alle Elektronen, von einem wie tiefen Energie- 
zustand sie auch herrühren und in einen wie hohen Energiezustand sie überführt 
werden, werden mit großer Wahrscheinlichkeit schließlich immer zu angeregten 
Elektronen am unteren Rande des ersten unbesetzten Energieniveaus des ZnS- 
Kristalls führen. Es sind also alle Elektronen bei den anorganischen Leuchtstoffen 
mit diesem Bande gekoppelt und kommen daher für die Emission über den zuge- 
setzten Aktivator in Frage. Der Unterschied gegenüber den organischen Leucht- 
stoffen besteht nicht nur darin, daß das Leuchten bei diesen Zinksulfiden von 
Fremdatomen herrührt, sondern auch darin, daß bei den organischen Stoffen nur 
ein Teil der Elektronen, nämlich die Doppelbindungselektronen, mit der beob- 
achteten Emissionsbande unmittelbar verknüpft ist. Geht man von dieser Vor- 
stellung aus, so sollte man erwarten, daß die Ausbeute bei den organischen Leucht- 
stoffen etwa !/, (eben weil nur !/, der Elektronen mit dem fraglichen Übergang ver- 
knüpft ist) der Lichtausbeute bei ZnS beträgt. Daß die gemessene Ausbeute größer 
ist, ist wohl wieder so zu verstehen, daß die übrigen Elektronen der organischen 
Leuchtstoffe mit diesen Doppelbindungselektronen auch eine gewisse Koppelung 
aufweisen. 


> 
r 


In der vorausgegangenen Diskussion haben wir uns nur mit den Eigenschaften 
des freien Moleküls beschäftigt. Es weisen aber eine ganze Reihe von Erschei- 
nungen darauf hin, daß die Ausbeute an Licht bei diesen Substanzen ganz wesentlich 
durch den Zustand der Substanz als Ganzes bedingt ist, also auf der Wechsel- 
wirkung mit den Nachbarmolekülen beruht; denn die Ausbeute ist nicht nur vom 
Aggregatzustand, sondern, wie wir gesehen haben, auch vom Erstarrungszustand 
abhängig. Wir meinen daher, daß auch der große Unterschied zwischen den Aus- 
beuten von Naphthalin und Phenanthren durch die Struktur des Kristalls und 
nicht durch die des Moleküls bedingt ist; denn auch die Untersuchung anderer Sub- 
stanzen z. B. von Anthracen, zeigt, daß die Ausbeute nicht irgendwie eine eindeutige 
Funktion der Struktur des einzelnen Moleküls ist. Wir glauben also, daß die ge- 
messene maximale Ausbeute im wesentlichen auch noch dadurch bedingt ist, daß in 
den einzelnen Kristallen eine mehr oder minder starke, von Kristall zu Kristall 
verschiedene Auslöschung (Killer-Wirkung) vorhanden ist. Es wäre also in die- 
sem Sinne durchaus denkbar, daß z.B. die Fluoreszenz-Ausbeute des Naphthalins 
sich durch geeignete Behandlung des Kristalls noch wesentlich steigern ließe, und 
auch die von uns gefundene maximale Ausbeute des Phenanthrens braucht noch 
nicht die theoretisch maximal mögliche zu sein. Dabei wird die auslöschende 
Wirkung durch die benachbarten Eigenmoleküle bedingt. Daß eine solche Aus- 
löschung durch die Moleküle der eigenen Art besonders stark ist, geht aus den wohl- 
bekannten Versuchen von Wawilow!?2) hervor, der als erster bei der optisch ange- 
regten Fluoreszenz dieser Stoffe in Lösung mit wachsender Konzentration der art- 


12) Wawilow, Z. Physik 31, 750 (1925); Z. Physik 48, 249 (1927); Wawilow u. 
Lewschin, Z. Physik 35, 920 (1926); Wawilow u. Pheofilov, C.B. U.R.SS, 34, 
2220 (1942); Wawilow, J. Physics U. R. SS 7, 141 (1943); Wawilow, C.R. U. R. 8S 


42, 381 (1944); 45, 47 (1944), 
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eigenen Moleküle eine starke Auslöschung beobachtet hat. Wir glauben, daß die 
von uns beobachtete Abhängigkeit der Ausbeute von den äußeren Bedingungen 
durch eine ganz analoge Auslöschung hervorgerufen wird. 


Eine solche Auslöschung wird immer dann eintreten, wenn ein angeregtes Mole- 
kül mit einem unangeregten arteigenen Molekül in Wechselwirkung tritt, und der 
Abstand der Moleküle dabei so gering wird, daß bei dem Energieaustausch zwischen 
den beiden Molekülen die Potentialkurven der gegenseitigen Wechselwirkung der 
beiden Moleküle — des Zustandes, in dem eins der Moleküle angeregt ist und des- 
jenigen Zustandes, in dem beide unangeregt sind — sich genügend nahe kommen. 
Dann wird nach dem Frank-Condon-Prinzip, hier auf die Translationsenergie 
der beiden Moleküle angewendet, die Elektronenanregungsenergie mit erheblicher 
Wahrscheinlichkeit in kinetische Energie der Moleküle umgewandelt werden. 

Am größten wird eine solche Auslöschung sein, wenn eine möglichst geringe 
Symmetrie in der Anordnung der Moleküle gegeneinander vorhanden ist und die 
Abstände zwischen den Molekülen starken Schwankungen unterworfen sind. Daher 
ist es durchaus verständlich, daß nach unseren Beobachtungen mit der Verflüssi- 
gung der Substanz ein starkes Absinken der Fluoreszenz-Ausbeute eintritt. Daß 
die Auslöschung in dem kristallinen Zustand gering ist, besagt, daß die Moleküle 
in diesem Zustand nur sehr selten durch Wärmeschwingungen einander genügend 
nahe kommen, um eine solche Umwandlung der Elektronenenergie in kinetische 
Energie der Moleküle zu bewirken. Dies ist infolge der symmetrischen Lage der 
Moleküle im Kristall durchaus verständlich. Eine Auslöschung kann hier nur 
dann erfolgen, wenn eine genügend starke Abweichung von diesen symmetrischen 
Anordnungen auftritt. In der Flüssigkeit kommen solche starken Annäherungen 
der Moleküle während der Zeit der Lebensdauer des angeregten Zustandes natür- 
lich häufig vor;,daher ist bei diesen die Auslöschung relativ groß. 


Die in Abschnitt 7 beschriebenen Versuche zeigen, daß aber auch in der aus der 
flüssigen Schmelze erstarrten Masse noch eine zusätzliche Auslöschung vorhanden 
ist. Um eine solche Auslöschung zu erklären, genügt es wegen der Energiewanderung 
der Anregungsenergie anzunehmen, daß in der erstarrten Masse nur wenige Mole- 
küle vorhanden sind, die derart stark miteinander gekoppelt sind, daß bei ihnen 
ein Übergang der Anregungsenergie i in kinetische Energie bewirkt wird. Allerdings 
ist die Energiewanderung von Molekül zu Molekül im Falle der vorliegenden or- 
ganischen Moleküle sehr viel geringer als in anderen bekannten Fällen. Denn hier 
kommt dem Übergang von dem tiefsten angeregten Zustand zu dem unangeregten 
Normalzustand, wie aus den Absorptions- und Emissionsspektren zu entnehmen 
ist, nur eine geringe Übergangswahrscheinlichkeit zu. Da die angeregten Mole- 
küle sehr schnell ihre Kernschwingungsanregung an das Gitter abgeben, bleibt im 
Kristall schließlich nur der elektronenmäßig angeregte Zustand ohne Kernschwin- 
gung übrig, und nur der Übergang dieses Zustandes in den unangeregten Normal- 
zustand ist für die Energiewanderung maßgeblich. Diesem Übergang entsprie :ht 
aber wegen seiner geringen Übergangswahrscheinlichkeit nur eine sehr geringe 
Dipolwechselwirkung. Trotzdem wird ein soicher Energieaustausch stattfinden, 
weil auch noch Wechselwirkungsenergien höherer Ordnung auftreten. Daß diese 
Koppelungen nicht sehr groß sind, sieht man schon daran, daß die Banden dieser 
Fluoreszenzspektren, wenigstens bei tiefen Temperaturen, außerordentlich scharf 
und gegenüber den Molekülen in Lösungen oder sogar im Dampfzustand nicht ver- 
breitert erscheinen. Wir werden also nur mit geringen Koppelungsenergien, die 
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für die Ausbreitung der Anregungsenergie auf die benachbarten Moleküle maß- 
geblich sind, zu rechnen haben. Trotzdem. wird diese Koppelungsenergie noch groß 
genug sein, um innerhalb der Lebensdauer der Moleküle, die zwischen 10-8 sec und 
10-7 sec liegen dürfte, sehr viele Übergänge zu benachbarten Molekülen zu bewirken. 
Wir werden also auch hier mit einer Ausbreitung der Anregungsenergie über den 
Kristall zu rechnen haben, ähnlich wie man sie bei anderen Stoffen annimmt. Trifft 
die Anregungsenergie (Exciton) bei ihrer Ausbreitung durch den Kristall auf einen 
Molekülkomplex, der so beschaffen ist, daß bei ihm eine Umwandlung der Anregungs 
energie in kinetische Energie mit erheblicher Wahrscheinlichkeit möglich ist, so wird 
Auslöschung erfolgen. Man wird also die geringere Ausbeute der erstarrten Schmelze 
gegenüber dem sublimierten Zustand durch die Annahme zu erklären haben, daß 
sich der ideale Kristallzustand in der Schmelze noch nicht völlig ausgebildet hat. 
In ihr ist entweder noch eine große Zahl von Molekülen vorhanden, die eine geringe, 
aber gegenüber dem normalen Kristall vergrößerte Wahrscheinlichkeit für Aus- 
löschung aufweisen, oder in ihr existieren relativ wenige Molekülkomplexe mit 
einer großen Wahrscheinlichkeit für Auslöschung; trifft die Anregungsenergie bei 
der Diffusion durch den Kristall auf diese relativ wenigen Molekülkomplexe, 
so wird Auslöschung erfolgen. Rechnet man mit einer Koppelungsenergie, 
deren » von der Größenordnung 1012 sec! ist, so würde es schon genügen, 
wenn etwa jedes 10*te Molekül völlige Umwandlung der Anregungsenergie in 
kinetische Energie hervorruft, um die experimentell gefundene Auslöschung zu 
erklären. - 


Die starke Auslöschung im flüssigen Zustand beruht bei dieser Betrachtungs- 
weise darauf, daß hier immer genügend viele Molekülkomplexe vorhanden sind, 
die eine sofortige Energieübertragung ermöglichen. 


Die Tatsache, daß durch Pulverisieren der in der Schmelze erstarrten Substanz 
die ursprüngliche Lichtausbeute wieder hergestellt wird, erklärt sich bei dieser Auf- 
fassung Uaraus, daß diejenigen Stellen, die eine starke Auslöschung bedingen, durch 
das Pulverisieren zerstört werden. Es ist ja bekannt, daß durch Reiben und Er- 
schütterungen eingefrorene instabile Zustände zerstört werden. 


Die sehr starke auslöschende Wirkung arteigner Moleküle im Gegensatz zu 
artfremden Molekülen, z. B. in Lösungen, ist auf die besondere Resonanzwechsel- 
wirkung zwischen artgleichen Molekülen zurückzuführen, die eine größere An- 
näherung und Anziehung der Moleküle hervorruft, so daß eine Annäherung der 
fraglichen Potentialkurven ermöglicht wird. Es könnte hierbei auch noch der von 
Förster!®) diskutierte Punkt eine Rolle spielen, daß bei zwei nahe benachbarten 
Molekülen der angeregte Zustand in zwei Zustände aufspaltet, von denen der tiefere 
Zustand meistens eine geringere Emissionswahrscheinlichkeit und damit eine 
längere Lebensdauer hat. 


d) Ausbeute bei x-Teilchenanregung 


Schließlich hätten wir noch die Beobachtung zu diskutieren, daß die Licht- 
ausbeute bei Anregung durch x-Teilchen, bezogen auf gleiche absorbierte Energie, 
immer geringer ist als bei Anregung durch f- und y-Strahlen. Diese Erscheinung 


18) Förster, Ann. Physik 2, 55 (1948). 
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ist bei allen bisher untersuchten Substanzen mit alleiniger Ausnahme der Alkali- 
sulfide beobachtet worden. Nun ist der einzige wesentliche Unterschied zwischen 
Anregung durch schnelle Elektronen und Anregung durch «-Teilchen der, daß 
die Elektronenerregungsdichte bei «-Teilchenanregung sehr viel größer ist als 
bei Anregung durch schnelle Elektronen. Rechnet man sich die Wahrscheinlichkeit 
aus, mit welcher ein x-Teilchen bei seinem Weg durch den Kristall ein Elektron an- 
regt, so findet man, daß fast jedes Molekül, an dem das «-Teilchen vorbeifliegt, 
angeregt wird, während bei Elektronenanregung zwischen den einzelnen ange- 
regten Molekülen viele unangeregte liegen. Man kann dies überschlagsmäßig 
folgendermaßen berechnen: Ein a-Teilchen von 5 - 106eV hat eine Reichweite 
von etwa 5 mg Substanz/cm?. Auf 5 mg/em? kommen im Phenanthren etwa 
5- 104 Molekülschichten. Da das die Schicht durchsetzende «-Teilchen auf dieser 
Strecke etwa 10% primäre und sekundäre Anregungen erzeugt, wird praktisch zu- 
nächst jedes Molekül, das in der unmittelbaren Nähe der Bahn des a-Teilchens 
liegt, auch angeregt. Es wird also bei x-Teilchenerregung die Umgebung der 
erregten Moleküle vom normalen Zustand weit entfernt sein. Wir möchten an- 
nehmen, daß die Wechselwirkung dieser stark angeregten Moleküle miteinander 
wie Störstellen im Kristall wirken und eine zusätzliche Auslöschung im Kristall 
hervorrufen, die durch die geringe Ausbeute bei «-Teilchenerregung bewirkt wird. 


. 


für physikalische Chemie und Elektrochemie. 


(Bei der Redaktion eingegangen am 10. August 1948.) 
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Untersuchungen über die Bewegung des Brennflecks 
auf der Kathode eines een 


N 


(Mit 24 Abbildungen) 


Inhaltsiibersicht 


1. Mit Hilfe von Zeitlupenaufnahmen eines Hg-Niederdruckbogens wird bei 
höheren Stromstärken eine Aufteilung der kathodischen Entladung in einzelne 
Brennflecke gefunden. Dabei beträgt die statistisch schwankende Stromstärke 
eines einzelnen Flecks im Mittel etwa 6,7 A. (Kap. IIa.) 

2. Die Unruhebewegung eines einzelnen Flecks entspricht rein formal der 
Brownschen Bewegung eines kleinen Teilchens in einem ebenen Gas. Das für die 
Molekularbewegung gültige Gesetz, daß das mittlere Verschiebungsquadrat des 
Teilchens der Beobachtungszeit proportional sein muß, läßt sich auch auf die zu- 
fällige Brennfleckwanderung anwenden. Da das Verschiebungsgesetz auch noch 
bei den kleinen Beobachtungszeitintervallen von 4 - 10-4 s die Verhältnisse richtig 
wiedergibt, muß man annehmen, daß der zeitliche Abstand der den Brennfleck 
ungeordnet umlenkenden Vorgänge um mindestens eine Zehnerpotenz geringer ist 
als diese kürzeste Beobachtungszeit. (Kap. IIb.) 

3. Auf den Brennfleck wird ferner bei geeigneter Form der Entladungsstrecke 
eine Richtkraft ausgeübt, die ihn in eine der Verkürzung der Bogensäule ent- 
sprechende Richtung treibt. (Kap. IIc.) 

4. Es wird vermutet, daß diese Richtkraft mit dem Steenbeckschen Mini- 
mumprinzip ursächlich zusammenhängt und daß sich die Häufigkeitsverteilung 
des Brennflecks über den Kathodenquerschnitt durch einen Boltzmann -Ansatz 
beschreiben läßt, dessen Exponent durch das Verhältnis einer mit der Brennspan- 
nung monoton steigenden Funktion zur thermischen Unordnungsenergie ge- 
geben ist. (Kap. IIla.) 

5. Statistische Messungen der Häufigkeitsverteilung auf einer ebenen Kathode 
im Zusammenhang mit der Brennspannungsverteilung ergeben, daß diese Annahme 
(Abs. 4) nicht zu Recht besteht. (Kap. IIIb e). 

6. Der an einer linearen (in einer Richtung stark ausgedehnten) Kathode ge- 
messene Zusammenhang zwischen der örtlichen Häufigkeitsverteilung des Brenn- 
flecks und den zugehörigen Brennspannungswerten läßt sich sehr genau durch einen 
Boltzmann-Ansatz wiedergeben. (Kap. IIId e.) 

7. Die zwischen den Exponenten des Boltzmann-Ansatzes, der Diffusions- 
konstanten des Brennflecks und seiner Driftgeschwindigkeit nach der Theorie der 


1) Binustetion TH Berlin 1943. Zur Veröffentlichung in Bd. XXII der Wiss. Ver- 
öffentl. Siemens eingereicht am 14. 12. 1943, dessen Herausgabe infolge der Kriegsereig- 
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statistischen Schwankungserscheinungen bestehenden Zusammenhänge erweisen 
sich für den Fall der linearen Kathode als befriedigend erfüllt. (Kap. IV.) 

8. Ohne Zuhilfenahme weiterer Annahmen ist eine Deutung der bisherigen 
Meßergebnisse durch das Minimumprinzip nicht möglich. (Kap. V.) 

9. Durch Untersuchung des Einflusses von transversalen Magnetfeldern auf 
die Driftgeschwindigkeit des Brennflecks wird auf den Feldeinfluß des Entladungs- 
stromes geschlossen. Berücksichtigt man einen von J. Stark gefundenen Effekt, 
bei dem der Brennfleck vom transversalen Magnetfeld im elektrodynamisch 
falschen Sinne abgelenkt wird, so lassen sich die Meßergebnisse allein durch die 
Wirkung des magnetischen Eigenfeldes der Entladung deuten. (Kap. VI.) 


I. Einleitung 


Der kathodische Brennfleck des Quecksilberdampfbogens ist bereits Gegenstand 
zahlreicher experimenteller und theoretischer Untersuchungen gewesen, die haupt- 
sächlich zur Klärung der Vorgänge im Brennfleck selbst beitragen sollten. Vor- 
liegende Arbeit beschäftigt sich mit einer anderen den Brennfieck betreffenden 
Frage, bei der auch die Stromspannungsverhältnisse und die Abmessungen der Ent- 
ladungsstrecke eine Rolle spielen. 

Wahrscheinlich unter dem Einfluß der heftigen Verdampfungsvorgänge im 
Brennfleck, dessen Stromdichte bei einigen tausend A/cm? liegt, führt der Fleck 
auf der die Kathode bildenden Quecksilberoberfläche eine sehr lebhafte Bewegung 
aus. Diese Erscheinung, die man z. B. an Quecksilberdampf-Glasgleichrichtern 
beobachten kann, sofern keine besonderen Maßnahmen getroffen sind sie zu ver- 
meiden, ist gewöhnlich vollkommen regellos. Durch Anordnung von Molybdän- 
oder Wolframblechen in der Quecksilberoberfläche ist eine Fixierung des Brenn- 
flecks an das Halteblech möglich, die hauptsächlich eine Benetzung des Bleches 
mit Quecksilber zur Voraussetzung hat. Verzichtet man auf diese im Gleich- 
richterhau wegen der besseren Konstanz der Ionisierungs- und Dampfströmungs- 
verhältnisse häufig verwendete Maßnahme, so verursacht der umherlaufende 
Fleck eine kräftige Wellenbewegung des Kathodenquecksilbers. 

Die auffälligen Bewegungserscheinungen des nicht fixierten Brennflecks sind 
bisher noch nicht genauer untersucht worden. Sie‘stellen den Gegenstand vor- 
liegender Arbeit dar. Da sie nicht ausschließlich statistischer Natur sind, sondern 
bei geeigneter Form des Entladungsgefäßes durch eine Überlagerung von unge- 
ordneten und gerichteten Kräfte zustande kommen, ergibt sich nach Beschreibung 
der Vorgänge in Kapitel II die Aufgabe, eine Trennung der ungeordneten und gerich- 
teten Größen vorzunehmen. Dies geschieht in den Kapiteln III und IV. In den 
Kapiteln V und VI wird eine Deutung der Beobachtungsergebnisse versucht. 


I Zeitlupenaufnahme des Brennflecks 
Den Ausgangspunkt für die Untersuchungen lieferten einige vom Verfasser 
bei Aufnahmen des Brennflecks mit der Filmkamera beobachtete besondere Vor- 
gänge. Bei den Aufnahmen fand die Zeitlupe von Zeiss-Ikon Verwendung, mit der 
eine Bildfrequenz von -3000/s erreicht werden kann. Abb. 1 gibt das benutzte 
Versuchsgefäß wieder, das eine Aufnahme des Brennflecks senkrecht zur Queck- 


silberoberfläche gestattet. Die Auswertung des Filmstreifens ergab folgendes: wr 
im 
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a) Aufteilung des Brennflecks in Einzelflecke 


Der vom 220 V-Gleichstromnetz über ohmsche Widerstände mit der Haupt- 
anode betriebene Bogen zeigt, daß sich der Kathodenbrennfleck in mehrere Einzel- 
flecke aufteilt. Schon A. Güntherschulze ?) hat gefunden, daß dieser Effekt 
statistischer Natur ist und bei Stromstärken über 10 A deutlich zu werden be- 
ginnt. Wie aus den Filmaufnahmen hervorgeht, macht sich die Erscheinung bereits 
bei Stromstärken von wenigen Ampere w(4) 


bemerkbar, wenn auch hier die Wahr- 
scheinlichkeit einer Aufteilung gering ist. | 
Durch Auszählen der Fleckzahl x je r GH | 
Filmbild erhält man bei Division des Ent- GG 

ladungsstromes J durch x eine mittlere FA 1-157 | 


Stromstärke je Einzelfleck i = I/x. Dieser 
Wert zeigt auch bei konstant gehaltenem 
ZT innerhalb der verschiedenen Aufnahmen 
erhebliche Streuungen. Als Beispiel gibt 
Abb. 2 die durch Auswerten von je 
180 Filmbildern ermittelte Häufigkeits- 
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Abb. 1. Gefäß für Filmaufnahmen Abb. 2. Statistische Häufigkeit der An- 
ur des Brennflecks : zahl der Brennflecke je Einheit des Ent- 
ladungsstroms 


x 
verteilung von x/J für einen Entladungsstrom von J = 16,3 und 75,7 A wieder. Auf- 
getragen ist die Wahrscheinlichkeitsfunktion w(x#/I) gegen x/I, wobei derart normiert 
wurde, daß das Wahrscheinlichkeitsintegral über alle x/J den Wert 1 erhält. Man 
erkennt aus dem Auftreten der großen x/I- bzw. der kleinen i-Werte, daß sich der 
Fleck auch bei geringen Entladungsströmen gelegentlich aufteilt. Die Häufig- 
keitsmaxima liegen für die beiden Entladungsströme an gleicher Stelle x/I bzw. 
i= I/x. Letzteres stimmt ungefähr mit den arithmetischen Mittelwerten von 4 
überein, für die nachstehende Werte der Zahlentafel 1 gefunden wurden. Der in 
der Zusammenstellung angegebene Quecksilberdampfdruck pug wurde durch 
Messung der Temperatur an der Kondensationsgrenze ermittelt. 


2) A. Güntherschulze, Die Vorgänge an der Kathode des Quecksilbervakuum- 
lichtbogens. Z. Physik 11, 74 (1922). 
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Zahlentafel 1 


Pus I i 
Tk rT A A 
33-108 16,3 6,20 aus 180 Aufnahmen berechnet 
83-10 75,7 6,78 „ 180 ” = 
1,8- 10-2 17,9 7,16 » 90 
18-1072 75,7 6,81 we 


Unter Beriicksichtigung der der Anzahl der Beobachtungen entsprechenden 
Gewichte berechnet sicn für das anscheinend druck- und stromunabhingige i 
ein Mittelwert aus allen Messungen von 6,7 A. Dieser Wert ist unter Umständen 
mit einem systematischen Fehler behaftet, der durch eine gewisse Unsicherheit 
bei der Auszählung der Brennflecke zustande kommt. Durch die stark bewegte 
Quecksilberoberfläche treten leicht Reflexe auf, die für Brennflecke gehalten 
werden können. Um die Täuschungsmöglichkeit gering zu halten, wurden bei den A 
Aufnahmen Öffnungsverhältnis und Belichtungszeit so klein gewählt, daß die 
Brennflecke auf dem Filmstreifen gerade noch erkannt werden konnten. Zur 
Verwendung kam Agfa-ISS-Umkehrfilm. Die Bildfrequenz betrug rund 3000/s. 
Dabei wurde die vor dem Filmbild befindliche Schlitzblende auf eine wirksame 
Belichtungszeit von 1/30000 s eingestellt und das Öffnungsverhältnis auf 1:17,6 
herabgesetzt. Die lichtschwächeren Reflexe dürften damit zum größten Teil unter- 
drückt worden sein. Bei der Auswertung wurden im übrigen jeweils zwei Flecke 
als einer gewertet, wenn ihr Mittelpunktsabstand bei einem scheinbaren Fleck- 
durchmesser von etwa 0,4 mm kleiner als 0,5 mm war. Hierdurch entstehende 
Zählfehler wirken dem durch Reflexe verursachten entgegen. 

Erwähnt sei noch die Tatsache, daß bei gleichzeitigem Brennen mehrerer 
Flecke ein einzelner Fleck plötzlich verschwinden kann, ohne sich einem anderen 
bis zum Zusammenlaufen nähern zu müssen. Das ist verständlich, da es sich hier 
um parallel geschaltete Entladungsstrecken handelt, bei denen einzelne Teilströme 
zugunsten anderer bis zum Erreichen einer Stabilitätsgrenze abnehmen können und 
verlöschen. Dagegen ist das Neuentstehen eines Fleckes nur durch einen Teilungs- 
vorgang möglich, weil auf der freien Kathode die für eine spontane Neuzeugung 
erforderliche hohe Trägerdichte fehlt. Damit erklärt sich die aus den Aufnahmen 
ersichtliche Tendenz der Brennflecke, in Gruppen zusammenzubleiben, ohne Zu- 
hilfenahme elektrodynamischer Vorgänge. 


b) Verhalten eines einzelnen Brennflecks 

Das anscheinend völlig ungeordnete Umherlaufen des Brennflecks auf der 4 
Kathode legt den Gedanken nahe, zu prüfen, ob sich der Brennfleck wie ein Teil- 
chen verhält, das eine Brownsche Bewegung ausführt. Sollte das der Fall sein, 
so müssen die in Zeitintervallen At auftretenden Verschiebungen Ax der Gesetz- 


mäßigkeit Az? ~ At genügen ®). Voraussetzung für die Gültigkeit dieser Beziehung \4 
ist, daß At so groß gewählt wird, daß die zur Zeit t vorhandene Geschwindigkeit a 
zur Zeit t + At nach Größe und Richtung völlig durch eine ausreichende Anzahl | 
neuer Umlenkungen abgeklungen ist. Nach oben hin ist At dadurch begrenzt, daß . a 
die zugehörigen Ax-Werte zur Vermeidung von Wandberiihrungen nicht die 4 


Größenordnung des Kathodendurchmessers erreichen dürfen. 


3) Siehe z. B. G. Joos, Lehrbuch der theoretischen Physik (1932) S. 496. 
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Eine Nachprüfung des Gesetzes wurde durch Auswertung eines Filmstreifens 
durchgeführt, bei dem der Entladungsstrom 5,5 A betrug. Bei diesem Strom ist 
nach dem Vorherigen praktisch mit nur einem Brennfleck zu rechnen. Der Bogen 
brannte diesmal, durch Widerstände stabilisiert, gleichzeitig nach beiden Hilfs- 
anoden hin. (Vgl. Abb. 1.) Die Auswertung erfolgte in der Weise, daß für ein ge- 
gebenes At, z. B. für das durch den zeitlichen Abstand der nacheinander folgenden 
Filmbilder bedingte kleinste 5.105 (die höchste Bildfrequenz war bei 
dieser Aufnahme nicht ganz erreicht), die von Bild zu Bild auftretenden Fleck- 
verschiebungen Ax gemessen und von diesen Werten die vorzeichenunabhängigen 
Quadrate gebildet wurden. Der arithmetische Mittelwert der Quadrate ist die zu 
bestimmende Größe Az? Die zu der willkürlich festgelegten x-Richtung senk- 

rechten Verschiebungen Ay wurden in 

gleicher Weise ausgewertet. Da keine Vor- 
zugsrichtungen vorhanden sind, muß ihr 

mittleres Quadrat gleich dem der z- 
Richtung sein und darf für die Aus- 
| wertung mit verwendet werden. Der Kürze 
| halber sei der sich durch Mittelung der 
| . Werte für beide Richtungen ergebende Be- 


Os 4 trag weiterhin mit Ax? bezeichnet. Als 
größtes Zeitintervall wurde At= 0,1s 
” gewählt, ein Wert, bei dem die Wand- 
| | störungen noch zu vernachlässigen sind. 

0 4 é 
” yar 03 5 In Abb. 3 ist das Ergebnis der Auswertung 
wiedergegeben. Aus Griinden der Dar- 
Abb. 3. Zeitliche Abhängigkeit des 
mittleren Verschiebungsquadrates stellungsgenauigkeit ist nicht Az? in 


Abhängigke#t von At, sondern VAx: als 
Funktion von VAt aufgetragen. Die aus etwa je 100 Verschiebungen berechneten 
MeBpunkte.liegen innerhalb des zu erwartenden Genauigkeitsbereiches auf emer 
durch den Nullpunkt gehenden Geraden. 

Auf Grund dieses Befundes stellen wir also fest, daß sich ein einzelner Brenn- 
fleck ideal ungeordnet bewegt. Eine Bahngeschwindigkeit des Flecks läßt sich 
aus den Filmaufnahmen nicht ermitteln, da ja die Giiltigkeit des Verschiebungs- 
gesetzes bis zu den Zeitintervallen von 4 - 10-4 s herab besagt, daß innerhalb dieser 
Zeit noch viele Umlenkungen erfolgen müssen. Sollte sich überhaupt eine mittlere 
Bahngeschwindigkeit angeben rn so müßte sie nach vorliegenden Messungen 


noch erheblich größer sein als V. Ar: min/Atmin = 2 10? em s!. A. Günther- 
schulze?) rechnet mit Werten von 10% em s=!. 

Abschließend kann man also unter Verzicht auf große Feinheiten die Bewegung 
eines einzelnen Brennflecks rein formal so betrachten, als handelte es sich um ein 
kleines Partikel, das in einem ebenen Gas mit gegebener Temperatur ständigen 


thermischen Stößen ausgesetzt ist. 
e) Einfluß der Lage der Anode 


Bei den Aufnahmen mit höheren Stromstärken wie unter a) wurde beobachtet, 
daß die Brennflecke nach einer Brenndauer von wenigen Sekunden stets die der 
Hauptanode (s. Abb. 1) zugekehrte Seite der Kathodenfläche erreichten, wo sie 
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dann in unmittelbarer Nähe der Glaswand hin- und hertanzend weiterbrannten. 
Dagegen ließen sich keinerlei Vorzugsrichtungen erkennen, wenn die dicht an der 
Kathode liegenden Hilfsanoden verwendet wurden. 


III. Statistische Untersuchungen über die örtliche Häufigkeitsverteilung iR 
des Brennflecks auf dem Kathodenquersehnitt 
a) Theoretische Vorbemerkungen 


Die merkwürdige Erscheinung, daß sich der Brennfleck im Sinne einer Ver- 
kürzung der Bogensäule zu bewegen sucht, legt den Gedanken nahe, nach einem 
Zusammenhang mit der von M. Steen beck*) bei Untersuchungen am Luftlicht- 


"bogen im schwerefreien Raum gefundenen Bogenzugspannungen zu suchen. 


Unter der Annahme, daß die auf den Brennfleck wirkende Richtkraft durch 
den Gradienten eines Potentials @ (x, y) beschrieben werden kann, muß für die 
über genügend lange Zeiten zu messende Aufenthaltswahrscheinlichkeit des Brenn- 
flecks an der Stelle (x, y) der Kathodenfläche der Boltzmann-Satz gelten: 


_ G(x, ¥) - (0,0) 
kT 


= . 
Die Größe n, mit der wir. die Aufenthaltswahrscheinlichkeit ausdrücken wollen, 
soll hierbei folgende Bedeutung haben: In gleichmäßigen Zeitabständen sei durch 
sehr viele, kurzzeitige Momentaufnahmen der Ort des Brennflecks ermittelt und 
durch je einen Punkt im Koordinatensystem (x, y) eingetragen. Die Größe n ist 
dann die Anzahl der Punkte auf der Flächeneinheit. , ist die Punktdichte an einem 
willkürlich gewählten Bezugspunkt. Da wir über die Unruhebewegung des Brenn- 
flecks nur aussagen können, daß sie sich im beobachteten Bereich mit der Theorie 
der statistischen Schwankungserscheinungen beschreiben läßt, erscheint es nicht 
ratsam, eine fiktive thermische Energie mit kT' zu bezeichnen, was nur für mole- 
kulare Stöße auf eine definierte Masse einen Sinn hätte. Wählen wir an Stelle von 
kT das Zeichen E, so wird, indem wir gleichzeitig das Potential für den Bezugs- 
punkt (0,0) gleich Null setzen, 
_ 
n=me . (1) 


Nehmen wir an, daß die Unruheenergie über den Kathodenquerschnitt kon- 
stant bleibt, so ist fiir vorgegebenen Entladungsstrom und Druck die unbekannte 


Potentialfunktion bis auf den Faktor E durch das meßbare Verhältnis = zu be- 


0 
stimmen. Damit ergibt sich also zunächst die Aufgabe, zu untersuchen, von welchen 
Parametern g abhängt. 

Einen Hinweis für die Lösung des Problems scheint das Steenbecksche 
Minimumprinzip zu liefern 5), dem folgender Gedanke zugrunde liegt: Enthält ein 
gleichstromdurchflossener Stromkreis eine Entladungsstrecke, so bildet sich die 
Entladung im stationären Fall derart aus, daß ihre Brennspannung den kleinst- 
möglichen Betrag annimmt. Wenn auch dieser Regel nur bei Erfüllung besonderer 


4) M. Steenbeck, Untersuchungen am Luftlichtbogen im schwerefreien Raum, 
Z. techn. Physik 18, 593 (1937). 

5) M. Steenbeck, Eine Prüfung des Minimumprinzips für thermische Bogensäulen 
an Hand neuer Meßergebnisse. Wis. Veröff. Siemens-Werken, XIX, 59 (1940). 
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Voraussetzungen exakte Gültigkeit zukommt®), so hat ihre Anwendung doch in 
sehr verschiedenen Fällen bei der Berechnung stationärer Entladungsformen gute 
Übereinstimmung mit der Wirklichkeit ergeben. Als Arbeitshypothese wollen wir 
daher zunächst annehmen, daß qm durch irgendeine mit der Brennspannung U 
monoton steigende Funktion gegeben ist. Zur Ermittlung dieses Zusammenhanges 
ist die Abhängigkeit n (U) zu messen. 

Um die experimentell unbequeme, gleichzeitige Messung von Ort und Brenn- 
spannung des Kathodenflecks zu umgehen, wurde ein statistisches Verfahren ange- 
wandt, das, sofern die Brennspannungserhöhung die alleinige Ursache für die durch 
die Boltzmann -Verteilung entstehende Fleck- 
dichteverminderung in anodenferneren Ka- 
thodengebieten ist, eine zeitliche Trennung 
der Messung von Brennspannung und Fleck- 
verteilung gestattet. Es beruht auf folgender 
Überlegung: Denken wir uns auf der Kathode 
Linien konstanter Fleckdichte » = const ein- 
getragen. Da die Brennspannung U für ge- 
gebene Geometrie, konstanten Strom und Druck 


n=const. 


U = const. 
Abb. 4. Zur Messung der örtlichen 
Dichteverteilung 


lediglich eine Funktion von 7 sein soll, gehört 
zu jeder dieser Linien auch ein bestimmter 
Wert U. In Abb. 4 ist eine derartige Linie 


n = const, U = const aufgezeichnet. Die Anode 
befindet sich links von ihr, so daß U im Sinne der Pfeilrichtung zunimmt, während 
n in dieser Richtung abnimmt. Bezeichnen wir mit w (n)-dn die Wahrschein- 
lichkeit dafür, daß bei einer Messung der Brennfleck in den Dichtebereich 
n..n-+ dn fällt, und mit g (U) dU die Wahrscheinlichkeit dafür, daß die Brenn- 
spannung im Bereich U ..U + dU liegt, dann können wer die Wahrscheinlich- 
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Abb. 5. Zuordnung der statistischen Häufigkeiten von n und U 


keit dafür, daß ein Brennfleck rechts von der Linie auftritt, auf zweierlei Weise 


n oo 
ausdrücken. Sie beträgt nämlich einmal f w(n)dn und andererseits f g(U)dU. 
6 


fw (n) dn = fav) du. (2) 
Ü 


0 


Gleichsetzung liefert 


Gelingt es also durch getrennte Messungen die beiden Funktionen w(n) und g(U) 
zu bestimmen, so lassen sich durch graphische Integration nach Gl. (2), wie in 
6) R. Rompe u. W. Weizel, Über die Bedeutung des Steenbeckschen Minimum- 
prinzips. Z. Physik 120, 31 (1942). 
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Abb. 5 dargestellt, die Wertepaare n und U einander zuordnen. Damit erhält man 
punktweise die gesuchte Beziehung n (U). 

Die Ermittlung der Funktion w (n) geschieht grundsätzlich nach folgendem 
Verfahren: Durch eine große Anzahl N von Momentaufnahmen sei die Häufigkeits- 
verteilung ‘des Brennflecks auf dem Kathodenquerschnitt festgelegt. Teilt man 
nun die Kathodenfläche in Zonen ein, derart, daß in diesen die Fleckdichte zwischen 
den Grenzen n und n + An liegt, so ist w(n) An als die Wahrscheinlichkeit für 
das Auftreten des Brennflecks in einer Zone gegeben durch das Verhältnis der 
Anzahl der Punkte in diesem Bereich zur Gesamtzahl N. Bezeichen wir den 
Flächeninhalt des Zonenbereiches (n .. + n) mit f (n, An), so wird 


w(n) - An = In) 


Über die praktische Durchführung des Verfahrens wird bei der Besprechung der 
Meßergebnisse berichtet. 

Die Messung der statistischen Verteilung der Brennspannungswerte erfolgt 
zweckmäßig in einer Schaltung, die gleich das für die Auswertung nach Gl. (2) 


co 
benötigte Integral [g(U)dU liefert. Das ist auf folgendem Wege möglich: 
U 


Eine periodisch wirkende Schalteinrichtung gibt in gleichmäßigen Zeitabständen 
für sehr kurze Meßzeiten eine Registriervorrichtung frei, die nur dann anspricht, 
wenn die Brennspannung U einen bestimmten, beliebig wählbaren Wert U, über- 
schreitet. Nennen wir die Gesamtzahl der von der Schaltung untersuchten Brenn- 
spannungswerte Z und die hierbei erfolgenden Registrierungen z (U,,), so erhalten 
wir aus dem Verhältnis beider das zu bestimmende Integral 
1, 


Voraussetzung für die Brauchbarkeit des Verfahrens ist, daß die Meßzeit für 
einen Brennspannungswert so kurz gewählt wird, daß der Brennfleck während dieser 
Zeit praktisch stillsteht, d. h., daß U als konstant anzusehen ist. Dagegen ist der 
zeitliche Abstand zwischen den einzelnen Messungen im allgemeinen so groß zu 
wählen, daß keine statistische Nachwirkung eintritt. Bei einer gegenüber der Nach- 
wirkungszeit (rs Laufdauer des Brennflecks über Kathodenquerschnitt) großen 
Gesamtbeobachtungszeit, wie sie in dieser Arbeit sowohl bei Bestimmungen der 
Funktion w (n) als auch für die Brennspannungsverteilung vorliegt, tritt jedoch 
durch Nachwirkung keine Fälschung der Ergebnisse ein. 


b) Versuchsanordnung I (ebene Anordnung) 


Zur eindeutigen Messung der örtlichen Häufigkeitsverteilung des Brennflecks 
muß man mit geringen Bogenströmen arbeiten, um die Ausbildung von mehreren 
Brennflecken zu verhindern. Einfache seitliche Anbringung der Anode bringt dabei 
keine für unsere Zwecke ausreichende Fleckdichteveränderung über den Kathoden- 
querschnitt. Bewährt hat sich dagegen folgende, in Abb. 6 dargestellte Anordnung, 
bei der ein wirksames Ansteigen der Brennspannung und entsprechendes Abnehmen 
der Fleckhäufigkeit dadurch erfolgt, daß der Brennfleck unter eine die Kathode 
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in 1cm Höhe über dem Quecksilberspiegel teilweise bedeckende Glasplatte laufen 
kann. Abb.7 gibt schematisch das verwendete eiserne Versuchsgefäß, das einen 
raschen Wechsel der Kathodenanordnung gestattet, und die optische Anordnung 
wieder. Als Anode dient ein die Beobachtung von oben durch das Schauglas nicht 
behindernder Eisenring, der so weit von der Kathode entfernt ist, daß Bedeckungs- 
schwankungen der anodischen Entladung keinen Einfluß auf die Vorgänge an der 
Kathode haben können. Gezündet wird der Bogen durch die bei technischen Gleich- 
richtern gebräuchliche Quecksilberspritzeinrichtung, mit deren Hilfe gegen eine 
in den Kathodenraum hineinragende Molybdändrahtzündelektrode ein kurzzeiti- 
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Abb.6. Kathodenaufbau der Abb. 7. Versuchsaufbau I 


Versuchsanordnung I 


ger Abreißbogen gezündet wird. Das Gefäß besitzt zur besseren Festlegung des 
Dampfdruckes einen einfachen wasserdurchflossenen (20° C) Kühlmantel und eine 
nicht dargestellte Umwälzpumpe mit Rückkühler. Ein Warmluftstrom gegen das 
vakuumdicht eingeschmolzene Schauglas verhindert ein Beschlagen mit Queck- 
silber. Das Gefäß liegt ständig an der Diffusionspumpe. Um den Bogenstrom gegen- 
über den Brennspannungsschwankungen konstant zu halten, wird mit einer hohen 
treibenden Spannung (440 V) gearbeitet. 

Die Fleckdichteverteilung wird mit der in Abb. 7 dargestellten optischen An- 
ordnung gemessen. Über einen Spiegel wird der dem Brennfleck zur Verfügung 
stehende kreisférmige Kathodenausschnitt in etwa natürlicher Größe auf der 
Mattscheibe einer großen Porträtkamera abgebildet. Eine von einem Zeitrelais 
gesteuerte, elektromagnetische Verschlußvorrichtung gibt das Objektiv alle 3s 
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für die Dauer von 1/25 s frei. Zur weiteren Herabsetzung der Belichtungsdauer 
auf 1/1500 s rotiert vor dem Objektiv eine von einem Synchronmotor mit 
1500 U/min angetriebene Schlitzscheibe. Damit erscheint das Bild des Brennflecks 
in Abständen von 3s auf der Mattscheibe praktisch punktförmig. Eine photo- 
graphische Registrierung hat sich nicht bewährt, weil einerseits das helle Licht 
der unter der Kathodenglasplatte stark eingeengten Bogensäule die Aufnahmen 
verschleiert und andererseits Reflexe des Brennflecks innerhalb der für die Aus- 
wertung erforderlichen großen Punktzahl irrtümlich mitgezählt werden könnten. 
Daher wurde die Mattscheibe mit durchscheinendem Koordinatenpapier bedeckt 
und der dort aufblitzende Leuchtfleck von Hand eingetragen. Nach einiger Übung 
ist es infolge der Trägheit des Auges ohne weiteres möglich, den Ort des Brenn- 
flecks auf + 1mm genau anzugeben und etwaige Reflexe vom Brennfleck zu 
unterscheiden. Die Bedeutung der in Abb. 7 gestrichelt eingezeichneten Beleuch- 
tungsvorrichtung und der 
in der Kamera eingebauten . 
unten erläutert werden. % eee 
Mit der in Abb. 8 “iter Impuls - 
schematisch dargestellten ger30 < 
Schaltung wird die Ver- | 
teilungsfunktion der Brenn- Zählwerk 
spannungswerte gemessen. 
Über die meßbar ver- 
änderliche Gegenspannung Abb. 8. Prinzip der Versuchsschaltung I 


U, liegt die Brennspannung 
nach einstufiger Verstärkung am Gitter 1 eines Doppelgitterstromtores. Gitter 2 wird 
in gleichmäßigen Zeitabständen mit kurzzeitigen Spannungsimpulsen beschickt, 
deren Gesamtzahl vom Zählwerk 2 angezeigt wird. Eine Zündung des Stromtores 
und damit ein Ansprechen des Zählwerkes 1 kann nur im Augenblick eines 
Spannungsimpulses auf Gitter 2 erfolgen, und zwar auch nur dann, wenn die Brenn- 
spannung einen gewissen, durch die Gegenspannung gegebenen Wert überschreitet. 
Da es bei den Messungen lediglich auf die durch das Umherlaufen des Brennflecks 
bedingten Brennspannungsänderungen, nicht aber auf einen langsamen, durch 
Anodenfallschwankungen entstehenden Spannungsgang ankommt, wurde für die, 
Messung die Potentialdifferenz zwischen Anode und einem Punkt verwendet, der 
sich auf den zeitlichen Mittelwert des Anodenpotentials einstellt. Hierfür dient ein 
RC-Glied im Verstärkereingang, dessen Zeitkonstante so groß sein muß, daß die 
durch die Fleckbewegung entstehenden Spannungsschwankungen ausgeglichen 
werden. Eine Zeitkonstante von 24 s war hierfür ausreichend. 

Abb. 9 zeigt die Einzelheiten der Schaltung. Die Gegenspannung wird mit der 
Gittervorspannung der Verstärkertriode eingestellt. Ihr Arbeitswiderstand liegt 
in der Kathodenleitung und ist zur Ausschaltung von höherfrequenten Störspan- 
nungen mit einem kleinen Kondensator überbrückt. Das ce I 
arbeitet mit Anodenwechselspannung auf den Gesprächszähler 1 (Zählangabe = 
Die Zündimpulse auf das zweite Gitter werden der Kondensatorentladung des 
Stromtores JJ entnommen, das von einem mit einer rotierenden Schalteinrichtung 
in Reihe liegenden Steuersatz alle 0,5 s im Spannungsmaximum der Anodenwechsel- 
spannung des Stromtores J gezündet wird. Nicht in die Schaltung eingezeichnet 
ist das Zählwerk 2, das, mit dem rotierenden Schalter gekuppelt, die Anzahl der 
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Spannungsimpulse auf Gitter 2 des Stromtores J registriert (Zählangabe = Z). ausg' 


Die verwendeten Gleichspannungen sind durch einen Glimmlichtstabilisator und einge 
die Wechselspannungen für die Röhrenheizung mit magnetischen Spannungs- unte 
reglern konstant gehalten. gren: 

Zur Inbetriebnahme der Anordnung wird der Verstärkereingang zunächst kurz- Zone 
geschlossen und die negative Gittervorspannung so eingestellt, daß das Zähl- Zone 


Nergerät S20V 


Ahb. 9. Versuchsschaltung I 


werk 1 gerade nicht anspricht. Das ist für die Dauer einer Meßreihe mit einer Ge- 
nauigkeit von + 0,01 V möglich. Bezeichnen wir diese Spannung mit U,, und die 
Differenz zwischen Momentanwert der Brennspannung und dem zeitlichen Mittel- 
wert U— U mit U’, so entspricht einer Vergrößerung der negativen Gittervor- 
spannung auf den Betrag U, ein Wert 


Beschränken wir uns darauf, an Stelle von n (U) die physikalisch gleichwertige Ab- Pur 
hängigkeit n (U’) zu ermitteln, so haben wir also statt Gl. (4) zu schreiben 


z(U U 0) 
= | g(U) av’. 
U’ 


c) Ergebnisse der Messungen an Versuchsanordnung I 


er Der Gang der Auswertung der Beobachtungsergebnisse wird am einfachsten 
= an Hand eines Beispieles erläutert. Für einen Entladungsstrom von 5,0 A wurde 
die in Abb. 10 wiedergegebene Häufigkeitsverteilung über den Kathodenquer- 
schnitt gefunden. Die Kathodenfläche ist in quadratische Felder unterteilt und 
mit der Angabe der je Feld beobachteten Fleckzahlen bezeichnet. In den kleineren 
Randfeldern bedeuten die Zahlenwerte die auf die als Flächeneinheit anzusehenden 
Rechteckfelder bezogenen Punktzahlen. Die Gesamtzahl der bei dieser Messung 
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ausgezählten Punkte beträgt N = 7223. Durch die in Abb. 10 strichpunktiert 
eingezeichnete Symmetrieachse wird die Kathode in zwei spiegelbildliche Hälften 
unterteilt, deren Fleckzahlen von 3730 und 3493 innerhalb der statistischen Fehler- 
grenzen übereinstimmen. Als Beispiel für die Einteilung der Kathodenfläche in 
Zonen, innerhalb derer die Punktdichte n zwischen nm und n + An liegt, ist eine 
Zone mit n = 30 und An = 15 durch stärkere Umrandung gekennzeichnet. Die 
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81 | | 52) 46) | 1B | 12) \ 
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- 
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10. Fleckdichte-Verteilung für 5,0 A 


von dieser Zone bedeckte Fläche f beträgt 32,5 Einheiten und die Anzahl der 


Punkte nf in der Zone 1207. n ergibt sich zu oe = 37,2 Punkte/Flächen- 
einheit und w(n)-An berechnet sich nach Gl. (3) zu 0,167. Nachstehende Zahlen- 
tafel 2 enthält in den ersten drei Spalten die für alle Zonen ermittelten Werte 7, 


Zahlentafel 2 


7 w(n) An 2 w(n) An Uw v 
0 
7,6 0,062 0,062 | 2,8 1,05 - 10-3 

22,3 0,188 0,250 1,25 3,09 - 10-3 ‘4 
37,2 0,167 0,417 0,45 5,15 - 10-3 xi 
51,9 0,189 0,606 —0,53 7,18 - 10-3 
66,8 0,168 0,774 —1,4 9,24 - 10-3 
82,0 0,136 0,910 —2,4 1,13 - 10-2 
93,9 0,090 1,000 —4,5 1,30 - 10-2 - 

1,000 
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= 258 BE. 
w(n) An und die für den Anschluß an die Messung der Brennspannungsverteilung U’. 
n 
zu berechnende Summe N w(n) An. bei 
0 
Eine Betrachtung der Zahlenwerte in Abb. 10 vermittelt wegen der bei kleinen span 

Punktzahlen erheblichen statistischen Streuungen von n nur einen rohen Eindruck Sa 

von der Häufigkeitsverteilung der Fleckdichte über den Kathodenquerschnitt. 

In Abb. 11 ist nun versucht worden, diese Schwankungen durch ein Interpolations- erha 
den 
relat 
letz: 

für ı 
9 A — 
3% 
hag 
= 
= aa 
Ku 

Abb. 11. Fleckdichte-Verteilung nach kreuzweiser Mittelung der Meßwerte der Abb. 10 Kol 
mit 

verfahren auszugleichen, um die Kurvenscharen rn = const einzeichnen zu können. in f 

Hierbei wurde folgendermaßen vorgegangen: Trägt man die jeweils in einer senk- 

rechten Spalte der Abb. 10 stehenden Zahlen in Abhängigkeit von ihrem Abstand = 

von der Symmetrieachse in ein Koordinatensystem ein, so kann man unter Berück- > | 

sichtigung der Symmetrie durch diese Punktfolge eine abgleichende Kurve legen, 
aus der für jedes Kästchen ein verbesserter Zahlenwert für n abgelesen wird. Ganz Abl 

ähnlich verfährt man für die waagerechten Spalten. Damit erhält man für jedes i. 

Kästchen eine zweiten Zahlenwert, der mit dem ersten gemittelt wird. Die auf 2 

diese Weise entstehenden n-Werte sind in der oberen Hälfte des Kathodenkreises * 5 

in Abb. 11 dargestellt. Da sie sich diesmal monoton von Kästchen zu Kästchen te 

ändern, lassen sich aus ihnen bequem die Kurvenscharen n = const konstruieren, Bez 
die in der unteren Bildhalfte eingezeichnet sind. ee 
Die zum Entladungsstrom von 5,0 A gehörige a nn aan, wird ied 

Die 

in Abb. 12 w wiedergegeben. Aufgetragen.ist — = vor 
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ilung U’. Für jeden der Meßpunkte beträgt Z = 500. Die Kurve schneidet die Ordinate 


bei Z — 0,5. Das bedeutet, daß die Zählschaltung die Momentanwerte der Brenn- 


einen spannung einwandfrei mit dem arithmetischen Mittelwert vergleicht. Durch Gleich- 


-. setzung der Werte By w(n) An in Zahlentafel 2 mit rn entsprechend Gl. (2) und (5) 
0 


‘lons- erhält man aus der Kurve die zu jedem  gehérigen U’-Werte der Spalte 4. Um 
den von der Anzahl der Beobachtungen abhängigen Wert 7 in hiervon freien, 


relativen Einheiten darzustellen, bilden wir n =v (Spalte 5). Damit stellen die 


letzten beiden Spalten Jas Ergebnis dieser Untersuchung dar. Abb. 13 gibt die 
für unser Beispiel mit / = 5 A und die für 3,5 und 8 A in gleicher Weise gemessenen 
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Abb. 12. Brennspannungsverteilung Abb.13. Zusammenhang zwischen Aufent- 
für 5,0 A haltswahrscheinlichkeit und Brennspannung 


Kurven » (U’) wieder. Die in der 5 A-Kurve stärker markierten Meßpunkte sind 
Kontrollwerte, die durch eine direkte Messung des Zusammenhanges n (U’) er- 
mittelt wurden. Dies war durch eine einfache Abänderung der Versuchsanordnung 
nen, in folgender Weise möglich: 


Mit Hilfe der in Abb. 7 angedeuteten Beleuchtungsvorrichtung fällt bei Frei- 
aii gabe des Objektivs ein Lichtimpuls auf die in der Kamera befindliche Photozelle. 
Diese arbeitet an Stelle des Kontaktapparates und des Steuersatzes iiber einen ein- 
stufigen Verstärker auf das Steuergitter des Stromtores // der Schaltung in 
rw. Abb. 9. Damit sind beide Apparaturen gekuppelt. Mit Hilfe eines Kopfhörers, 
der in den Stromkreis des Zählwerkes 1 geschaltet wurde, konnte beim Aufleuchten 


auf 

Pen: des Brennflecks auf der Mattscheibe festgestellt werden, ob die zugehörige Brenn- 
wee spannung den durch die Kompensationsspannung gegebenen Wert iiberschreitet. 
ee: Hierdurch war eine Nachprüfung der durch statistische Zuordnung gefundenen 
’ 


Beziehung n (U’) durch direkte Messung ermöglicht. Eine saubere Zuordnung von 
ind n und U’ über den ganzen Kathodenquerschnitt erfordert bei diesem Verfahren 
jedoch eine so große Anzahl von Meßreihen, daß hierauf verzichtet werden mußte. 
Die Grenzkurven für n und U’ werden nach dem Kathodenrande hin unscharf. 
Infolgedessen wurden die an der 5 A-Kurve vermerkten Werte nur durch Berück- 
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sichtigung der in der Nahe der Symmetrieachse liegenden Punkte bestimmt. Sie 
liegen innerhalb der Meßgenauigkeit befriedigend auf der statistisch gefundenen 
Kurve. 

Die Tatsache, daß die Meßpunkte in der logarithmischen Darstellung nicht 
auf einer Geraden liegen, besagt, daß das Potential m der Gl. (1) in keinem linearen 
Zusammenhang mit der Brennspannung stehen kann. Es gelingt nun zwar, die 
Meßergebnisse für konstanten Strom beispielsweise durch eine Gleichung der 
Form » = 1% e-kı(U - Vo)" genügend genau wiederzugeben, in der », den zu einem be- 
liebigen Bezugswert U, gehörigen v-Wert und k, und k, Konstanten bedeuten, doch 
bleiben dabei k, und k, vom Bezugspunkt abhängig, was für eine Boltzmann- 
Verteilung nicht der Fall sein darf. Ohne den Meßwerten Zwang anzutun, schien 
es vorerst nicht möglich, eine der Boltzmann-Verteilung entsprechende Dar- 
stellung zu finden. Da ferner die Annahme einer über den Kathodenquerschnitt 
konstanten Unruheenergie Z wegen eines denkbaren Einflusses der je nach der Lage 
des Brennflecks unter der Glasscheibe verschiedenen Dampfströmungsverhältnisse 
auf die Bewegungsunruhe nicht sicher zu begründen war, wurden die Versuche 
mit einer anderen Anordnung fortgesetzt. Über den Einfluß des Entladungsstromes 
sei an dieser Stelle nur festgestellt, daß mit zunehmendem Strom die Kurven 
v (U’) steiler verlaufen, d. h. daß die Fleckhäufigkeit nach den Gebieten steigender 
Brennspannung zu stärker absinkt als bei geringerem Strom. ee 


d) Versuchsanordnung II (lineare Anordnung) fs = q 
In der zweiten Versuchsanordnung fand das in Abb. 14 dargestellte Entladungs- 
gefaB Verwendung. Das horizontal liegende zylindrische Glasrohr von 1 m Linge 
enthält als Kathode längs einer Mantellinie einen etwa 4,5 cm breiten Quecksilber- 
streifen, auf dem der Brenn- 


HAünIwasser 4 fleck praktisch nur in 


= Achsenrichtung umherlaufen 


kann. Anoden- und Katho- 
denzuführung befinden sich 
an demselben Rohrende. Zur 
Festlegung des Quecksilber- 
dampfdruckes wird der Kühl- 
mantel mittels Umwälz- 

| Phorozelle | 7 pumpe und Rückkühler von 
Wasser konstanter Tempe- 
ratur durchströmt. Ein Luft- 
strom verhindert bei höhe- 
ren Wassertemperaturen das 
Weichwerden der Versuchsgefäß und Wassermantel verbindenden Pizeinkittung. 
Der Quecksilberdampfdruck wird hierdurch nicht beeinflußt; denn die Temperatur 
der Kondensationsgrenze am nichtwassergekühlten Gefäßende war stets in Über- 
einstimmung mit der Wassertemperatur. Während in der Anordnung I die Brenn- 
spannung von den zwei Ortskoordinaten (x, y) des Brennflecks abhing, ist hier die 
Brennspannung lediglich eine Funktion von x. Wegen des konstanten Säulen- 
gradienten ist diese linear und leicht zu ermitteln, so daß als wesentliche Aufgabe 
nur die Messung der Häufigkeitsverteilung des Brennflecks in Richtung der Rohr- 
achse übrigbleibt. Dies geschieht mit Hilfe einer Photozelle (vgl. Abb. 14), die, längs 
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Abb. 14. Versuchsaufbau II 
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‘des Rohres verschiebbar angeordnet, immer dann vom Licht der positiven Säule 
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getroffen wird, wenn sich der Brennfleck rechts von ihr befindet. Die Photozelle 
arbeitet wieder über ein Verstärkerrohr auf das eine Gitter eines ein Zählwerk be- 
tätigendes (Zählangabe = z) Doppelgitterstromtores, dessen zweites Gitter von 
einer Stromtorkippschaltung in gleichmäßigen Zeitabständen mit kurzzeitigen 
Spannungsimpulsen positiv gemacht wird. Bedeutet Z die Anzahl der Spannungs- 


impulse während der gesamten Meßzeit, so gibt si die Wahrscheinlichkeit dafiir 


an, daß sich der Brennfleck rechts von der an der Stelle x befindlichen Photozelle 
aufhält. 

Die Einzelheiten der Schaltung sind in Abb. 15 dargestellt. Bei Belichtung der 
Zelle sinkt der Ruhestrom der Verstärkgrröhre von 5mA praktisch auf Null. 
Damit wird das eine © -. 
Gitterdes Doppelgitter- 
stromtores freigegeben. 
Die Kippschaltung zur 
Erzeugung der Span- 
nungsimpulse fiir das 
zweite Gitter wirkt fol- 
gendermaßen : Durch die 
Aufteilung der zwischen 
Anode und Kathode des 
Eingitterstromtores lie- 
genden Wechselspan- 
nung über einen ohm- 


schen Spannungsteiler Gent 
liegt am Gitter eine Impulsschaltung 


kleine Wechselspan- | Zus, Abb. 15. Versuchsschaltung II 
nung, die das Rohr zum 

Zünden bringt. Der auf den Kondensator von 1 uv F kurzzeitig fließende Ladestrom 
erzeugt am Kathodenwiderstand von 1 kQ einen Spannungsabfall, der auf das zweite 
Gitter des Doppelgitterstromtores wirkt. Eine Neuzündung des Eingitterrohres 
kann erst dann erfolgen, wenn sich der Kondensator über den die Kippfrequenz 
bestimmenden Widerstand von 200 kQ nahezu entladen hat. Durch das Anlegen 
des Gitters an Wechselspannung zündet das Rohr im Scheitel der Wechselspannung. 
Damit ist die Höhe der Spannungsimpulse genügend unempfindlich gegen geringe 
Zündwinkelschwankungen. Überraschend genau läßt sich die Kippfrequenz ein- 
halten. Sie betrug über mehrere Wochen ohne jegliche Veränderungen in der 
Schaltung 150 Impulse/min mit einer Schwankung um + 1%, so daß an Stelle eines 
Zählwerkes für Z die Uhr verwendet werden konnte. Über die Ansprechempfind- 
lichkeit der Zählschaltung ist noch zu bemerken, daß eine Registrierung immer erst 
dann erfolgt, wenn der Brennfleck bereits so weit nach rechts gelaufen ist, daß das 
Licht der Säule auf die Photozelle fällt, die je nach Entladungsstrom und Druck 
erst in einigen cm Abstand hinter dem Brennfleck beginnt. Wie die Meßergebnisse 
zeigen werden, ist dieser Umstand ohne Bedeutung, da der Nullpunkt für die Zäh- 
lung von x vollkommen willkürlich gewählt werden kann. Störendes Streu- und 
Nebenlicht wurde durch geeignete Verkleidung von Photozelle und Versuchsgefäß 
ferngehalten. Die Zündung des Versuchsgefäßes erfolgt am ungekühlten Gefäß- 
ende durch einen von außen in Höhe des Quecksilbermeniskus angebrachten 
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leitenden Belag, der mit einem hochfrequenten Spmmngened beschickt wird, 
Durch Verwendung eines Thermostaten für das Kühlwasser, einer selbsttätigen 
Zündung und verschiedener Sicherheitseinrichtungen konnte die Apparatur ohne 
Beaufsichtigung Tag und Nacht laufen. Von Zeit zu Zeit wurde die Photozelle 


versetzt und so für gegebenen Strom und Dampfdruck . in Abhängigkeit von x 


mit hoher Genauigkeit gemessen. Die zugehörigen Brennspannungswerte wurden 
in Sonderversuchen teils oszillographisch, teils mit einem gewöhnlichen Dreh- 
spulinstrument bgim Durchlaufen angebrachter Meßmarken ermittelt. Als trei- 
bende Spannung für die Entladung wurden wieder 440 V verwandt, um den Strom 
möglichst unabhängig von der Höhe der Brennspannung konstant zu halten. 


® 
e) Ergebnisse der Messungen an Versuchsanordnung II 


In Abb. 16 ist die für die drei verschiedenen Entladungsströme bei einem 
konstant N Quecksilberdampfdruck von 4,3 - 10-3 Torr gemessene Ab- 


hängigkeit > (x) logarithmisch dargestellt. Die Meßwerte liegen innerhalb von 
2 Zehnerpotenzen sehr gut auf geraden 


er: - Linien. Da die Aufenthaltswahrscheinlich- 
~ u keit des Brennflecks namentlich für den 
1 N Br Strom von 8,16 A in den anodenfernen 

Gebieten sehr gering ist, waren fiir jeden 
1 \ N Meßpunkt mehrere Beobachtungstage er- 
0" N N forderlich. Auf die sonst bei statistischen 
1 \ N Messungen meist durchgeführte Voraus- 
Aa 123% A berechnung des mittleren Meßfehlers aus 
=. N der Anzahl der Registrierungen mußte 
q | \ verzichtet werden, weil hier die Zählungen 
x \ \ nicht nachwirkungsfrei sind. Befindet sich 
o'r X nämlich der Brennfleck einmal rechts 
s- von der Photozelle, und spricht das Zähl- 
se \ werk an, so ist wegen der im Mittel ge- 
ringen Laufgeschwindigkeit des Flecks 
T IN die Wahrscheinlichkeit, dafür, daß sich bei 
o der nächsten Freigabe der Zählschaltung 
& \ der Fleck immer noch rechts von der 
au = Zelle aufhält, erheblich größer, als wenn 
# I 1=816 A er erst von links in dieses Gebiet hinein- 
Zr 3 laufen müßte. Eine einwandfreie Voraus- 
a Er berechnung des statistischen Fehlers aus 
der Zahl der Registrierungen hätte einen 
Abb. 16. Ortliche Häufigkeitsverteilung zeitlichen Abstand zwischen den Meß- 
des Brennflecks impulsen erfordert, der gegenüber der zu- 


nächst unbekannten Nachwirkungsdauer 
(Größenordnung 1 min) groß sein müßte. Die geringe Streuung der MeBpunkte 
beweist jedoch, daß die Anzahl der Beobachtungen ausreichend war, und daß 
wir die gefundenen Werte darstellen können durch einen Ausdruck der Form 
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Da andererseits au die Wahrscheinlichkeit dafür angibt, daß sich der Brennfleck 


im Gebiet rechts der Marke x, d. h. im Bereich x... oo aufhält, gilt 
== | n(x)de, (7) 


x 


worin n wiederum die Fleckdichte und N die Gesamtzahl der Flecke bedeutet. 
Durch Differentiation beider Gleichungen nach x erhalten wir 


— Cd, p)g(l, p) — (a) 
oder | 
n(x) = (8) 
Da der Bezugspunkt für die Zählung von x auf dieses Ergebnis keinerlei Einfluß 


hat, können wir, diesmal die gefundene Abhängigkeit zum Boltzmann-Satz in 
Beziehung setzen. Demnach müßte g x gleich dem Verhältnis von potentieller 


zu Unordnungsenergie 2 sein (vgl. Gl. (1)). In nachstehender Zahlentafel 3 


sind die aus den Messungen berechneten g- Werte in Abhängigkeit vom Entladungs- 
strom J und Quecksilberdampfdruck p dargestellt. Die Zusammenstellung ent- 
hält außerdem die zugehörigen, aus Brennspannungsmessungen an verschiedenen 
Stellen x bestimmten Spannungsgradienten €. 


Zahlentafel 3 


p g | € 
3,94 A 4,3 - 1073 Torr | 5,36 - 10”? cm! | 0,2 er: 
5,6 - 10-4 1,28 - 10-1 '01 
D 571A | 1,6 - 1073 | 1,00 - 1071 0,1 188 
be - 1073 | 0,94 + 107% 0,246 
1,16 - 10”? 1,02 - 1071 | 0,238 
r 8,16A  4,3- 10-3 1,89 - 1071 | 0242 


Man erkennt, daß g praktisch nur von J abhängt. Bei einem Entladungs- 
strom von 5,71 A ist g in einem Druckbereich von 1,6 - 10-3... 1,16 - 102 Torr 
konstant. Erst bei 5,6 - 10-4 Torr wird ein etwas höherer Wert gefunden. Bevor 
wir die Ergebnisse zu deuten versuchen, sei zunächst noch eine weitere Meßmethode 
mitgeteilt, mit der die gefundenen g-Werte durch ein vom bisherigen völlig ver- 
schiedenes Verfahren geprüft werden konnten. 


IV. Untersuchung der Brennfleckdiffusion 


Wie die unter IIb beschriebenen Messungen gezeigt haben, ist das mittlere 


Verschiebungsquadrat der zufälligen Brennfleckwanderung Az? proportional dem 
zugehörigen Beobachtungszeitintervall At. Für die Proportionalitätskonstante 
liefert die Theorie der Schwankungserscheinungen ?) den Faktor 2D, d.h. 


7) Siehe z. B. G. Sub, Lehrbuch der theoretischen Ph: Physik (1932), S. 496. 
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D bedeutet hierin die Diffusionskonstante des betrachteten Teilchens im ther- 
misch bewegten Medium, die sich aus dem Anstieg der gemessenen Geraden der 
- Abb. 3 zu 8,0 em? - s-1 berechnet. Gl. (9) gilt nur für den Fall, daß keine Richt- 
kraft auf den Fleck wirkt. Wir nehmen nun an, daß sich die durch die Richtkraft K 
erzeugte Drift v, durch den Beweglichkeitsansatz der Form 


u=b-K (10) 


beschreiben läßt, worin 6 die nur von Strom und Druck abhängige Beweglichkeit 
bedeutet. Zwischen Diffusionskonstante, Beweglichkeit und thermischer Energie 
besteht folgender allgemeingültiger Zusammenhang’): 


D=b.-E. (11) 
Da nach Gl.(1) und (8) „-7 x ist, andererseits X = 7 "ach Gl. (10) und (11) durch 5 


gegeben ist, bietet die ins der Drift und ie Diffusionskönstanten eine Mög. 
lichkeit, die gefundenen g- Werte zu überprüfen. v, kann in einfacher Weise dadurch 
bestimmt werden, daß der Brennfleck mit Hilfe eines Magneten bis zum Rohr- 
ende abgelenkt wird. Läßt man ihn unbehindert zurücklaufen und stoppt die 
zum Durchlaufen von angebrachten Meßmarken gehörigen Zeiten ab, so ergibt sich 


v, als Mittelwert von = . Zur Ermittlung von D muß Gl. (9) für den Fall der 


Drift korrigiert werden. Hierzu dient folgende Überlegung: Der Ort x des Teil- 
chens werde alle t s gemessen. Zu jedem Zeitintervall 7, das, wie bereits 
bemerkt, so groß sein muß, daß innerhalb dieser Zeit viele Umlenkungen erfolgen 
können, gehört eine Verschiebung in x-Richtung Az, die ihrem Vorzeichen nach 
positiv oder negativ sein kann. Az setzt sich zusammen aus einem ungeordneten 
Anteil £ und dem Driftanteil v, 


Ax = € UT . 
Dann beträgt nach z Beobachtungen, die in der Zeit zr erfolgen, i: Br 
= Ax,= SE, +27. (12) 
Durch Quadrieren ergibt sich: 


Axt. = ze +25 E, 


Im Mittel werden die Summen des zweiten und dritten Gliedes wegen der 
gleichen Häufigkeit positiven und negativen Vorzeichens null, d.h. 


oder, wenn wir zt = At setzen und das zugehörige Verschiebungsquadrat cinta 
mit Az? bezeichnen, 


£2 
Az? = = At + At. 


£ ist aber identisch mit dem Az?/At der Gl. (9) und gleich 2 Dm RN Wir 
erhalten damit die gesuchte Beziehung > 
4 Az? =2D-At+v?-At. (13) 
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Durch Mittelbildung in Gl. (12) folgt für &,, wie bereits vorweggenommen, 

Der Gang der Untersuchung ist nun folgender: Nach Bestimmung von 1% 
werden für verschiedene, während einer Meßreihe konstante At wie in IIb die 
Verschiebungen Ax gemessen, Dies kann wegen der zur Verfügung stehenden 
großen Beobachtungsstrecken mit Stoppuhr und Längenmaßstab erfolgen. Für 
kurze Zeiten werden die Ax-Werte aus dem bekannten Spannungsgradienten über 
die oszillographierte Brennspannung errechnet. Nach Berechnung von Ax? wird 
Ax®—v} - At? in Abhängigkeit von At linear aufgetragen. Der Anstieg der von 
der Theorie geforderten Geraden liefert die ge- cm: 


Abb. 17 gibt als Beispiel die für einen Ent- | =’ 

ladungsstrom von 5,71 A und für einen Druck ™™ 

von 4,3 - 10-® Torr gemessenen Werte wieder. ©, Ya 

beträgt 0,84 cm -s~1, die Meßpunkte liegen be- ™ 


friedigend auf einer Geraden, und wir erhalten 
für D 12,8cm?-s“. Dieser Wert stimmt mit 
dem am Fleck ohne Richtkraft durch Auswertung 


| /=57A p= 43-107 Torr 


der Filmaufnahmen bei 5,5 A gefundenen Anstieg so wen 

(vgl. Abb. 3), aus dem sich D zu 8 cm? - s~ be- i 

rechnet, in Anbetracht der vielen Fehlermöglich- 0 

ahlentafel 4 sin | nun die für drei verschiedene des mittleren Verschiebungs- 

Stromstärken bei einem Druck von 4,3 - 10-3 quadrates 


Torr auftretenden v,/D-Werte den unter gleichen 

Bedingungen beobachteten g-Werten gegenübergestellt. Die Übereinstimmung ist 
zwar nicht vollständig, doch bleibt das Verhältnis von g zu v,/D mit etwa 1,4 viel 
besser konstant, als man hätte erwarten können. Die gefundenen g-Werte scheinen 
damit genügend gesichert. 


Zahlentafel 4 


| g 
I | | D | | 9 
3,94 A 0,33 cm + 8,6 cm? - 3,83 - 102 5,36 - 10”? em! L 1,40 
5,71 0,84 6,56 - 10-2 9,40 - 10-2 1,43 
8,16 1,98 14,0 | 1,41 - 1071 1,89 - 10-1 1,34 — 


V. Versuch einer Deutung der bisherigen Ergebnisse 


Die Verwendung des Minimumprinzips als Arbeitshypothese legt den Gedanken 
nahe, folgende Deutung der Meßergebnisse zu versuchen. Bei Verlängerung der 
Säule um die Längeneinheit ist eine Arbeit zu leisten, die durch das Leistungs- 
produkt € - J gegeben ist. Aus Dimensionsgriinden ist diesem Produkt eine Zeit r 
hinzuzufügen, der man vielleicht die Bedeutung der mittleren Trägerlebensdauer 
beimessen könnte. Damit wird € - J - ein Maß für die bei der Säulenverlängerung 
zu leistende Ionisierungsarbeit. Durch die Vergrößerung der Leiterschleife wächst: 
aber ferner der Energieinhalt des magnetischen Feldes. Die zu leistende Arbeit 
beträgt im Sinne einer Blaswirkung auf die Leiterschleife nach außen —!/, dL/dz - I®, 
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fiir Versuchsanordnung II also — const - I. = Der Boltzmann-Ansatz müßte / 
demnach lauten 


Nimmt man nun versuchsweise an, daß E nur wenig von J abhängt, so müßte 
sein 
= ¢,:T—c,:J/E, mit c, und 
Wie hier nicht näher ausgeführt, ergab die Nachprüfung dieser Gleichung keine 
Übereinstimmung mit den Meßwerten. Das erste Glied zeigt trotz des zu erwartenden 
starken Druckeinflusses auf die Trägerlebensdauer keinerlei Druckabhängigkeit. 
Das zweite hatte ein anderes Vorzeichen, so daß nach einer anderen Deutungs- 
möglichkeit gesucht werden mußte. 


Das falsche Vorzeichen des zweiten Gliedes obiger Gleichung läßt vermuten, 
daß hier der von J. Stark®) an einer Quecksilberdampf-Niederdruckentladung 
entdeckte anomale Einfluß von transversalen Magnetfeldern auf den Brennfleck 
eine Rolle spielen könnte. Dabei handelt es sich um eine Ablenkung des Brenn- 
flecks in eine Richtung, die der nach elektrodynamischen Überlegungen zu erwarten- 
den Auslenkung entgegengesetzt gerichtet ist. Dagegen wird die positive Säule 
gewöhnlich im elektrodynamisch richtigen Sinne abgelenkt. Die Erscheinung ist 
später von mehreren Verfassern ®) untersucht worden, ohne daß bisher eine völlige 
Klärung der Verhältnisse erreicht werden konnte. Insbesondere liegen meines 
Wissens keine Angaben darüber vor, wie die in elektrodynamisch falschem Sinne 
auf den Brennfleck wirkende Kraft mit der magnetischen Feldstärke zahlenmäßig 
zusammenhängt. Die nachstehend beschriebenen Untersuchungen sollen die Frage 
prüfen, ob das magnetische Eigenfeld des Entladungsstromes unter Berücksichti- 
gung des anomalen Effektes etwa ausreicht, um die beobachteten Richtkräfte 
zu deuten. 


VI. Untersuchung des Einflusses von transversalen Magnetfeldern 
auf die Driftgeschwindigkeit des Brennflecks 


Das Versuchsgefäß der Anordnung II wurde in ein von einem Helmholtz- 
Spulensatz erzeugtes, homogenes Magnetfeld von einigen Oersted gebracht, das 
senkrecht zur Gefäßachse und parallel zur Quecksilberoberfläche gerichtet war 
und die ganze Gefäßlänge erfaßte. Die Auslenkung der positiven Säule an die 
Gefäßwand erfolgte unabhängig von Druck und Stromstärke stets im elektrody- 
namisch zu erwartenden Sinne, während die Driftgeschwindigkeit des Brennflecks 
bei einem die Zugwirkung der Entladung genügend überkompensierenden Feld 
elektrodynamisch falsch gerichtet war. 


8) J. Stark, Induktionserscheinungen am Quecksilberlichtbogen im Magnetfeld. 
Physik. Z. 4, 440 (1903). 

*) E. Weintraub, Investigation of the Arc in Metallic Vapours on an Exhausted 
Space. Philos. Mag. 7, 95 (1904); M. N. Minorsky, La rotation de l’are électrique dans 
un champ magnétique radial. J. Physique 9, 127 (1928); R. Tanberg, Motion of an 
Electric Arc in a Magnetic Field under Low Gas Pressure. Nature 124, 371 (1929); U. 
Tiberio, Die Wirkung des magnetischen Feldes auf den elektrischen Bogen im Hg- 
Dampf. (Text italienisch) L’Elettroteenica 17, 485 (1930). rte ae 
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Um nun den Einfluß des Eigenfeldes zu ermitteln, wurde die Driftgeschwindig- 
keit v in Abhängigkeit vom Zusatzmagnetfeld H untersucht. Die Vorzeichen on 
durch folgende Regeln bestimmt: 


Wenn Feld so gerichtet wie Eigenfeld des Entladungsstromes: 7 >0, 
wenn Drift so gerichtet wie Eigendrift des Brennflecks: v> 0, 


wenn Kraft so gerichtet wie Eigendrift des Brennflecks: K>0. 
Als MeBbeispiel ist in Abb. 18 v (#) fiir einen Strom von 5,8 A und einen Druck 


von 2,5-10-* Torr wiedergegeben. 


raden der Form 


Dieser Zusammenhang erweist 
sich im untersuchten Bereich von 
1,5: 103 bis 2,5 - 10? Torr als 
druckunabhängig. Dagegen sind 
v, und c, wie aus Abb. 19 und 20 
zu entnehmen ist, stark von der 
Stromstärke abhängig. Die Eigen- 
drift v, ist mit dem v, der Gl. (14) 
identisch. Es ist nun zwar nicht 
gelungen, direkt von den Meß- 
werten ausgehend, etwa unter 
Verwendung eines parabolischen 
Ansatzes für v, (/) nach Abb. 19 
und eines linearen Ansatzes für 
c(I) nach Abb. 20, einen ver- 
ständlichen Zusammenhang die- 
ser Größen zu erhalten. Doch 


Die Meßpunkte liegen sehr gut auf einer Ge- 


v=u+cdH. (15) 


8cm-s" 


V 


n 


7=58A 
p=25-107? Torr 


/ —H 


-4 -2 Fi 0 4 60ersted 


% 4 
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lassen sich die Ergebnisse unter £ 


Zugrundelegung folgender Vor- 


stellung in ein Bild einordnen, 
das den Erfahrungen nicht wider- 
spricht: Der Beweglichkeitsansatz der Gl. (10) ist auf jeden Fall zulässig, weil 
die Driftgeschwindigkeit über die ganze Rohrlänge konstant ist. Die Beweglich- 
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Abb. 19. Stromabhängigkeit der Eigendrift- 


geschwindigkeit 


Abb. 18. Driftgeschwindigkeit im transversalen 


Magnetfeld 


keit b ist dabei stromabhängig. Unter 
völligem Verzicht auf das Minimum- 
prinzip setzen wir nun 


v()=b(I)-(«P+PßIH). (16) 

Das erste Glied in der Klammer 
stellt dabei die durch das Eigenfeld 
hervorgerufene Kraftwirkung dar, die 
hier so behandelt wird, als erfolge sie 
im Sinne der Elektrodynamik, aber 


mit umgekehrten Vorzeichen. Für 


die Zulässigkeit dieses Ansatzes inner- 
halb des Untersuchungsbereiches 
spricht das lineare Ansteigen von v 
mit H. x, dem die Bedeutung von 
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'1/,-dL/dx zukommt, wird nur ‘in erster Näherung als konstant anzusehen sein, 
weil mit zunehmendem Magnetfeld die Säule an die obere Gefäßwand gedrückt 
x wird und weil damit dL/dx wegen des wachsenden Abstandes von Säule zu 
Kathodenrückleitung etwas zunehmen muß. Als zweiter Anteil erscheint in 
der Klammer ein Glied, das die Kraftwirkung des überlagerten Zusatzfeldes 
wiedergibt. Dabei bedeutet ß eine Konstante, die durch die Gestalt der Ent- 
_ladungsbahn in Brennflecknähe gegeben ist. Durch Vergleich mit Gl. (15) er- 
halten wir 
(1) = b (I) « I? (17) 16 Dersted - 
c(I)=b(I)ßI. (18) 4 
Division beider Gleichungen liefert 7 
= (19) << 
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Abb. 20. Stromabhängigkeit von c= 77 _ Stromabhängigkeit der Beweglichkeit 


In Abb. 21 ist nun v (Z)/c (I) J in Abhängigkeit von J aufgetragen. Die Streu- 
ung der aus den Meßpunkten berechneten Werte rührt hauptsächlich von der 
Unsicherheit der Bestimmung von v, her. Der eingezeichnete Kurvenzug ist aus 
den ausgleichenden Kurven der Abb. 19 und 20 berechnet. Wie zu erwarten, steigt 
a/ß mit wachsendem J mit Sättigungscharakter wegen der Zunahme von d.L/dx 
an. Nehmen wir ß als praktisch unveränderlich an, so können wir nach Gl. (18) 
die Stromabhängigkeit der Beweglichkeit angeben: 


Die hiernach ermittelten Werte 8 5 sind in Abb. 22 in Abhängigkeit von J darge- 
stellt. Man erkennt eine wegen des wachsenden Fleckdurchmessers bzw. beginnen- 
der Fleckaufteilung zu erwartende Abnahme der Beweglichkeit mit steigendem 
Strom. 

Die Kenntnis der Stromabhängigkeit der Beweglichkeit gestattet über die 
Messung der Diffusionskonstanten eine Berechnung der Stromabhängigkeit der 
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Unordnungsenergie E nach Gl. (11). D(Z) wurde wie bisher durch Ermittlung 4 
der mittleren Verschiebungsquadrate bestimmt. Die Driftgeschwindigkeit brauchte a 
jedoch diesmal nicht beriicksichtigt zu werden, weil ein Magnetfeld verwandt ae 

wurde (H = — v,/c), das v, kompensierte. Abb. 23 gibt die gemessenen Werte D (I) 2: 
und Abb. 24 das der Unordnungsenergie proportionale Verhältnis DJ/c wieder. “ 

Daß E mit J ansteigt, erscheint wiederum verständlich, wenn wir die Unruhe- = 
bewegung auf den heftigen Verdampfungsvorgang zurückführen. 4 
Die Größe v,/c I läßt sich auch nach Gl. (15) und (16) als das von 1 A erzeugte, u 
auf den Brennfleck wirkende Eigenfeld des Entladungsstromes deuten. Es liegt 
im Größenbereich von 0,1 Oe. Ein gleiches Feld wird aber in der Mitte einer langen 
Leiterschleife mit 4cm Mittelpunktabstand, d.h. dem halben Rohrradius als 
Schleifenweite erzeugt. (Das magnetische Erdfeld hat bei den Untersuchungen 
keine Rolle gespielt, da die Horizontalkomponente genau senkrecht zu H stand.) 
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Zusammenfassend läßt sich also sagen, daß mit der zugrunde gelegten Vorstel- a 
lung die Verhältnisse in nicht unverniinftiger Weise beschrieben werden können. 
Es besteht nun noch die Frage, ob es einen Sinn hat, die zuletzt völlig außer acht 
gelassene Forderung des Minimumprinzips zu berücksichtigen. Die von M. Steen- 
beck gefundene Bogenzugspannung im Lichtbogen soll dem Zug entsprechen, um 
welchen bei starken Elektrolyten der osmotische Zug infolge der elektrostatischen as 
Kräfte verringert wird. Für die Zugspannungen im nichtisothermen Plasma er- 


konstanten 


€ 70 
Abb. 23. Stromabhängigkeit der Diffusions- 


4 6 


Abb. 24. Stromabhängigkeit der Unruhe- 
energie des Brennfl 


gibt sich nach R. Rompe und M. Steenbeck”): 


Da bei uns die mittlere Trägerkonzentration N höchstens 101? cm? beträgt, 
erhalten wir mit einer Elektronentemperatur 7’ von schätzungsweise 15000° K 
eine Zugspannung von 1 - 10-* dyn cm”, d.h. für die Säule mit etwa 50 cm? 
Querschnitt eine Zugkraft von 5 - 10-2? dyn. Dagegen beträgt die vom magneti- 
schen Eigenfeld auf den senkrecht zum Brennfleck gerichteten Teil der Entladungs- 
bahn ausgeübte Kraft bei 6 A ungefähr 1 dyn. Bedenkt man, daß diese nach außen 


10) R.Rompeu.M.Steenbeck, Der Plasmazustand der Gase. Ergebn. exakt. Naturw. 
Bd. XVIII, 257 (1939). 
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gerichtete dehneipanindtn Kraft durch den anomalen Magnetfeldeffekt soweit 
iiberkompensiert wird, daB eine Bewegung nach der entgegengesetzten Seite 
zustande kommt, so muß man folgern, daß der Kohäsionszug demgegenüber gar 
keine Rolle spielt. 

Der starke Einfluß des Magnetfeldes des Entladungsstromes macht es verständ- 
lich, weshalb sich die mit Versuchsanordnung I gefundenen Ergebnsisse nicht durch 
den Boltzmann-Satz beschreiben lassen. Wahrscheinlich ist die Forderung eines 
auf der ganzen Kathode eindeutigen Zusammenhanges zwischen Fleckdichte und 
Brennspannung, auf der das statistische Meßverfahren beruhte, überhaupt nicht 
zulässig. Die annähernde Übereinstimmung der statistisch gefundenen Meßwerte 
mit den lediglich auf der Symmetrieachse der Kathode durch direkte Messungen 
kontrollierten Werten mag ein Zufall sein. Eine nachträgliche, rechnerische 
Berücksichtigung des Eigenfeldes ist wegen der komplizierten Geometrie und 
der Nähe von Eisenteilen nur sehr schwer durchführbar. Die stärkere Abnahme der 

* Fleckdichte bei größerem Strom erfolgt jedenfalls im Sinne unserer Überlegung. 
Gasdruckänderungen, deren Einfluß als unwesentlich erkannt wurden, können sich 
auf die Vorgänge nur insofern auswirken, als sie die Stromdichteverteilung und 
damit das magnetische Feld der Entladung im geringem Maße beeinflussen. (Bei 
erheblich höheren Drucken erfolgt allerdings nach M. N. Minorsky®) eine Um- 
kehr des anomalen Magnetfeldeffektes zu elektrodynamisch richtigem Verhalten.) 

Die vorliegenden Untersuchungen gestatteten lediglich die Bestimmung von 
relativen Werten. Um absolute Angaben für EZ und 5 zu erhalten, wäre eine 
Absolutmessung der auf den Brennfleck wirkenden Kraft erforderlich. Es wurde 
noch versucht, diese dadurch zu bestimmen, daß der Brennfleck an einem im Queck- 

i silber schwimmenden Schiffchen mit daran befestigtem, die Quecksilberober- 

fläche durchdringenden Molybdändraht fixiert wurde. Aus der Driftgeschwindig- 
keit dieses Schiffehens und dem gesondert zu bestimmenden Widerstandsgesetz 
sollte dann die Absolutkraft berechnet werden. Leider erwies sich die Hand- 
habung dieser Anordnung wegen der starken Wirkung der Oberflächenspannung 
als so schwierig, daß von den Messungen Abstand genommen werden mußte. 

Herrn Dr. M. Steenbeck danke ich bestens für die Anregung zu dieser Arbeit 
und sein stetes, förderndes Interesse. 


Berlin-Siemensstadt, Siemens-Schuckertwerke A.-G., Gerätewerk Siemensstadt. 


(Bei der Redaktion eingegangen am 16. September 1948.) 
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(Mit 3 Abbildungen) 
Inhaltsübersicht 


Die von Goos und Hänchen!) beobachtete Seitenverschiebung bei der Total- 
reflexion eines begrenzten Lichtbiindels wird im Anschlu8 an eine Arbeit von 
C. Schaefer und R. Pich?) berechnet. Es wird hierbei eine Näherungslösung 
der Maxwellschen Gleichungen verwandt, die der Begrenztheit des Bündels 
Rechnung trägt. Die Verschiebung des Strahles ergibt sich anschaulich aus der 
Betrachtung des Energieströmungsvektors im dünneren Medium. Das Resultat ist 
das gleiche — auch bezüglich der Polarisationsabhängigkeit des Effektes, — das 
K. Artmann®) auf anderem Wege erhielt. Die beiden Theorien geben in gleicher 
Weise die experimentellen Ergebnisse qualitativ wieder. Im Experiment wird 
aber keine Polarisationsabhängigkeit gefunden, während die beiden Theorien 
sie in quantitativer Übereinstimmung behaupten. 


30os und Hänchen!) haben einen interessanten Versuch angegeben, der die 
bei der Totalreflexion im dünneren Medium sich abspielenden Vorgänge nachzu- 
weisen und zu messen gestattet, ohne daß eine Energieabzapfung in diesem Medium 
stattfindet. Ihre Methode ist also eine rückwirkungs- 
freie. Sie ließen ein enges Lichtbündel von rechteckigem 
Querschnitt im Innern einer planparallelen Glasplatte 
an der Grenzfläche Glas—Luft unter einem Winkel 
reflektieren, der bereits im Gebiet der Totalreflexion 
und zwar in unmittelbarer Nähe des Grenzwinkels der 
Totalreflexion lag. (Abb. 1.) Die reflektierende Rück- 


seite der Glasplatte war dabei zur Hälfte versilbert, App. 1. Schematischer 
so daß die eine Hälfte des einfallenden Lichtbündels Strahlverlauf bei Total- 
an Glas (total), die andere Hälfte an Silber reflektiert reflexion 


wurde. 

Sie beobachteten dabei, daß die an Glas total reflektierte Bündelhälfte gegen- 
über der an Silber reflektierten eine seitliche Versetzung erfuhr. 

Goos und Hänchen führen diese Erscheinung (in etwas ungewohnter Weise) 
darauf zurück, daß die Reflexion nicht an der Grenzfläche F (Glas—Luft), sondern 
an einer dahinterliegenden, fiktiven Ebene F’ erfolge, deren gegenseitiger Abstand 
durch die sogenannte „Eindringtiefe‘‘ bestimmt sei. Nach ihrer nur schematisch 


1) F. Goos u. H. Hänchen, Ann. Physik 6. Folge 1, 333 (1947). 
2) C. Schaefer u. R. Pich, Ann. Physik 80, 245 (1937). 
®) K. Artmann, Ann. Physik (6) 2, 87 (1948). 
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gedachten Darstellung (Abb. 1) nimmt der an Silber reflektierte Strahl dabei den 
Weg PQR, der an Glas totalreflektierte den Weg PQST', woraus sich eine seitliche 
Verschiebung vom Betrage D ergibt. Die punktierte Linie im 2. Medium soll dabei | 
das Eindringen des totalreflektierten Lichtes ins 2. Medium andeuten, ohne etwa | 
verbindlich den Strömungsvektor der Energie zu kennzeichnen. Der leitende Ge- 
danke, der zur Anstellung des Versuches führte, war offenbar der, daß bei der Total- 
reflexion das Licht bis zu einer Tiefe in der Größenordnung einer oder mehrerer 
Wellenlängen in das 2. Medium eindringt, was zu der Vorstellung verführt, daß die 
Reflexion nicht an der geometrischen Grenzfläche F, sondern an einer virtuellen 
Ebene F’ erfolgt. Demgegenüber muß man wohl daran festhalten, daß als Ebene 
der Reflexion: immer noch die Grenzfläche zwischen Glas und Luft anzusehen ist, 
da ja in dem Sprung des Brechungsquotienten die Ursache für die reflektierte 
Welle liegt und an dieser Grenzfläche die für die Entstehung der reflektierten Welle 
maßgeblichen Grenzbedingungen zu erfüllen sind. 

Eine einmalige Reflexion ergibt noch keine meßbare Verschiebung. Goos 
und Hänchen ließen daher den Lichtstrahl in der Glasplatte 70mal und mehr 
reflektieren, um dadurch die entsprechende Multiplikation der Verschiebung zu 
erhalten. Da sie das Eindringen des Lichtes in das 2. Medium als Ursache für die 
Verschiebung des totalreflektierten Strahles erkannten, war es naheliegend, die 
Größe der Verschiebung zur ‚‚Eindringtiefe‘‘ in Beziehung zu setzen. Die ‚Ein- 
dringtiefe‘ z, ist dabei etwa diejenige Strecke im 2. Medium senkrecht zur Grenz- 
fläche, auf der die Amplitude der quergedimpften Welle auf den e-ten Betrag abge- 
sunken ist. Nach der elementaren Theorie der Totalreflexion, die mit unendlich 
ausgedehnten Wellen operiert, wird die Dämpfung durch die Exponentialfunktion | 


Vsin® 9, — sin? gg 
e wiedergegeben, wobei z den senkrechten Abstand von der 
Grenzfläche bedeutet, so daß sich für die Eindringtiefe z, ergibt. 


Ay 1 
1 
Goos und Hänchen finden die Seitenverschiebung D experimentell propor: 
tional zur Eindringtiefe: 


n, und n, sind die Brechungsquotienten in Glas (Medium 1) und Luft (Medium 2) 
Y, ist der Einfallswinkel, 9, der Grenzwinkel der Totalreflexion. A, ist die Wellen- 
linge des Lichtes im Glas und k eine Konstante, die Goos und Hänchen aus 
ihren Messungen zu 0,52 bestimmten. Sie haben damit empirisch bereits die rich- 
tige Formel gefunden. Die Verschiebung D ist danach abhängig von dem Einfalls- 
winkel und wächst stark an, sobald sich dieser dem Grenzwinkel der Totalreflexion 
nähert, was experimentell bestätigt werden konnte. Wenn der, Einfallswinkel 
nur noch um 13 Bogenminuten von dem Grenzwinkel der Totalreflexion abweicht, 
beträgt die Verschiebung pro Reflexion etwa 3—4 u. Merkwürdigerweise wurde 
aber der Verschiebungseffekt unabhängig von dem Polarisationszustand des be- 
nutzten Lichtes gefunden, so daß die Formel (2) und der Wert von 0,52 für die 
Konstante k erhalten blieb, gleichgültig ob bei dem Versuch senkrecht oder pe 


zur Einfallsebene polarisiertes Licht verwendet wurde. 
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2 


K. Artmannt) entwickelt in einer späteren Arbeit eine etwas andere Vor- 
stellung von dem beobachteten Verschiebungseffekt. Er faßt das einfallende, seit- 
lich begrenzte Lichtbündel in ähnlicher Weise, wie es wohl zum ersten Male von 
P.Debye°) durchgeführt wurde, als das Resultat der Überlagerung von unendlich 
vielen, unendlich ausgedehnten, ebenen Wellen auf, die in einem engen, räum- 
lichen Winkelbereich mit kontinuierlich variabler Neigung verteilt sind. Jede 
Elementarwelle erfährt bei der Totalreflexion einen Phasensprung gegenüber der 
einfallenden Welle, der — das ist das Entscheidende — von Welle zu Welle ver- 
schieden ist, da die Phasenverschiebung eine Funktion des Einfallswinkels ist. 
Die Zusammensetzung der gegeneinander phasenverschobenen, reflektierten Wellen 
gibt den beobachteten, räumlichen Verschiebungseffekt. 

Die theoretische Bestimmung der von Goos und Hänchen eingeführten 
Konstanten & führt im Gegensatz zum Experiment je nach der Schwingungs- 
richtung des einfallenden, linearpolarisierten Lichtes zu einem verschiedenen Be- 
trag. Schwingt der Vektor der elektrischen bzw. magnetischen Feldstärke senk- 
recht zur Einfallsebene des Lichtes, so ergibt sich für die Seitenverschiebung: 


1 n+? 
= 2 1 
bzw 
Dg. == (=) 1 . . (4) 
n,/) an, Vsin® 7, — sin? y, 
Die Konstante & hat also in den beiden Fallen die verschiedenen Werte = . 
1 
2 
bzw. (*) . u . Nimmt man fiir”, (Glas) den Wert 1,52, für n, (Luft) denWert 1, 


so ergeben sich theoretisch die numerischen Werte 0,21 und 0,48 gegeniiber dem 
experimentell bestimmten Wert von k = 0,52. 

Hierzu ist zu bemerken, daß das Ergebnis der Artmannschen Rechnung in 
einer Arbeit von C. Schaefer und R. Pich®), die weder von Artmann noch von 
Goos und Hänchen zitiert wurde, bereits vollständig enthalten ist, wenn auch 
der Effekt der Seitenverschiebung nicht ausgesprochen wurde. Die Deutung, die 
der Vorgang der Seitenverschiebung in der Schaefer-Pichschen Auffassung 
erfährt, ist noch dazu von besonderer Anschaulichkeit. Das Eindringen in das 
dünnere Medium bei der Totalreflexion eines begrenzten Lichtbündels (Abb. 2) 
vollzieht sich, wie bereits in 2 Arbeiten von Picht?) und Noether®), vor allem 
aber in der Arbeit von Schaefer-Pich klargestellt wurde, in der Weise, daß im 
Randgebiet*R, die Energie in das 2. Medium eindringt, dann unter alimahlichem 
Richtungswechsel in die zur Grenzfläche parallele Richtung umbiegt, diese Rich- 
tung fast über den ganzen Querschnitt des Strahles beibehält, um im Randgebiet 
R, wieder in das 1. Medium zurückzukehren. An der Stelle R, entsteht deshalb 
eine Defizit der reflektierten Energie infolge desVerlustes der eindringenden Energie, 
an der Stelle R, ein Überschuß infolge der zurückflutenden Energie. (Im Gebiete 
R, werden also tatsächlich mehr als 100°, reflektiert.) Diese Energiewanderung 


5) P. Debye, Ann. Physik (5) 30, 755 (1909). ee 
6) C. Schaefer u. R. Pich, a. a. O. 
*) J. Picht, Ann. Physik (5) 8, 433 (1929). ee EIER 


8) F. Noether, Ann. Physik (5) 11, 141 (1931). 
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& ist nun die Ursache für die Verschiebung des Intensitätsschwerpunktes des total- 
reflektierten Lichtes. Das ist aber bereits der beobachtete Effekt der 
Strahlversetzung. Wegen der Diskrepanz zwischen Experiment und. Theorie 


a bezüglich der Polarisationsabhängigkeit des Verschiebungseffektes erscheint eine 
= 2. unabhängige, theoretische Herleitung nicht überflüssig, zumal das Resultat der 
b beiden theoretischen Darstellungen — auch im Punkte der Polarisationsabhängig- 
a keit — genau das gleiche ist. 


Die elementaren Formeln der Totalreflexion gelten bekanntlich nur für eine 
unendlich ausgedehnte Welle. In jedem praktischen Fall hat man es aber mit 
einer seitlich begrenzten Welle zu tun, für die die Maxwellschen Gleichungen 
nicht mit dem üblichen Ansatz zu befriedigen sind. Schaefer und Pich benutzen 
deshalb einen variierten Ansatz, 
der die Maxwellschen Glei- 
chungen genähert löst. 

Eine ebene Welle (Abb. 2) 
falle aus dem Medium 1 
kommend unter dem Einfalls- 
winkel g,, der größer als der 
Grenzwinkel der Totalreflexion 
Y, sein soll, auf die ebene Grenz- 
fläche zweier Medien 1 und 2 
auf (n,> n,). Die ‚‚Breite‘‘ des 
Lichtbündels, quer zum Strahl 
Abb. 2. Die Energieströmung bei Totalreflexion in der Einfallsebene gemessen, 

eines begrenzten Bündels sei B,—f,. Alle Feldgrößen seien 

der theoretischen Einfachheit 

halber von der Koordinate y, die senkrecht zur Zeichenebene steht, unabhängig, 

so daß man sich das Lichtbündel in der y-Richtung als unendlich ausgedehnt 

vorzustellen hat. Die Amplitude des Lichtes soll in der Einfallsebene in der Rich- 

tung ö quer zum Strahl eine Funktion des Ortes sein, derart, daß die Amplitude 

in einem mittleren Bereich ß,—ß, einen konstanten Betrag hat und in den beiden 

Randgebieten R, und R, allmählich auf den Wert 0 absinkt. Statt des Ansatzes 

für die unendlich ausgedehnte Welle (linearpolarisierte Welle, €, senkrecht zur 
Einfallsebene) 


Medium 2 


E,=4,:A(ß) - 


= cos, Ve, € 6 


Ver sing, €, 


ergibt sich fiir unsere (in einer Dimension) seitlich begrenzte Welle ein kompli- 


zierterer Ausdruck mit den unterstrichenen Zusatzgliedern: fee 


cos Ve, €,—ising, Ve; & . ei? 9) 
= sing, Ve E,+ cosg, Vex ei? Dar (16) 
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Hierbei bedeuten: 9, den Einfallswinkel, A (8) die Amplitudenfunktion, die 
den Wert der Amplitude abhängig vom Ort 8 im Lichtbündel in relativem Maße 
angibt. (A (ß) liegt immer zwischen 0 und 1.) 8 mißt den senkrechten Abstand vom 
Mittelstrahl und unterscheidet sich nur durch einen konstanten Faktor von der 
Koordinate 6, die den gleichen Abstand in demselben Maß mißt, in dem x und z 
gemessen werden. 


9 . 
B= -b; b= xcosp, + zsing,. 


A, bezeichnet den Maximalbetrag der Amplitude, der voraussetzungsgemäß im 
Gebiete gilt. ist das übliche Argument einer ebenen Welle = - 7 


\Yaı 
und J die Koordinate in der Fortpflanzungsrichtung: / = 2 sing,— z cosqy. 
ay ist eine Konstante: | x) | < 1, die dafür sorgt, daß die Zusatzglieder klein bleiben. 
a, bezeichnet den Absolutbetrag des Zusatzgliedes. In diesen Zusatzgliedern tritt 


charakteristischerweise immer nur der Differentialquotient a auf. Die Welle 


unterscheidet sich also nur in den ,,Randgebieten‘‘ von der Form der unendlich 


ausgedehnten Welle. Die eingefiihrten Konstanten sind gerade hinreichend, um die 
Grenzbedingungen an der Grenzfläche beider Medien zu befriedigen. Die Zusatz- 
glieder selbst sorgen fiir die Erfiillung der Grenzbedingungen an der Begrenzung 
des Biindels gegen den strahlungsfreien Raum. Wichtig fiir die Giiltigkeit der 
Näherungslösung ist eine einschränkende Bedingung, die zum Ausdruck bringt, 
daß die Änderung der Amplitude,in den Randgebieten hinreichend langsam vor 
sich gehen soll. 


dA 
= 9 
oder wenn man fp = = "X, b einsetzt, 
| | 
<IAI- (12) 


Das heißt: Die Änderung der Amplitude, ausgedrückt durch den 1. Differential- 
quotienten der Amplitudenfunktion nach dem Argument (+) soll auf einer Strecke 


quer zum Strahl von der Größenordnung einer Lichtwellenlänge klein gegen die 
Amplitude selbst sein. Und der 2. Differentialquotient soll wiederum klein gegen 
den 1. sein. Rechnet man mit den neuen Feldgrößen den Poyntingschen Vektor 
an der Grenzfläche Glas—Luft aus (gemittelt über eine Periodendauer), dann er- 
gibt sich nur für die einfallende Welle der gleiche Ausdruck wie bei einer unendlich 


ausgedehnten Welle: 
to+T 


f S dt = _ Ye, E} -At(p). (13) 


to 


E, ist der maximale Amplitudenbetrag für die einfallende Welle. E, ist identisch 
mit dem A, der Formel (8). Der neue Buchstabe E, wird bei Schaefer-Pich 
an dieser Stelle nur zur Unterscheidung von der Amplitude R, der reflektierten 
Welle eingeführt. 


Ve 
7 ; 
(5) | 
(6) 
(7) 
pli- 
(8) 
| 
(9) 
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Während bei einer unendlich ausgedehnten Welle der Mittelwert der ins 2. Me- 
dium einströmendem Energie 0 ist (wohlgemerkt nur der Mittelwert und nicht der 
Momentanwert), trifft dies bei einer begrenzten Welle nur noch.in dem mittleren 
Gebiet M (Abb. 2) zu, in dem die Amplitude des einfallenden Lichtes konstant ist, 
nicht aber in den beiden Randgebieten R, und R,. An dem einen Rand strömt 
Energie nahezu streifend in das 2. Medium ein, um am anderen Rande wieder voll- 
ständig in das 1. Medium zurückzufluten. In dem Zwischengebiet zwischen den 
beiden Rändern fließt die am Rande eingetretene Energie im 2. Medium parallel 
zur Grenzfläche. Die ‚„Totalität‘‘ der Reflexion ist also gewährleistet. Die Rand- 
gebiete erweisen sich als Quelle und Senke für die in das 2. Medium übertretende 
Energie. Daraus ergibt sich dann zwanglos die von Goos und Hänchen gefun- 
dene Seitenversetzung: 

Schaefer und Pich berechnen für den Poyntingschen Vektor der reflek- 
tierten Welle den Ausdruck: 


| = Van (14 


B, ist die bereits erklärte Koordinate # für den reflektierten Strahl, &, entspricht 
bei der reflektierten Welle der Konstante «, bei der einfallenden Welle 


w = ) sin? g,— sin?g, . 
Die reflektierte Energie läßt sich nun formal genau wie die einfallende schreiben, 
wenn an Stelle des Amplitudenfaktors A (f,) ein veränderter Faktor A’ (,) ge- 
schrieben wird, indem dem Umstand Rechnung’getragen wird, daß die reflektierte 


Amplitude in den Randgebieten um einen kleinen Betrag 4A verkleinert bzw. ver- 
größert wird uc der in das 2. Medium einströmenden Energie E 


|S,| = Va E 0 3, Va {A(B,)— _ 
Vex E, {A? (B,) —2A( (B,) AA ( (B,) + 4? A (B, )} 


oder unter Yenssikiäniueng des kleinen 3. Summanden in der RE. 
Klammer 


| 5; Ver Bi {42(B,)— 2A (B,) 4A (,)}- (16) 


Um die Vorstellung zu fixieren, wollen wir z. B. die in Abb. 3 angedeutete 
Intensitätsverteilung quer zum Strahl annehmen. (Es sei aber darauf hingewiesen, 
. daß die gezogenen Folgerungen un- 

Ab) abhängig von dem speziellen Verlauf 
der Intensität sind, wenn nur die ein- 

eifalende Wele schränkende Bedingung (12) erfüllt 
refektene Nele} q ist.) Die Amplitude A (8) soll also 
Ad 4 zwischen ß, und ß, konstanten Betrag 
4A_/V \ haben und in den Randgebieten R, 
L \ bzw. R, allmählich auf 0 fallen. Die 
by Bedingung für die Kleinheit des 
Abb. 3. Intensitätsverteilung im Querschnitt 1, Differentialquotienten dA besagt 


des einfallenden und des totalreflektierten 
Lichtstrahls £ im wesentlichen, daß die „Breite‘ der 
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Randgebiete R, und R, noch groß gegenüber der Wellenlänge sein soll und daß die 
Änderung der Steigung von GE ebenfalls klein sein soll gegenüber der Steigung 


selbst. Wir wollen also etwa im Randgebiet einen nahezu linearen Amplituden- 
anstieg annehmen mit Ausnahme der Stellen ß,, 8, und ß,, ,, an denen mit Rück- 
sicht auf die Stetigkeit des 1. Differentialquotienten von A die Ecken im Ampli- 
tudenverlauf abgerundet werden müssen. 

Die beschriebene Intensitätsverteilung sei die der einfallenden Welle, bzw. die 
der reflektierten Welle in dem fiktiven Fall, daß auch in den Randgebieten eine 
100proz. Reflexion stattfände, ein Fall, der in dem Experiment von Goos und 
Hänchen durch die Reflexion an der Silberschicht nahezu realisiert ist. (Zwar 
beträgt dabei die Reflexion nicht 100%, sondern weniger, sie ist aber — worauf es 
hier ankommt — an allen Stellen der Grenzfläche, also auch in den Randgebieten 
konstant.) Es sei nun fg diejenige Stelle im Randgebiet, an der die Lichtamplitude 
den Schwellenwert erreicht, der von dem Meßinstrument (Auge oder photographische 
Platte) gerade noch wahrgenommen werden kann. 

Die Lichtamplitude an der Stelle Bg sei Ax « Bx ist also diejenige Stelle, die als 
Kante des Schwärzungsbildes erscheint. Wenn nun infolge der an das 2. Medium 
erfolgenden Energieabgabe die reflektierte Amplitude den durch die punktierte 
Kurve angedeuteten Verlauf annimmt, dann wandert der Schwellenwert der 
Amplitude A, nach der Stelle Bx. Für diese Verschiebung Bg — fx gilt aber, 
falls sie nur klein ist, angenähert: 


d4_ AA 
| 
wobei AA die Amplitudenänderung an der Stelle Bx bedeutet. Ersetzt man fk— Px 


nach der Gleichung: x — Bre = (6g—5x) D (wobei D in dem 


gleichen Maß gemessen wird wie die ursprüngliche Koordinate x yz), dann ist D 
identisch mit der von Goos und Hänchen bestimmten Verschiebung. Es gilt also: 


(17) 


AA = "2% (18) 
B, 
Setzt man diesen Wert in den Ausdruck fiir | S ‚| ein, so ergibt sich: 
= 2 dA (B,) 
Und der aii mit (14) . es 
| S,| 82 Ver B,) - A(B,) 
w Ay 4 w 


Da für die ganze Betrachtung immer ein Einfallswinkel in der unmittelbaren Nähe 
des Grenzwinkels der Totalreflexion angenommen wurde, kann sin p, durch n,/n, 
ersetzt werden, so daß sich in Übereinstimmung mit der Artmannschen Formel 
ergibt: 
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Diese Formel gilt, wie die ganze Herleitung, für den Fall, daß die elektrische Feld- 
stärke des linear polarisierten Lichtes senkrecht zur Einfallsebene schwingt. Rech- 
net man die Schaefer-Pichschen Formeln für den komplementären Fall durch, 
daß die magnetische Feldstärke senkrecht zur Einfallsebene schwingt, dann er- 
gibt sich — ebenfalls in Übereinstimmung mit Artmann: 


Artmann versucht die Diskrepanz zwischen Theorie und Experiment beziiglich 
der Polarisationsabhingigkeit des Verschiebungseffektes mit dem Hinweis ver- 
ständlich zu machen, daß in seiner Rechnung noch nicht die Beugungserschei- 
nungen an dem das Bündel begrenzenden Spalt berücksichtigt seien und spricht 
die Hoffnung aus, daß eine strengere Rechnung mit Berücksichtigung der Beugung 
die Theorie wieder in Übereinstimmung mit dem Experiment bringen könnte, 
Verfasser teilt nicht diese Meinung; denn die der Rechnung von Schaefer und 
Pich zugrunde gelegte Welle, bei der der Intensitätsverlauf innerhalb der Flanken 
in weiten Grenzen beliebig gewählt werden kann, hat eine gewisse Ähnlichkeit mit 
einer Welle mit Beugungsfransen, so daß für diese Welle die Beugung gewisser- 
maßen schon berücksichtigt ist, wobei man sich allerdings darüber klar sein muß, daß 
die Frage erst noch untersucht werden muß, ob die einschränkende Bedingung (12), 
die die Art der Intensitätsverteilung in den Flanken reguliert, in Übereinstimmung 
gebracht werden kann mit dem Intensitätsverlauf, wie er sich an den Begrenzungen 
eines Lichtbündels in konkreten Beugungsfällen ergibt. Zum mindesten erscheint 
es näherliegend, die augenblicklich noch bestehende Unstimmigkeit zwischen Theorie 
und Experiment in der Schwierigkeit der Messungen begründet zu sehen und zu 
hoffen, daß sich diese bei noch weiterer Verfeinerung der Meßmethodik wird be- 
seitigen lassen. 


Köln, II. Physikalisches Institut der Universität. 


(Bei der Redaktion eingegangen am 3. September 1948.) 
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